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METHODS OF GEOMETRY OF DIFFERENTIAL 

EQUATIONS IN ANALYSIS OF THE INTEGRABLE 

FIELD THEORY MODELS 

ARTHEMY V. KISELEV 

Abstract. In this paper, we investigate the algebraic and geo- 
metric properties of the hyperboUc Toda equations Uxy = exp(if u) 
associated with nondegenerate symmetrizable matrices K. A hi- 
erarchy of analogs to the potential modified Korteweg-de Vries 
equation ut — u^xx + u^ is constructed, and its relation with the 
hierarchy for the Korteweg-de Vries equation Tt = Txxx-'rTT^ is es- 
tablished. Group-theoretic structures for the dispersionless (2-1-1)- 
dimensional Toda equation Uxy = exp(— Uzz) are obtained. Geo- 
metric properties of the multi-component nonlinear Schrodinger 
equation type systems ^t = i'^xx + */(|^|) ^ (multi-sohton com- 
plexes) are described. 
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'You boil it in sawdust: you salt it in glue: 
You condense it with locusts and tape: 
Still keeping one principal object in view — 
To preserve its symmetrical shape.' 

The Hunting of the Snark, Lewis Carroll. 

The Toda equation (jHSl) and the Toda type equations associated with 
the semisimple Lie algebras ([10]) play an important role in the mod- 
els of modern conformal field theory. We recall that the Toda equa- 
tions appear in gravity theory (J2l E]), in the Yang-Mills field the- 
ory ([HZ!), in differential geometry ([23 E]), and in the classifica- 
tion problems for nonlinear partial differential equations ( |l()()j ). Also, 
the Toda equations are related with the integrable dynamical systems 
([22]), the Frobenius manifolds, and the associative algebra structures 
([in]). The detailed study of the following systems that appear in the 
above-mentioned theories is reduced to the study of the Toda equa- 
tions: these systems are the antiselfdual vacuum Einstein equations, 
the Yang-Mills equations, the Gauss-Mainardi-Codazzi equations for 
complex curves in the Kahler manifolds, the dynamical equations for 
the Laplace invariants of differential equations, the Korteweg-de Vries 
equation, and the WDVV (Witten-Dijkraaf-H. Verlinde-E. Verlinde) 
equation, respectively. 

The algebraic approach towards the description of properties of the 
two-dimensional hyperbolic Toda equations 

u^y = exp{Ku) 

was developed in the papers by A. N. Leznov and M. V. Saveliev f|66j). 
V. G. Drinfel'd and V. V. Sokolov ([22]), B. A. Dubrovin (TB]) et al. 
In these papers, the Toda equations are interpreted as the flat con- 
nection equations on semisimple complex Lie algebras (or the Kac- 
Moody algebras) such that K is the corresponding Cartan matrix. 
Then the Toda equations are said to be associated with the Lie al- 
gebras (respectively, the Kac-Moody algebras). These equations are 
exactly integrable ([HE])- In the fundamental paper |E2], the integrable 
Drinfel'd-Sokolov hierarchies were assigned to the Toda equations asso- 
ciated with the Kac-Moody algebras. The Drinfel'd-Sokolov equations 



SYMMETRIES AND CONSERVATION LAWS IN INTEGRABLE MODELS 3 

are the bi-Hamiltonian Korteweg-de Vries type equations. Still, the 
Lie-algebraic approach does not provide the exhaustive description of 
geometric properties of the Toda equations themselves. Indeed, the 
structure of generators of the Noether symmetries Lie subalgebra for 
the Toda equation, the existence of recursion operators for its sym- 
metry algebra, or the relation between conservation laws for the Toda 
equation and the Hamiltonian structures for the Korteweg-de Vries 
equation was not revealed. In particular, until recently it remained 
unnoticed that all these properties of the Toda equation are preserved 
for the general case of the Toda equations u^y = exp{Ku) associated 
with nondegenerate symmetrizable matrices K, which are not assumed 
to be the Cartan matrices. 

The homological methods developed by I. S.Krasil'shchik, V. V. Ly- 
chagin, A. M. Vinogradov ([IDlinSinniinil), and their scientific school 
are a powerful tool in the analysis of algebraic and geometric struc- 
tures for differential equations. Owing to essential achievements of this 
theory, it was quite natural to apply the homological methods to the 
study of the Toda equation and related systems. 

In this article, we carry out the detailed investigation of the geomet- 
ric properties of the hyperbolic Toda equations. Also, we construct 
new Hamiltonian evolution systems associated with the Toda equa- 
tions and establish their nontrivial links with other important classes 
of the mathematical physics equations such as the Korteweg-de Vries 
equation. 
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Introduction 

The discrete one-dimensional Toda lattice ( 

g„ = exp(g„_i -g„) -exp(g„ -g„+i), n e Z, qn = qn{r) (1) 

is a model nonlinear integrable equation in mathematical physics. The 
equation 

e'^ ■ Qrr = exp [q{z - e) - q{z)) - exp [q{z) - q{z + e)) (2) 

is a continuous analog of Eq. ((T)). Here £ > and q = qij.z). Equa- 
tion is then obtained by using the relation qnij) = qij^ne). The 
equation 

Mt-t- = i-D^ o exp(-u), u = u{t, z), 

is the limit for the family of equations © as e ^ +0. Here D^ is the 
total derivative with respect to z and ^^(r) = g„_i — g„ = u{T,ne). In 
the paper ^^1 Dubrovin proved that Eq. (0) is reduced to the nonlinear 
Schrodinger equation ii e = i. In this paper, we consider the multi- 
component analogs (multi-soliton complexes, see P) 

^t = i^^, + ifm)'^. (3) 

of the nonlinear Schrodinger equation. Here \l/ is the r/i-component 
vector, i = ^/^, and / e C~(M). 

The passage from one-dimensional equation (Q) to the two-dimen- 
sional Toda equations 

Uccy = exp{Ku) (4) 

was described in the paper |HE]- By [^1], the acceleration d'^/dr'^ with 
respect to the time t G M is replaced by the d'Alambertian d^/dxdy; we 
also assume (x, y) E C^. The Toda fields u^{x, y) are counted by the 
superscript j G [1, r]. Suppose (cJj) is the rood system of a semisimple 
complex Lie algebra q of rank r and 

2(ai, aj] 



K 






l"jl 



is its Cartan matrix. In [001, equations (JH) were supposed to be associ- 
ated with the Cartan matrices by default. In the sequel, we extend the 
existing picture and analyse the properties of the Toda equations (JU 
such that i^ is a nondegenerate symmetrizable matrix (see precise def- 
initions below). 

Let the matrix K in Eq. (0)) be the Cartan matrix of the Lie algebra 
of type Ar. Then we set 

u\x,y) = u{x,y,z)\^^.^ 

for each 1 < j < r. The continuous limit of Eq. (^ as r — > oo and 
e: — > -|-0 is the dispersionless Toda equation 

u^y = exp(-M2^). (5) 
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This equation appeares in many problems of mathematical physics, 
for example, in the gravity theory ([H], see also the paper ^87j and 
references therein). 

In this paper, we investigate algebraic and geometric properties of 
the Toda equations (jlHHI) and the nonlinear Schrodinger equation (jH)). 
Also, we describe the relation between the Toda equation (j3)) and the 
hierarchies of the Korteweg-de Vries equations and (fTHjl (see below). 
We apply modern cohomological methods and algorights in our study 
of the geometric properties of all these equations. By using the alge- 
braic approach, we reject the idea of a decisive break-through and the 
inevitably immense calculations in coordinates. We operate with the 
notions of homological algebra in the category of infinite prolongations 
of differential equations instead. 

The results exposed in this article are also found in the sequence of 
pubhcations |Ilj-|Hn]. 

The paper is organized as follows. In the introduction, we fix nota- 
tion, formulate necessary definitions, and give a brief list of the final 
results. 

In Chapter 1, we describe important properties of the symmetry 
algebra as well as conservation laws for the Toda equations. In Sec. ^ 
we review some known properties of the hyperbolic Toda equations 
associated with the semisimple Lie algebras. Also, we define the Toda 
equations associated with nondegenerate symmetrizable matrices. In 
what follows, we study these equations that require less restrictions 
upon the matrix K. In Sec. |2] within Chapter 1, we obtain the Noether 
symmetries of the Lagrangian for the Toda equations. In Sec. 01 we 
construct a continuum of recursion operators for the Toda equation. 

In Chapter 2, we construct the commutative hierarchy 21 of analogs 
for the potential modified Kortweg-de Vries equation. This hierarchy 
is identified with a commutative subalgebra of the Noether symmetries 
algebra for the Toda equations. Also, some aspects of the Hamiltonian 
formalism for the hyperbolic Toda equations themselves are discussed. 
Then we relate the symmetry hierarchy 21 with the higher Korteweg-de 
Vries equations. 

In Chapter 3, we illustrate the methods of geometry of partial differ- 
ential equations. Namely, we investigate the properties of the disper- 
sionless Toda equation and the multi-component analog of the nonlin- 
ear Schrodinger equation, which is related with the former equation. 

In Chapter 4, we consider Backlund transformations for the Toda 
equation associated with the algebra 5 [2(C) (the Liouville equation) 
and related systems. Also, we obtain one-parametric families of these 
transformations. Examples of integrating Backlund transformations in 
nonlocal variables are given in Chapter 4, zero-curvature representa- 
tions are pointed out, and the relations between all these structures 
are analysed. 
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In Appendix El several geometrical methods of solving the boundary 
problems for equations of the mathematical physics are listed. 

Basic definitions and notation. First let us introduce some notions 
from the geometry of partial differential equations. We follow [U HUl 
Ea Ea ill and the articles UHlllZj. 

Differential equations and their symmetries. Let S = {F"(a;, u, p) = 
0, 1 < a < r} be a partial differential equation of order k imposed by 
the smooth functions F°' on n independent variables x = *(x^, . . ., x"), 
r dependent variables u = *(m^, . . ., u'^), and the derivatives 

p = {pi\p>^ = 9I"ImV9(x^)*^ . . . d{xy\ 

a = {ii,..., in}, |cr| = zi H \- in < k}. 

Consider the trivial r-dimensional smooth fibre bundle 

TT : M*- X M" ^ M". 

The space J''{n) of kth jets of the bundle vr is the union |J^ J^, where 
J^ is the set of all equivalence classes [s]^ of sections s in the bundle vr 
such that the sections s are tangent with order k at the point a; G M". 
Define the sequence of smooth fibre bundles 

by the formula 

We denote by vr^ the smooth fibre bundle tt^: J^{ti) -^ ffi" defined 
by 7rfc([s]a,) = X. Also, by J-'ki^!') we denote the ring of smooth (C°°) 
functions on J'=(7r). Finally, we set J'-ooItt) = C~(M"). 

We treat the variables x, u, and p as coordinates in the jet space 
J^(7r). By the Hadamars lemma (|3H|); two sections s, s' G r(7r) are 
equivalent in J^(vr) iff their partial derivatives at the point x coincide 
up to order k. 

A submanifold 

^ = {F" = 0} C J^Ti) 

is a differential equation of order not greater than k imposed on r 
functions of ra independent variables. An equation S C J^i'n:) is regular 
if the mapping 



TTfcl^l^^M 



is a surjection. 



Definition 1. Consider the sections s of the fiber bundle vr such that 
[s]^ = for some point 9 G J^{t^)- Consider the tangent planes to the 
graphs rj of the fcth jets of s. The linear span of all these tangent 
planes is the Cartan plane Cg = Cg at the point 6. The union C of the 
mappings 6 ^-^ Cg hj all 6 G J^{t^) is the Cartan distribution. 
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Definition 2. The subset 

^(') = {^fe+^ = [sr',sGr(vr)| 

I jkis){x) is tangent to S at the point 9k = [s]^ with order > /} 

in J^^\Ti) is the /th prolongation of an equation £ C J''{tt). The 
inverse hmit 

S^ = proj hm ^('^ 

with respect to the projections 

is called the infinite prolongation of the equation S. In what follows, we 
omit the superscript oo if it is clear from the context that the infinite 
prolongation S°° but not the differential equation S is considered. 

The infinite prolongation of the empty equation {0 = 0} ~ J^iT^) is 
called the infinite jets space J°°(7r). The projections 

and 

VToo: J°°(7r)^M" 
are defined by the formulas 

T^oo,kiPoo) = 0k 

and 

respectively, where 

Ooo = {x, Ok eJ\TT)\ken}e J°°(7r). 

The algebra J^{t^) of smooth functions on J°^(7r) is defined by using 
the following procedure. We set 

J^{Ti) = [Jj^k{-n)^ ke{-oo}un. 

k 

The module A*(7r) of differential i-forms on J°°{7r) is defined by the 
relation 

A;>0 

The restriction C^ of the Cartan distribution C onto the equation 
S^ is an n-dimensional Frobenius distribution. This distribution de- 
fines the decomposition of the tangent space to S°^ into the horizontal 
and vertial subspaces. The horizontal component is generated by the 
restrictions of the total derivatives 
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onto £°° . These restrictions are denoted by Z)j. The dual description 
of the Cartan distribution Cs in the language of differential forms is 
also useful. The de Rham differential on 8^ is representable as the 
restriction onto 8°^ of the sum composed by the horizontal differential 

n 

d/j = ^ da;* ® Di 

i=l 

(this is the lifting of the differential on the base of the fibre bundle tt) 
and the Cartan differential 

dc = d - dh. 

Therefore the space A'(£^) of differential /-forms on the equation £°° is 
the direct sum 

A'(^)= QK\£)®aK{£) 
i+j=i 
of the horizontal i-forms A*(£^) and the Cartan j-forms C^A{S). The 
Cartan forms uj^ = dc(M^) constitute a basis in C^A{S). Here -u^ are 
coordinates on S^. 

The horizontal differential dh generates the horizontal de Rham com- 
plex 

^ J^(7r) ^ A\n) ^ . . . ^ A"(7r) ^ 

of the space J°°(7r). The cohomologies of this complex are called 
the horizontal cohomologies and are denoted by H^{7i). We denote 
by H^{£) the cohomologies of the restriction of this complex onto the 
equation £°° . From the definition it follows that the equivalence classes 
\q\ G H'^^^{£) are conservation laws for the equation £. 

Definition 3. (1) The evolutionary derivation is the operator 

where ip^ G C°°{J^{7r)) for some k and D^^ is the composition 
of total derivatives D^ that corresponds to the multiindex a. 
We denote by x the jF(7r)-module r(7r) (<?>jr_^ •^(^)- Also, we 
denote 

x = Hom^(^)(x,A"(7r)). 

(2) Suppose A^ is the nonlinear differential operator determined 
by a function ip G C°°(J'^(7r)). The operator i^ that acts by 
the rule 

Uiv) = 9v(^) 
is called the universal linearization operator of A^. In coordi- 
nates, we have 

T^-D .1 
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Any Lie field X that preserves the Cartan distribution C can be 
decomposed to the sum X = 3^ + Y, where Y E C and B^ is the 
evolutionary vector field. Any infinitesimal transformation of the space 
J^{tt) can be extended up to the Lie field. In coordinates, the lifting 
rules are the following. The field 



X = ^a,A + 5^3,^._„, 



i,j 



is assigned to the field 



'^ dx' ^-^ ^ du^ 

« i 

An example is found in Eq. (J35|) on page 12^1 

Definition 4. A symmetry of an equation E'^ is a vr-vertical vector 
field X such that X preserves the Cartan distribution C^ = IJee^"" ^^- 

Theorem 1 (,10,). Suppose S C J^{t^) is an equation 

{F'=0,...,F^ = 0}, 

such that 

vroo,o(^°°) = J'i-rr). 
Then the symmetry Lie algebra sym£^°° is isomorphic to the Lie algebra 
of solutions to the system of the determining equations 

3^{F) = onS 

or, equivalently, of the equations 

owing to the definition of the linearization operator L Here ^p E ><Loo- 
The bracket 

endowes the algebra sjm.S°° of solutions ip with the Lie algebra struc- 
ture. 

In what follows, we identify the generating sections 
of symmetries 

(we denote u^. = Dcr{u^)) of differential equations with the symmetries 
9^ themselves (see j221lini)- Roughly speaking, the component y?* of a 
generating section ^p measure the velocity u\ of the dependent variable 
u^ along the "integral trajectories" of the field 9^. 
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Definition 5. Suppose (^(m, . . ., u^) is a symmetry of a differential 
equation S, wliere a is a multiindex and 

is the derivative of the dependent variable u. Assume that the symme- 
try if has the flow 

Ar'. u{x,0) t— i> u{x,t). 

This flow is defined on solutions of the equation Ut = ^ and maps 
the solutions s{x) = u{x, t)|^^q of equation S to solutions of the same 
equation at points r > 0. A solution s{x) of the equation £ is called 
(p-invariant if s{x) is a stationary solution of the evolution equation 

Ur = (piu,. . .,Ua). 

Therefore, the search of the y^-invariant solutions for a given equation 
£ = {F = 0} is reduced to solving the system {F = 0, ip = 0}. 

No we specify an important class of Hamiltonian equations within 
the set of evolution equations ut = f(t, x, u, p). Here we formulate 
the definitions of the Poisson bracket and the Hamiltonian operators 
class (see [1011211112] for further detailes). 

Definition 6. Let A G CDiff(x(7r), x(7r)) be an (r xr)-matrix operator 
in total derivatives: 

A=\\A'^\, A'^ = ApD^, 

and let £i, £2 ^ H"{7!') be two Lagrangians. By definition, set the 
Poisson bracket (the variational bracket) on H"'{tt) by the formula 

{c„ c,u = (E(/:i), A(E(A))> = [E ^ ■ ^" (if) H ' (^) 

where (■, ■) is the natural coupling k{t\) x x(7r) -^ W\'r^ and [■] denotes 
the equivalence class of differential forms. 

Definition 7. Suppose the operator A is subject to the assumptions 
of the previous definition. Then the operator A is called Hamiltonian 
if the Poisson bracket defined in Eq. ® endows H'"'{7t) with the Lie 
algebra structure over M such that the followin relations hold 

{£i,£2}a + {A,A}a = 0, (7a) 

{{£1, AJA, AJA + {{A, CsJA, Ci}a + {{£3, Ci}a, C2JA = 0. (7b) 

The bracket {•, -Ja is the Hamiltonian structure. 

Two Hamiltionian operators Ai, A2 are compatible if XAi + fiA2 is 
a Hamiltonian operator again for any A, ;U G M. 
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Condition ^ is satisfied iSA + A* = 0. Refer [10112311121112 for the 
suitable criteria^ which check whether Jacobi's identity ()7b|) holds for 
a given operator A. For example, any skew-symmetric C-differential 
operator A with constant coefficients is Hamiltonian. 

Definition 8. The evolution equation 

ut = A{E^{n)) (8) 

is a Hamiltonian evolution equation assigned to the Hamiltonian Ti G 
H"'{it) by the Hamiltonian operator A. 

Example 1 ([10]). The Korteweg-de Vries equation 

Tt = -(^T^xx + 3TTx, (3 = const, (9) 

is hamiltonian with respect to the pair 

Bi = Dx, B2 = -l3Dl + 2T- Dx + T, (10) 

of the compatible Hamiltonian operators. Indeed, we have 

Tt = B^o EtHI^T^ + \T^) dx) = 4 o Er(iTMx). 

Conservation laws. In this subsection, we formulate important defini- 
tions and statements concerning conservation laws for differential equa- 
tions. Also, we recall the relation between symmetries and conservation 
laws for the Euler equations. 

Definition 9. A conservation law 
[v] e H''-\£) = {uje A"-i(^) I d,(^) = 0}/ 

/ {u e A^-\S) I Lu = d,7, 7 e A"-2(^)} 

for an equation S is an equivalence class of horizontal (n — l)-forms 
rj G A"^^{S) closed on £, 

dhT] = V(F) dx. 

Here 

n 

dh = Y^ dx' (g) Di 

i=l 

is the restriction of the horizontal differential dh on S, D^ is the total 
derivative with respect to x*, V is an operator in total derivatives, and 
the equivalence is the factorization by exact forms dh'y, 7 G A"~^(£^). 
Representatives rj of the equivalence classes [rj] G H'^~^{£) are called 
conserved currents for the equation £. 

The of notion of the Poisson bi-vectors is a useful intrunient in the Hamihonian 
operators theory (see |2(ill42| ). For example, the first Hamiltonian structure for the 
Korteweg-de Vries equation is 1 A -D^; and the second structure is 1 A B2, see 
Eq. H10|l . A bi-vector A is Poisson if and only if A satisfies the equation lA, A] — 0, 
where |-, •] is the Schouten bracket. A pait of the Poisson bi-vectors Ai and A2 
is compatible if |Ai,^2l — 0. In the papers [H21E3|, the equation |A,A] — was 
considered in a more general situation such that A is not necessarily a bi-vector. 
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Example 2 (|Sni)- The horizontal 2-form 

r] = Uxz exp(-M2;,) dx Ady + {\u].^ - u^x) dx A dz 
is a conserved current for the dispersionless Toda equation 

Uxy = exp(-M^^). 
Indeed, the relation 

Dy{\'^lz - Uxx) = Dz{uxz exp{-Uzz)) (H) 

holds on that equation. 

Definition 10. A regular equation S = {F = 0} is called i-normal if 
the condition 

implies V = 0. 

Proposition 2 (JD])- Letn be the number of independent variables x^ , 
. . ., x^ and let £ be a regular equation. Consider a coordinate neigh- 
bouhood Q{9°°) C £°° of a point 9°° G £°° and suppose that there is 
a set {v}of the internal coordinates v on Q such that the total deriva- 
tives Di{v) can be expressed in terms of these coordinates {v} for any 
i, 1 < i < n. Then the equation £ is i-normal. 



If £ is an £- normal equation, then the compatibility complex ([H 
inn]) for the equation £ has length 2, and the equation satisfies the 
assumptions of the '2-line' theorem. 

Remark 1. The Maxwell, the Yang-Mills, and the Einstein equations 
are not ^-normal since there is a nontrivial dependence between the 
equations that is provided by the gauge symmetry pseudogroup. 

Now we need a test that checks whether a given equation is ^-normal. 
We have 

Example 3 (fOlinSI)- The evolution equations are £- normal. 

Conservation laws [q] for the ^-normal equations are described by 
their generating sections ip = V*(l) G x. We see that 

d,r/=(V(F),l) = (F,V*(l)) + d,7. 

Then the coupling 

(V*(l),F) = d,(r/-7) (12) 

is an exact horizontal n-form. Suppose t] is trivial and therefore V = 0. 
Then, obviously, then ip = 0. 

Lemma 3 (IHII)- Let £ = {F = 0} be an i-normal equation and 
assume that H''{£) C if"(f) {e.g., H''{£) = 0). Suppose that the 
generating section ip of a current rj is zero. Then the conserved current 
f] is trivial. 
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In what follows, we shall use the following remarkable theorem many 
times. 

Theorem 4 ([94j). Let S = {F = 0} be an ^-normal equation in the 
fibre bundle tt: M™ x R" -^ M". Then the generating sections ip of 
conservation laws [rj] satisfy the equation 

tpi^v) = 0, (13) 

where i*p is the operator formally adjoint to ip and i*p is the restriction 
of i*p onto the equation S. 

Proof. Apply the Euler operator E to both sides of Eq. (fT^ . Thence 
we obtain 

= E((^, F)) = £^^,^^(1) = iU^/j) + i;iF) = (14) 

by the Leibnitz rule. Now restrict Eq. ()14|) onto the prolongation of F = 
and obtain the determining equation (fT^ imposed on the generating 
sections and satisfied on the initial equation S. D 

The generating sections ipri are often named gradients of the con- 
servation laws [rj] (see j221) by the following reason. The generating 
sections for evolution equations lie in the image of the Euler operator 
(that is, the "gradient" 6/6u) applied to the corresponding conserved 
density. 

Lemma 5 ([HHl IHS!)- Let £ = {ut = f{t, x, u, ui, . . .)} be an evolution 
equation. Assume that 

n 

T] = riodx+ ^(-l)*~^?7i dt Adx^ A-- ■ Adx' A-- ■ A dx'^ 

i=l 

is a conserved current for S: dh{f]) = 0, that is, 

A(r/o) + 5^5.«(r/i)=0. 

i 

Then its generating function is 

V^, = E(r^o)=C(l). 

The proof of Lemma El is straightforward. 

The Euler-Lagrange equation Se-l assigned to a Lagrangian C G 
i7"(7r) is 

Se-l = {G = £„(£) = 0}. (15) 

We recall that the Helmholz condition ic = ^h i^ valid for the image 
G of the Euler operator E. 

The correlation between conservation laws [77], their generating sec- 
tions ipr^, and the Noether symmetries ifc ^ sym C of an Euler-Lagrange 
equation £ = {E(£ = 0} is guided by the following version of the 
Noether theorem. 
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Theorem 6 ("S]). Let £ = {E(£) = 0} be the Euler- Lagrange equation 
for a Lagrangian L. Then the evolutionary derivation 9^ is a Noether 
symmetry of the Lagrangian C: 

Q^iC) = 

if and only if (f is the generating section of a conservation law [rj]: 
dfiT] = on S. 

Proof. Suppose ip is the generating section of a conservation laws [rj] 
for the equation £: 

d,(r/) = (l,n(F)) = (n*(l),F) + d,(7), 

where the couphng (-, ■) takes values in the space A"(7r) of horizontal 
n-forms. Thence, {ip, F) = d/i(?7 — 7) is an exact horizontal n-form. We 
have F = E(£) = £^(1) by the initial assumption, whence we easily 
obtain 

(y., F) = {^, £1(1)) = (4(y.), 1) + d,P = (B^(L), 1) + d,p. 

Therefore, (f is a Noether symmetry. The second implication in the 
theorem's statement, the sufficiency, is proved by inverting the reason- 
ings, n 

Coverings. Let £^ be a differential equation in the fibre bundle M" x 
M™- -^ M" and S°° be the infinite prolongation of S. The n-dimensional 
Cartan plane Cg C Tq{£'^) is defined at each point 9 G £°° . The Cartan 
distribution C^ on £°° is Frobenius: 

[Ce,Ce\(lCe. 

In local coordinates, Cs is defined by n vector fields Di, . . . , Z)„, where 
Di is the restriction of the total derivative with respect to the iih 
independent variable x* onto £°^ . 

Definition 11 (POl)- The equation £'^ with the n-dimensional Car- 
tan distribution C and the regular map r are called a covering over the 
equation £°° if at any point 9 G £°° the tangent map r^, e is an isomor- 
phim of the plane Cq onto the Cartan plane Cr{e) on the equation £°° 
at the point t{9). The equation £ is called the covering equation. The 
dimension of the fibre in the bundle r is the dimension of the covering. 

In coordinates, the structure of a covering is realized in the following 
way. Locally, we can realize the manifold £ and the mapping r: £ ^ 
£°° as the direct product £°° x W (W C R^ is an open subset, < 
N < 00) and the natural projection £°° xW ^ £°°, respectively. Then 
we describe the distribution C on £ by the vector fields 
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where Xij G C°^{S) are the coefficients of the r- vertical fields on S 
and Si, . . . , and sn are the Cartesian coordinates in M^. Then, the 
Frobenius condition [C, C] C C of integrability for the distribution C is 
equivalent to the system of relations 

[Di,Dj] = 0, z,j = l,...,n, 

which are equivalent to the equations 

Di{Xjk) = Dj{Xik) 

that hold for alH, j = 1, . . . , n, < A; < A^. 

The coordinates Sj are called nonlocal variables. In coordinates x*, 
M^, and Sj the rules Di{sj) = Xij to differentiate the nonlocal variables 
Sj as well as the initial equation S°° define the covering equation S. 

Example 4. Again, consider the Korteweg-de Vries equation Q. Ex- 
tend the set of variables t, x, Tj = Di(T) with the "nonlocality" 
s = J T dx: 

Sx =T, St = —pSxxx + 2^x- (16) 

We see that the equation that covers © is the potential Korteweg-de 
Vries equation. 

A symmetry of the covering equation £ is a nonlocal symmetry of the 
equation £°°. Suppose that the field X is a symmetry of an equation 
S°° and t: S -^ S°° is a covering. Two radically different cases are 
possible: 

(1) the symmetry X of the equation £^°° can be extended up to a 
symmetry X of the covering equation S, and 

(2) the converse: any lifting of the symmetry X is not a symmetry 
of the covering equation. 

In the second case, the field X generates the one-parametric family of 
equations St that cover the equation S°° for any t. 

Recursion operators. In this subsection, we describe an explicit method 
(in2l) of constructing recursion operators for differential equations. The 
Cartan generating forms approach developed by I. S. Krasil'shchik is 
the key-point of our reasonings. 

Let £ = {F = 0} be a determined system of differential equations, 
that is, the number m of dependent variables in this equation is equal 
to the number r of equations. Suppose (p G sym£^°° is a symmetry of 
this equation and consider an r-component column 

such that its components 

are the Cartan 1-forms that vanish on the Cartan distribution Cs on 
£°° and the coefficients at these forms are functions on S°^. Obviously, 
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the component-wise action 3^ _i tu = </)' is again an element of the 
module x of evolutionary derivations. The condition for Lp' to be a 
symmetry of the initial equation E is 

4"M = 0, (17) 

where ip is the restriction of the linearization operator onto C^k.{E\ 
Now we formulate the rule that assigns differential operators R in 
total derivatives to the columns (the generating forms) uj. Suppose the 
components tu* of a generating form u = *(ci;^, . . ., u"^) are 



iu' 



Then we obtain 

We see that the component- wise action of the Cartan forms cj* is equiv- 
alent to the action of the operators 






on the components of the initial generating section ip = ^{<f^, . . ., f^). 
Hence we deduce that the formal differential recursion operator R is 
the (r X r)-matrix, 

R=\\J2 <'^- 

II ' 'a 

Still, we observe that "usually" equation (fT7|) admits only the trivial 
solution 

t / \ T \ 

^0 = [^0, ■ ■ ■ , ^0) 
that corresponds to the identity recursion operator 

id: ip \-^ if' = if. 

The explanation is that the recursion operators for known equations 
of mathematical physics usually contain summands that involve the 
negative powers D~} of the total derivatives D^- 

Example 5. The recursion operator -RkcIV = B2oB^^ for the Korteweg- 
de Vries equation is equal to 

i?Kdv = -P Dl + 2T + T,- D-\ (18a) 

The recursion operator -RpKdv for the potential Korteweg-de Vries 
equation (fTH|) is 

i?pKdv = -P Dl + 2s, - D-^ o s,,. (18b) 
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From the geometric viewpoint, the situation is hke follows. We ob- 
tain nontrivial solutions of Eq. ()17|) by extending the set of local co- 
ordinates {x^, . . ., x", M^) by the "nonlocal" variables s*. These non- 
localities are defined by the compatible rules of their derivation with 
respect to the independent variables. Also, to each s* we assign the 
Cartan form 

n 

dc(sO = ds^-^4,dx^ 

i=i 
which now belongs to some larger distribution C^A{£). 

Consider the restriction ip of the linearization operator ip onto the 
extended set of the Cartan forms. Substitute the new linearization 
operator in equation (fTTj) and assume that a nontrivial solution uj ap- 
peared. Now we see that the generating form uj is coupled with some 
prolongation 

d 

i 

of the initial evolutionary field B^. Possibly, this prolongation involves 
a larger (ingeneral, an infinite) set of nonlocal variables, see the pa- 
per 113]. The result 

of the action of uj satisfies the equation 

If{v) = 0. 

These sections '{)'{x, u^, s\ . . .) are called shadows of nonlocal sym- 
metries. In general, there are shadows that cannot be extended up 
to a true nonlocal symmetry by using a prescribed set of the nonlocal 
variables s\ An example is found in Proposition 1581 on pagel^M 

Example 6. Consider the following one-dimensional covering over the 
potential Korteweg-de Vries equation ([THjl . Introduce the nonlocal 
variable ( such that 

Cx = -Isl Ct = PsiS3 - iPsl - sl (19) 

where Sj = -D^(s) and we similarly denote Q = -D^(C)- One easily 
checks that the covering equation is the Krichever-Novikov type equa- 
tion 

Ct = -/9 Cs + iPCkr' - 2V2{-Cif'. (20) 

Then the equation 

^?iidv(^) = 

admits a nontrivial solution in the extended set of variables t, x, Sj, 
and (. Indeed, the Cartan generating form of the recursion operator 
for Eq. (Uni) is 

c^pKdv = -Pdc{s2) + Sl dc(s) - dc(C). (21) 
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Now we describe the relation between recursion operators R in to- 
tal derivatives, which are convenient in applications, and solutions u 
solutions of the determining equation 

It is sufficient to solve this problem for the Cartan forms dc(s*). Now it 
is clear that the problem is reduced to calculation of the linearization 
^gi with respect to the nonlocal variable s*. To do this, we use the 
following lemma. 

Lemma 7. Let s G ^(vr) he a function] fix a superscript i & [1, n] of 
the independent variable x\ Then the relation 

^A(s) = Ao4 (22) 

holds. 
Proof. Suppose (p E x, then we obtain 

^D,(s){v) = 'dipO Di{s) = A o B<p(s) = A o 4(V5), 
whence follows relation ()22|) . D 

Lemma [7| defines the rule 

4. = d;' o 4^(^,) (23) 

of calculating the linearization of a nonlocality s^ . Here i G [1, n] is 
arbitrary. Therefore, to each generating form uj = *(ci;^, . . ., uj"^), where 

we assign the matrix operator 

where the index i enumerates the rows and 1 < A; < n is arbitrary. 

Example 7. Recursion operator (jlSbjl corresponds to the generating 
Cartan form (PT|) . Operator p8a|) corresponds to the Cartan form 

^Kdv = -/?dc(T2) + 2Tdc(T) + ridc(s). 

Bdcklund transformations. The notion of a covering is very useful in 
description of Backlund transformations between differential equations. 

Definition 12 (dUl)- Let Si C J'^'^iii), i = 1,2, be two differential 
equations and Ti: S ^ £°° be coverings with the same total space S. 
Then the diagram 

B{8,n,8i) = {8, ^-^i— 8 -^^^ 82] (24) 

is called a Bdcklund transformation B{8,Ti,8i) between the equations 
8i. Diagram (j^ is called a Bdcklund autotransformation if 8^ = 

coo coo 
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In what follows, any Backlund transformation between equations £ 
and £' is a system of differential relations imposed on the unknown 
functions u and u' such that the following property holds. Suppose the 
function -u is a solution to the equation £ and both functions u and u' 
satisfy these relations, then the function u' is a solution to the equation 
£' , and vice versa. 

Example 8. Backlund autotransformation for the two-dimensional La- 
place equation A2f = is defined by the relations 

Remark 2. Let Tj : £j — ^ £°°, j = 1, 2, be two coverings and fi: £i -^ £2 
be a diffeomorphism that maps the Cartan distribution C^jD(£^i) into 
Cr2^{£2)- Then the diagram B{£i, ri, T2 o /x, £j) is also a Backlund 
transformation between the equations £j and these coverings ri and 
r2 o /i are called equivalent. 

Remark 3. Let t: £ ^ £°° be a covering and /i: £^ ^ £^ be a nontrivial 
diffeomorphism of manifolds that preserves the Cartan distribution, 
e.g., a discrete symmetry that cannot be restricted onto £°°. Then the 
diagram 

£ ^-^- £ -^^ £ -^-^ £ (25) 

is also Backlund autotransformation for £. In Chapter 4, we apply 
this contruction to the study of Backlund autotransformation for the 
Liouville equation Uxy = exp(2M). 

Now we give a brief exposition of the mail results obtained in the 
present article. 

Main results. In Chapter 1, we consider a problem standard in the 
geometry of differential equations. Namely, we describe relations be- 
tween symmetries, conservation laws, the Noether symmetries, and re- 
cursion operators for the Toda equations. The structures obtained in 
this chapter are essentially used in what follows. 

Let K = \\kij, 1 < i, j < r|| be a nondegenerate (r x r)-matrix and 
K~^ = \\k'^m be its inverse. Suppose there is the set a of real numbers 
{cj 7^ 0, 1 < ? < r} such that the symmetry condition Kij = Kji holds 
fot the matrix k, = \\K,ij\\ whose elements are 

Then we say that the matrix K is symmetrizahle. 

The hyperbolic Toda equations associated with a nondegenerate sym- 
metrizable (r x r)-matrix K are 

r 

^Toda =[F, = v^^y- exp (^ hju') =0, 1 < ^ < r } . (26) 
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In particular, suppose g is a semisimple Lie algebra of rank r, let {a,, 
1 < z < r} be the system of its simple roots, and K = \\kij = 2(aj, aj) ■ 
|cti|~^; 1 < "^5 J < ''"II be the Cartan matrix of the algebra g. Tne we set 



(Ui'i — ^-^2 



-2 



The Toda equation (PUJ) associated with this matrix K is 
said to be ([66]) associated with the Lie algebra q. 

The Toda equations ()26|1 are the Euler-Lagrange equations with re- 
spect to the action functional 



-^Toda = / -^Toda dx A dy 

with the density 

r r r 

i,j=l i=l j=l 

The canonical Hamiltonian structure is known for the Toda equations 
(see m- 

Suppose i^ is a nondegenerate symmetrizable matrix. Then the Toda 
equation (j^ associated with K admits at least one integral ( |l()()j ). 
that is, an expession dependent of u^ whose total derivative Dy vanishes 
on this equation. Namely, we have 

r r 

^=\Y. '^^i'<< - 5^ «* ■ ^L e ker Dy. (27) 

By definition, put 

T,^Di{T). 

The differential consequences Tj to the functional T generate the sub- 
space T C keri)y within the kernel of the total derivative Dy. Indeed, 
any smooth function Q supplies the functional 

Q(x, T) = Q{x, T,Tu..., T^) e ker Dy. 

We say that the nondegenerate symmetrizable matrix K is generic if 
the integral T in Eq. ^n\ is a unique solution to the equation Dy(T) = 
on the corresponding Toda equation (|^. 

Still, by a special choice of the matrix K one can achieve the situ- 
ation such that the functional T will not be a unique integral for the 
Toda equation pH|l . By jHOj, a necessary and sufficient condition for r 
nontrivial independent solutions Q^ of the equation 

Dy{n') = 

to exist is that K be the Cartan matrix of a semisimple Lie algebra q. 
The Toda equations associated with q are exactly integrable ([HSIIZZ!)- 
From now on, by il C ker Dy we denote the subspace in ker Dy, which is 
differentially generated by all solutions f2* to the equation DyiVf) = 0. 
The total number of these solutions is denoted hy q: 1 < i < q < r. 
We always assume Q^ = T. 



SYMMETRIES AND CONSERVATION LAWS IN INTEGRABLE MODELS 21 
^ . . r 

By A = I A* I, where A* = ^ k"^^ , we denote the vector of the con- 

i=i 
formal dimensions of the Toda fields exp(n) = *(exp(M^), . . ., exp{u^)). 
The transformation 

X H-> X{x), 

y ^ yiy), (28) 

u\x, y)^u' = u\X, y) + A^ log X'{x)y{y) 

is a finite conformal symmetry of the Toda equation £^Toda- The generat- 
ing section (p of the infinitesimal form of conformal transformation ()28|) 
is 

(^ = n(0(x)), 

where the vector-valued, first-order differential operator D is 

n = u., + A-D^ (29) 

and is an arbitrary smooth function. The structure of the Lie algebra 
sym £^Toda generators is described in the paper jTS] . We have 

(1) Suppose i^ is a generic nondegenerate (r x r)-matrix and the 
operator D is assigned to this matrix by using formula (jSHI)- 
Then any symmetry of the Toda equation (j^H|l is 

<f = n{<p{x,n)), (30) 

where (f) is an arbitrary smooth function that depends on a set 
of the integrals ft C ker Dy . 

(2) Suppose the matrix K is subject to additional constraints such 
that the Toda equation (J26|l admits q independent integrals 
n* G ker Dy, where 1 < z < g < r and 

n = {n] = Dim} 3 T. 

We also assume that the first-order operator D is the (r x 
g)-matrix and satisfies some additional restrictions (see [7H]). 
Then symmetries (p of the Toda equation are of the form ()3()|1 
with respect to the new notation. 

The Jacobi bracket on the symmetries ip G symi^xoda induces a 
bracket on the arguments of the operator D. Suppose 

</.i = n(0i(x,T)) 



and 



Then we have 



y.2 = n(02(x,T)). 



Wl,V2} = □(</){l,2}), 

where 

0(1,2} = 9^51 (02) - 9^^2(01) + Dx{4>l)4>2 -4>lDxi<p2) 
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Here 0{i,2} = (p{i,2}{x, T), since the evolution T^ of integral (jTTj) along 
any symmetry (p = d ((/>) is equal to 

f^ = -[3Dl{<P) +TDM + D,.(T ■ 0). (31) 

Here and below we use the notation 

r 
i=l 

Consider the operator applied to the function in the right-hand side of 
Eq. (jHIl) . We see that this operator is the second Hamiltonian structure 
i?2 for the Korteweg-de Vries equation 0. 

In this paper, we prove that the generating sections ip-r^ of conserva- 
tion laws [q] and the Noether symmetries ^pc for the Toda equation ^Hjl 
are related by the equation ipn = f^pc- Hence we obtain 

Theorem f theorem ITHl on page 14 1|) . (1) Suppose fi* is an integral 
for the Toda equation (j33|) such that Dy{Q^) = Vj(F). Then 
there is the operator Vj such that 

Vj = Vi O D* O K. 

In particular, consider the integral Q^ = T defined in Eq. ()58|) . 
Then 

Vi = l. 
(2) The Noether symmetries of the Toda equation are 

i 

where fi' G ker Dy are the integrals for the equation S^oda., 

is the Euler operator with respect to Q^, O is an arbitrary set of 
differential consequences to fi*, and Q is a smooth function. 

Example 9 ([SZl)- The Noether symmetries ipc of the Toda equation 
fj55|) associated with a generic nondegenerate symmetrizable matrix K 
are 

^c = DoEt{Q{x,T)) 
up to the transformation x ^^ y. 

In Sec. El we construct a continuum of recursion operators for the 
symmetry algebra of the Toda equations. These operators are param- 
eterized by arbitrary smooth functions. Although the structure of the 
symmetry algebra itself is known, see Eq. ()3()|1 . presence of the recursion 
operators gives us additional information about the Toda equation and 
permits to establish the relation between ^Toda and other mathematical 
physics equations. 
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Theorem f theorem 1^ on page I45|) . (1) Equation fjSSj) admits con- 
tinuum of local recursion operators 

R: sym£:Toda-^ sym£:Toda, 
which are 

id 

Here fij are arbitrary smooth functions and the linearizations 
ifii with respect to the integrals Q^ for the Toda equation are 

i^. = (..., J2 dnydut-D^,...). (32) 

"■ V ' 

fcth component 

(2) There is a continuum of nonlocal recursion operators for Eq. (j2EI)j 
which are constructed in the following way. Assign the nonlocal 
variables s* to the integrals Q* by the compatible differentiation 
rules 

4 = ff, 4 = 0. 

The linearizations Igi are defined by the formulas 



where i^i are calculated by Eq. ()32|) . Then the required recursion 
operators are 

R = DoJ2f,{x,s,ft)-D-'oi^., 

i 

where the functions fi are arbitrary. In general, these operators 
do not preserve the locality of elements pUjl of the symmetry 
algebra sym £^Toda ■ 

In Chapter 2, we construct the commutative Hamiltonian hierarchy 
21 of r-component analogs for the potential modified Korteweg-de Vries 
equation. This hierarchy is related with the Korteweg-de Vries hierar- 
chy for equation (jH)) that describes the dynamics of the integral (J2Z|)- 

First, we consider an example (PS|) with r = 1. Namely, we consider 
the scalar Liouville equation 

u^y = exp{2u) 

and construct the symmetry sequence for this equation. This sequence 
is identified with the hierarchy 21 of the potential modified Korteweg-de 
Vries equation 

1 3 

Also, we assign the hierarchy 03 of the Korteweg-de Vries equation 

J-t ^ 2''' XXX +3-/ J 3- 

to the initial Liouville equation. 
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Now suppose r > 1. The hierarchy 21 is constructed in Sec. ID by 
using the following procedure. Introduce the nonlocal variable s such 
that Sx = T and Sy = 1. Choose the initial function 0_i = 1. Then 
construct two sequences, 

see Eq. (jSIl), and 

We denote the first sequence (pi by 53. This is the commutative bi- 
Hamiltonian hierarchy of local higher symmetries for the potential Kor- 
teweg-de Vries equation 

St = —p Sxxx + 2'^x' y ~ / , '^i^ ■ 

i 

The second sequence ipi is denoted by 21. This is the required hierarchy 
of the higher analogs for the potential modified Korteweg-de Vries type 
equation 

Ut = (fii. 

The following theorem is valid. 

Theorem f theorem 1^ on page lUUj) . The elements (fi of the sequence 
21 C sym£^Toda compose a commutative Lie algebra: 

[2t, 21] = 0. 

Then we describe the elements yjfc G 21 such that A; < 0. We prove 
that the symmetry ip_i is defined by the Toda equation itself repre- 
sented in the Hamiltonian form 

Uy = Aio E„((a ■ exp{Ku)) dx), 

where 

A, = K~' ■ D7/ 

is the first Hamiltonian structure for the equation within the hierarchy 

21 and a = *(ai, . . . , a,,), k = \\aikij\\. 

Theorem. 1 (theorem 1771 on page l7H|l . The generators ipk of the com- 
mutative Lie algebra 21 are the Noether symmetries of the Toda equa- 
tion: ipk e sym/^Toda- 

2 (proposition EHl on page l76|) . Suppose k > is arbitrary. Then the 
kth term 0^ = E7^(/ifc da;) of the hierarchy ^ is a conserved density for 
the kth higher potential modified Korteweg-de Vries equation. 

Now we formulate the most remarkable relation between the hierar- 
chy 21 for the potential modified Korteweg-de Vries equation ^^ and 
the hierarchy ?B for scalar equation ()91|). 
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Theorem. 1 f proposition I3T] on pagel79|). Each Noether symmetry 

(^£ = D o Er(Q(a;, T)) G sym^xoda 

of the Toda equation associated with a nondegenerate sym,metrizable 
m,atrix K is Hamiltonian with respect to the Hamiltonian structure 
Ai = K,~^ ■ D~^ and the Hamiltonian Ti = [Q{x, T)]: 

ipc = AioEuin). 

2 (theorem 1^21 on page 17^ . The Hamiltonian [h^dx] for the kth 
higher Korteweg-de Vries equation is the Hamiltonian for the Noether 
symmetry (y9fc G 21 for any integer k > 0. 

In Chapter 3, we apply the methods of geometry of partial differential 
equations to the study of the properties of equations (0) and ©• 

In Sec. IHl of Chapter 3 we obtain the geometric structures for the 
?TT,-component analog 

of the nonlinear Schrodinger equation Q. By [HHI, this equation admits 
a commutative bi-Hamiltonian hierarchy of higher symmetries and an 
infinite set of the conserved densities in involution if / = id. Suppose 
/ is arbitrary, then this is not true. We compute the symmetry algebra 
for this equation if / is the homogeneous function of weight A; this case 
is actual for applications in physics. Also, we point out m'^ conserved 
currents 

r^^. = ^i^J dx + i {^l W - ^' %) dt 

for this equation. Here / is arbitrary. These currents generalize the 
known conservation laws of energy for an zth component \1/* in the 
multi-component nonliear Schrodinger equation. 

In Sec. El of the same chapter we consider the dispersionless Toda 
equation (0), 

and realize the following scheme for it. First, we find the Lie algebra 
of its classical symmetries and assign classes of exact solutions to these 
symmetries; then we construct conserved currents for this equation. 
The results of our calculations are rather cumbersome. They are found 
on pages IHTHH^ and also in the paper jHO] • 

In Chapter 4, recent geometric concepts in the theory of Backlund 
transformations and zero-curvature representations are illustrated for 
the scalar Liouville equation 

^Liou = {uxy = exp(2M)}, (33) 

which is the Toda equation associated with the Lie algebra g = s[2(C). 
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In Sec.^Jwe construct one-parametric families of Backlund (auto)- 
transformations £t for equation (jH^ . These transformations are 

e^* exp(u + u) 




^'=^-- =2e'smh(S--^- <^^' 



By definition, put 



Consider the scahng symmetry 

*^ fc>i '^ fc>i '= 

of the Liouville equation (jHHI) . Our reasonings are based on the fol- 
lowing property of X. This symmetry cannot be extended up to a 
symmetry of the covering equation Et- Recently I. S. Krasil'shchik 
described (jHI!) the scheme of generating one-parametric families of 
coverings over differential equations. We have 

Theorem (|SZI)- Let t: S —>■ S be a covering and At: S —> S be a 
smooth family of diffeomorphisms such that Aq = id and 

Tt = T o At: S ^ S 

is the covering for any t G M. Then the evolution of the Cartan con- 
nection form Urt is 

^ = lXt,uJ^ (36) 

where Xt is a Tt-shadow for any t G M and |-, ■] is the Frolicher- 
Nijenhuis bracket, see Eq. p65p on paae \lU(\ 

We prove that scaling symmetry p5|l is the Tf-shadow such that the 
evolution of the connection form Ut in the covering Tt: St ^ Suon, see 
Eq. (jHll), is given by Eq. (jHE))- The proof is straighforward. It is based 
on a useful identity in total derivatives. 

Theorem (lemma IHll on page I102|) . Let u{x) and f{u) be smooth 
functions and Dj. be the total derivative with respect to x. Denote 
Uk = D^{u{x)), k > 0, uo = u. Then the relation 



du„ 

m=l 



holds for any integer n > 1. 



In Sec. ^3 we consider the problem of constructing the pairs of solu- 
tions to the hyperbolic Liouville equation ()33|) and the wave equation 



Oxy U 
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such that these solutions are related by the Backlund transformation. 
Namely, we obtain the set of the nonlocal variables such that Backlund 
transformations are successfully integrated and each of the solutions 
related by Backlund transformation is obtained explicitly in terms of 
these variables. 

In Sec. El we analyse the correspondence between Backlund trans- 
formations and zero-curvature representations. We use two distinct 
representations of the Lie algebra q = s[2(C) generated by the ele- 
ments (e, h, /). The first representation is the representation of q in 
the traceless matrices: 

The second representation involves vector firlds on the straight line: 

P(e) = l-^, pih) = -2E.^, pif) = -E^.-^. 

We point out Backlund transformations assigned to the known zero- 
curvature representations, and vice versa. 

In Appendix El we list several geometrical methods of solving the 
boundary problems for equations of the mathematical physics. Con- 
structing solutions invariant w.r.t. symmetries of the problem is dis- 
cussed. A method based on the respesentation of the equation at hand 
in the evolutionary form is pointed out. The methods based on the 
deformation of the boundary problem along discrete or continuous pa- 
rameters are described. Among these methods, we note the direct 
iterations, the boundary conditions homotopy, the relaxation method, 
and the deformation of the initial equation. Then, comparative anal- 
ysis of the results of computer experiments in applying these methods 
is carried out. 
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Part I. The Korteweg-de Vries type equations associated with the 

Toda systems 

In Part I, we analyse the properties of the symmetry algebra sym£^Toda 
for the hyperbolic Toda equations £^Toda associated with nondegenerate 
symmetrizable matrices K. By using the canonical nonlocal recursion 
operator for the algebra sym£'Toda, we construct its commutative Lie 
subalgebra 21 generated by the local higher Noether symmetries of the 
latter equation. We identify this subalgebra with the Hamiltonian hi- 
erarchy of the higher analogs for the r-component potential modified 
Korteweg-de Vries equation associated with K. Then we relate the 
hierarchy 21 with the commutative bi-Hamiltonian hierarchy ?B of the 
scalar Korteweg-de Vries equation. The hierarchy 03 itself is a com- 
mutative Lie subalgebra of the Noether symmetries of the scalar wave 
equation. By using all these relations between the hyperbolic equa- 
tions and the evolution hierarchies of their symmetries, we prove that 
the hierarchies 21 and 03 share the same sequence of the Hamiltonians 
and demonstrate that the Toda equation itself is the first nonlocal term 
in the symmetry algebra 21. 

Chapter 1. Conservation laws and the Noether symmetries of the 

Toda equations 

In this chapter, we study the geometric properties of the two-di- 
mensional Toda equations u^y = exp{Ku) (in particular, of the Toda 
equations associated with the complex semisimple Lie algebras, see jHSl 
l83]). Namely, we analyse the relations between their conservation laws 
([HHl EHI), Noether 's symmetries of the Lagrangian functional (|HE]), 
and the recursion operators for the symmetry algebra of the latter 
equations. 

In Sec.Q we pass from the scalar Liouville equation to the hyperbolic 
Toda equations associated with generic nondegenerate symmetrizable 
matrices K. Further, we consider properties of the Lagrangian for 
these equations. Then we point out the minimal integral T G ker Dy 
and its differential span T C Dy. Finally, we assign generators of the 
Lie algebra sym£^Toda to elements fl of the kernel ker Dy. 

In Sec. m we construct a one-to-one correspondence between the 
generating sections of conservation laws and the Noether symmetries 
of the Lagrangian Z^Toda for the Toda equation. The description of the 
Lie algebra sym /^xoda C sym ^Toda of the Noether symmetries for the 
Toda equation associated with a generic matrix K is based on this 
correspondence. Also, we establish some properties which are common 
to all integrals f2* G ker Dy of the Toda equation associated with the 
Cartan matrix K oi a semisimple Lie algebra. 

Finally, in SecOlwe construct a continuum of the recursion operators 
for the symmetry algebra of the Toda equation. There are two types of 
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the recursion operators: they are either local or nonlocal with respect 
to the total dervatives. 

The exposition follows the papers jHUHnilSZI- 

1. The Toda equation 
The Liouville equation and its generalizations. The Liouville equation 

^Liou = {m?C + Mr,r, = exp(2M)} (37) 

is a model exactly integrable nonlinear differential equation that ap- 
peares in many parts of mathematics and mathematical physics. The 
study of this equation was initiated by Liouville ([69j) and H. Poin- 
care dHO])- One of the problems that they posed, the uniformization 
problem for algebraic curves (compact Riemannian surfaces), was in- 
vestigated later by Kazdan and Warner (jUj). Suppose the Lagrangian 



Cuon = -Ui^l + ul + exp{2u)) d^ A dr]] 

is regularized as described in (|101,) and then calculated on a classical 
solution of the Liouville equation. In this case, Cuou is known ( |101j ) 
to be the generating function for the accessor parameters that charac- 
terize the uniformization of a Riemannian surface of genus 0. Also, the 
Lagrangian Cuou for Eq. (jTTjl is the potential for the Weil-Peterson 
metric on the Teichmiiller space of marked Riemannian surfaces (ibid). 
Equation ()37j) plays an important role in modern field theory. In 
particular, the quantum Liouville field appeares as a quantum anomaly 
in the string theory (^21)- Further, consider the free self-dual Yang- 
Mills equation 

where F^^, is the stress tensor. Then, finding the A^-instanton solutions, 
which minimize the action for the latter equation, also leads (PZj) to 
Eq. dSZI). 

In Riemannian geometry, the Liouville equation (jTfj) is the Gauss 
equation for a conformal metric on the Lobachevsky plane (see |231 
and also [41J). 

Example 10. Consider two pointwise equivalent Riemannian metrics 

dsj = /j(x,?/)(dx^ + dy^), 

where fj > and j is either 1 or 2, on open twofolds of constant 
Gaussian curvature 

K, = -{2f,y'A\ogf, = const,. 

Here A is the Laplace operator. Namely, suppose /2 = /i ■ exp(2M), 
then u satisfies the equation 

Au = -K^h exp(2M) - K^h. (38) 
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Obviously, the two-dimensional elliptic Liouville equation corresponds 
to the case of the flat metric such that /i = 1 {i.e., K^ = 0) and the 
Lobachevsky plane metric {K2 = — 1). This interpretation of Eq. ()37|) 
is treated as a limit case {Ki < 0) in the superconductivity theory 
in description of the Abrikosov curls in the 2D model (|ll]). The 
Liouville equation appeares in the study of the Kadomtsev-Pogutse 
equations. The latter equation is a reduction of the general mag- 
neto-hydrodynamic system, where some detailes that are inessential 
in the stability problem for high-temperature plasma in TOKAMAK 
are omitted (3^). 

For any solution of the elliptic Liouville equation (j37j) there is a 
Lobachevsky plane model that is conformally equivalent to the Eu- 
clidean plane with the diagonal metric gij = 6ij. Consider the stan- 
dard Poincare model in the upper half-plane y > such that /2 = l/l/^- 
Then we get a particular solution u = — logy of Eq. (jTTjl . 



Suppose fi = 1, Ki = 0, and K2 = +1. Then formula (p8|) de- 
scribes the conformal equivalence between the Eucludean metric 6ij on 
the plane E^ and the metric g^j = exp{2U) 6ij on the two-dimensional 
sphere S"^. The function U{x, y) satisfies the equation 

U^x + Uyy + exp(2W) = (39) 

that differs from Eq. (jHTj) by the sign at the exponent (or, equivalently, 
by the multiplication of the independent variables by i = V— 1). We 
shall say that Eq. piljl is the Liouville scaA'^ -equation. 

Remark 4. Consider the orthogonal metric 

ds^ = a^dx^ + (3^dy^ 

on the sphere of raduis p. Then, Eq. ()39|) is a particular (conformal) 
case of the equation ((TT 



(a"'/3x), + {p-'ay)y + p~^aP = (40) 

that describes an orthogonal coordinate net on this sphere. Indeed, 
a = P = pexpiU) is a solution of Eq. (J4(H) if U{x,y) is an arbitrary 
solution of Eq. (jHOl)- Obviously, this solution is not unique. The pair 

a = sinV, /3 = pVy, 

where V satisfies the sine-Gordon equation 

Vxy = sin V, 

or the pair 

a = p, P = psinx 
are also solutions to Eq. (I^Uj). 
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Suppose £^ is a differential equation. Throughout this article we 
denote the number of dependent variables u^ (the unknown functions) 
by r and suppose that these r functions depend on n base coordinates 
(usually they are denoted by x^, . . . , x"). 

In this paper, we consider the generalizations of the two-dimensional 
scalar Liouville equation (jHTjl to the cases n > 2 and r > 1. We inves- 
tigate the properties of these generalizations by using modern methods 
and tools of the PDE geometry. 

We obtain a generalization of Eq. (jH7|) for ra > 2 by the following 
procedure. Recall that the Liouville equation is the condition for the 
Euclidean metric on E" to be conformally equivalent to the conformal 
metric on an n-dimensional manifold of constant scalar curvature 

seal = R = const . (41) 

We fix the value R = —2 (such that the Gaussian curvature is K = —1 
at n = 2) and preserve the correlation with the two-dimensional case 
of Eq. dSIl). We also put 

ds^ = exp(2M)^dx'=dx^ (42) 

k 

Relation ()41|) is a nonlinear PDE for the function u{x^, . . . , x"'). 

Theorem 8 ( 46 ). Condition fHT|) is 

{n - 1)Am + \[n - l)(n - 2)(gradM)2 = exp(2M), (43) 

•where A is the Laplace operator in the Euclidean space E" and the 
scalar product (■, ■) is also induced by the Euclidean metric. 



Proof. The scalar curvature R of metric ()42|) is given by the formula 

i? = ^exp(-2M)i?; 



see EHl- Then we have 



i,q 



R"" . = 3 T^ _ /5 r« . _L r« tp _ r« fp 

''qqi ^« qq ^<1 qi ' pi qq pq qii 

where 

are the Christoffel symbols. After a straightforward calculation of the 
sums in g, i, and p = 1 . . . n, we get PH|) . D 

The Toda equations. M. Toda ([Hni) considered the integrable nonlin- 
ear dynamical system described by Eq. (^, i.e., the one-dimensional 
lattice with the exponential interaction. Nowadays, vastest literature 
is devoted to the study of the properties of system ((H) and its various 
generalizations (see P^llUUllMl ITT] and references therein). In order to 
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obtain the two-dimensional Toda system out of the nonperiodic Toda 
lattice 

where the Hamiltonian H is 

i 

we replace the "acceleration" d^/dt^ with the operator d'^/dxdy. Dif- 
ferent stages of this replacement process and the discovery of the Toda 
equations associated with the Lie algebras (and the Kac-Moody alge- 
bras as well, see [321) are discussed in the papers [HZl IHSl IHHl 1111 HI El 
1231211123, see also [221 Ell- 
in what follows, we derive the generalizations of Eq. (jHTj) from a non- 
trivial viewpoint. The resulting number r of the dependent variables 
u^, . . ., u"^ is r > 1. Namely, we use the notion of the Laplace in- 
variants ( |l()()j ). To start with, we simplify the calculations by making 
the change of coordinates such that the Liouville equation acquires the 
hyperbolic form u^y = exp(M). From now on, we study the hyperbolic 
equations and their symmetries. The difference between the elliptic 
and the hyperbolic cases vanishes after the complexification. Still, the 
study of the symmetry properties of PDE does not require the presence 
of a complex structure. Therefore we assume all equations to be real 
in agreement with the definition on page El 

Following |l()()j , we calculate the principal Laplace invariants Hq and 
Hi for a hyperbolic equation 

we have / = exp(-u). Then we get 

H,"^ - D (^\ + ^ ^ + ^ = exp(n) 
^ \duyj dux duy du 



dof 



Hi= -D^(jr— ]+jr— jr— + Tr= exp M . 
OUxJ OUx OUy OU 



The definition of other Laplace's invariants follows from the equations 

D^y{\ogHi) = -H,_i + 2Hi - Hi+i, I e Z. (44) 

We recall that the quasilinear equation u^y = / is of the Liouvillean 
type if its sequence of the Laplace invariants is finite, i.e., there are 
some p > 1 and g > such that Hp = H_q = 0. One easily checks that 
the sequence H^ for the Liouville equation vanishes at once: if, = if 
i is neither nor 1. This observation motivates the above definition^. 
Now suppose —q<i<p and make the substitution Hi = exp(W*). 



We also note that the matrix Laplace invariants are defined for system (|46|l of 
hyperbolic difi'erential equations, see |100| . 
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Now, restrict equation (|I^ onto the graphs of jets of the sections lA 
in the jet bundle M'' x M^ — >■ M^. Then the total derivatives Z)j are 
mapped to the derivations d/dx^. We finally obtain the system 

Uly = 2 exp(W°) - exp(Wi), U^^ = - exp(W°) + 2 expiU'). 

Generally, we obtain the system of equations 

Ky= Yl ^^^ ^^P(^') 

j = -q+l 

such that the structure of the nondegenerate {p + q— 1) x (p + g— 1)- 

matrix K = ||/c,J| is 
II VII 

kii = 2, ki^i+i = ki,i_i = -1, kij = if |i - jl > 1. (45) 

We see that K is the Cartan matrix of the Lie algebra of type Ar-i. 
Now we shift the index i that enumerates the variables W and make 
the change ([11]) of the form U = K ■ u. Thus, we get the system 



<y = exp(-M*-^ + 2M*-u^+^), l<i<r, u° = m"+^ = 0. 

This is the two-dimensional Toda system associated with the type A^-i 
Lie algebra, which is denoted by g. Generally, we assign the Toda 
equations 

u^y = exp{Ku) (46) 

to a semisimple Lie algebra with the Cartan matrix K by following the 
geometric scheme (^El), which is discussed below. In Chapter 4, we use 
the contructions supplied by this scheme and analyse the relationship 
between the canonical zero-curvature representations and Backlund 
transformations for a certain class of differential equations. 

First, we fix some notation. Let g be a semisimple complex Lie 
algebra of rank r. By {oj, 1 < i < r} denote its system of simple roots 
and let K be the Cartan matrix of g. We thus have 

K = \\kij = 2{ai,aj) ■ \aj\'^, I < i, j < r\\. 

By K~^ = \\k^m denote the inverse matrix for K and denote its ele- 
ments by k'^K Suppose A, B E g. Then we assume that 

e = Adz + Bdz (47) 

is the fiat connection form in the principal fibre bundle G —>■ M, where 
G is the Lie group of the Lie algebra g. We thus have 

de + i[e,9] = o. (48) 

In terms of the extended total derivatives we have, then, 

[d + A,d + B]=0 ^ -dA + dB + [A, B] = 0. (49) 
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Suppose further that Hj are the Cartan generators and Ej, Fj are 
the Chevalley generators of g, respectively, for all j such that 1 < j < 
r = rankg. The commutation relations 

[Hi,Hj] = 0, [Hi,Ej] = kjiEj, 

[Hi,Fj] = —kjiFj, [Ei,Fj] = 5ijHi 

hold. We choose the following ansatz for the connection coefficients A 
and B: 

r 

r (51) 

Then from Eq. pHj) we get the relations 

db'f^ - Bai + aPJ)) = 0, (52a) 

r r 

a log ai = -J2^Ji K, d log b} = J^f^ji < (52b) 

for the coefficients of Hj, Ej, and Fj, respectively, for any j such 
that 1 < j < r. From Eq. ()52b|) we obtain 

r r 

al = Y,k'' d log h), bl = -J2k'' d log ai 
and from Eq. ()52ap we get the relation 

r 

J2k"^ 99 log {ai ■ b}) = ai ■ b) (53) 

i=i 
for any i, 1 <i <r. By definition, put 

r 

Ui = J^k"^ log {< ■ b}) . (54) 

Now substitute Eq. (jK^ in Eq. (jSHJ- We finally get the Toda equation 
(see Eq. (jlUj) ) associated with the Lie algebra q {^^)- In the sequel, 
we use the coordinates x and y as synonims for the complex coordi- 
nates z and z, respectively. Also, we suggest the following convention: 
all symmetries, conservation laws, and similar structures for the Toda 
equations are treated up to the discrete symmetry x *-^ y. 

Now suppose K = \\kij, 1 < "i, j < r\\ is a nondegenerate (r x r)- 
matrix and K~^ = \\k^^\\ is its inverse. Assume that there is an r- 
component vector a such that Oj 7^ for any 1 < i < r and the 
symmetry condition Hij = ai ■ kij = Kji holds for the matrix k = ||kjj||. 
Then we say that this matrix K is symmetrizable, see also jHl]. By 
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k we denote the operator of left multiplication by the nondegenerate 
matrix k. 

The hyperbolic Toda equations associated with a nondegenerate sym- 
metrizable (r x r)-matrix K are 

r 

^Toda = {f' = uly - exp(5^ h^u^) =0, 1 < z < r}. (55) 

i=i 

In particular, suppose g is a semisimple Lie algebra of rank r, let {««, 
1 < 2 < r} be the system of its simple roots, and denote the Cartan 
matrix of the algebra q by 



K 



2 fa,-, a 






Then we set a,- = |a,| ^; we thus have 



_ 2{ai, aj) _ 

l^ij I 12 I 12 ^ji- 



Lagrangian formalism for the Toda equation. The Toda equation £^Toda 
is the Euler-Lagrange equation in the following sence. Consider the 
action functional 



-^Toda = / -^Toda dx A dy 

with the density 

^Toda = -| X] X] 9^" l^ijU^^U^^u + a? ■ exp (^ kijU^ ^ 

Here u^ = D^(u^) and g^^ = (20) i^ ^^^ inverse of the metric tensor 
g^^ = ( 1 ) ^^^^ defines the fiat metric ds^ = dxdy on the base of 
the jet bundle n. In local coordinates, the density of the Lagrangian is 



r r r 



^Toda = "i X^ l^ijUWy - X^ «i ■ ^Xp (^ hjU^j . (56) 

i,j=l 4=1 j=l 

Assign the Euler-Lagrange equations 

E„(/:Toda) = I J2 '^^^^' \ = f^-F = (57) 

j 
to the Lagrangian ^Toda- Then these equations are equivalent to equa- 
tion ()55p since the matrices K and k are nondegenerate simultaneously 
due to the assumption Oj 7^ 0. 

Suppose £^ is a differential equation, ip G ker^^ is its symmetry and 
ip e ker£| is the generating section of its conservation law. Now we 
obtain the transformation rules for (p and ip with respect to the repa- 
rameterizations that preserve the equation's manifold S and the ideal 
S°° of its differential consequences. The following lemma is valid. 
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Lemma 9 ("Ff]). Suppose £ = {G* = 0,l<i<r}bea nonoverde- 
termined differential equation and let 

he a nondegenerate transformation of the relations that define the equa- 
tion S. Then the following two conditions hold: 

(1) The identities 

iG = A- ip, tG = iF- 'A, 

are fulfilled. Here A is the reparameterization matrix for the 
relations that describe E . 

(2) Suppose (fiG £ ker ^c* is a symmetry of the equation S and ipG ^ 
kei Iq is an arbitrary solution to Eq. (fT^ for £ = {G = 0}. 
Now, consider the reparameterization G = AF of the relations 
that define the equation S ^ {F = 0}. Then (fp = ifG is a 
symmetry of the equations again: (fp G keiip, while a solution 
il^G of Eq. (fT!?|l is transformed by the rule 

(f)G ^ipF = ^A-iIjg ^ ker tp. 

Proof. By using the definition of tlie linearization operator f^ for G = 
AF, we get 

eG = iAF = A- ip. 



Tfierefore, 



{A ■ ip)* =tpoA* = tpO 'A. 



If tlie matrix A is nondegenerate, then the condition Ig{'^g) = is 
equivalent to A ■ If{.Vg) = 0, whence we get (pp = ipG- From the 
assumption IpC'A ■ ■?/'g) = we deduce the formula ipp = ^A ■ ipG for 
the rules of solution transformation for the equation i*p{iljF) = 0. D 

Now we apply the Noether theorem (see Theorem IHl on page IT^ : 

Corollary 10. Let the assumptions of Theorem IHl and Lemma El hold. 
Let ip G ker i*p be the generating section of a conservation law for 
the Euler-Lagrange equation S = {F = 0}. Then there is the Noether 
symmetry (p G kerip of the equation S such that the relation ip = ^A-(p 
holds. 

The minimal integral for the Toda equation. Consider the Toda equa- 
tion associated with a nondegenerate symmetrizable matrix K. One 
easily checks that Eq. (J55|l admits at least one integral ( |innj ). that is, 
an expession dependent of ul whose total derivative Dy vanishes on 
this equation. Namely, we have 

r r 
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Both components T and T of the traceless energy-momentum tensor 
for the Euler equations ()57|) are well known to be of form (j3Hj) up to 
the complex conjugation: 

e = Tdx + f dy, 

see, e.g., [77]. In Sec. El of Chapter 2 we shall discuss some aspects 
of the Hamiltonian formalism for the Toda equations and derive ()58|1 
from density ()5fi|l of the Lagrangian. Meanwhile, we construct the 
differential span T of the minimal integral T G ker Dy. By definition, 
put 

T, ^ Di{T). 

The differential consequences Tj to the functional T generate the sub- 
space T C ker Dy within the kernel of the total derivative Dy. Indeed, 
any smooth function Q supplies the functional 



Qix,T)=Q{x,T,Ti,...,T^) EkeiD 



y 



We say that the nondegenerate symmetrizable matrix K is generic if 
the integral T in Eq. (J58|l is a unique solution to the equation Dy(T) = 
on the corresponding Toda equation (|KK|l . 

We emphasize that by a special choice of the matrix K one can 
achieve the situation such that the functional T will not be a unique 
integral for the Toda equation (j33|) . The criteria for the equality 
dimkeri)y = 2 to hold at r = 2 are found in the paper [HH]- From 
now on, by ri C kerD^^ we denote the subspace in kerDj^, which is 
differentially generated by all solutions fi* to the equation Dy{Q^) = 0. 
The total number of these solutions is denoted byg: 1 < i < q < r. 
We always assume Q^ = T. 

Example 11 (|^). Consider the Cartan matrix K = (_i "2) '^^ ^^^ 
Lie algebra s[3(C). Then the Toda equation (jHH|) associated with K 
admits two integrals. The first integral Q^ = T, which is defined in 
Eq. (jHEl), acquires the form 

^^ = i^lf - "X + i'^lf - «L - «L, (59a) 

while the second integral is 

n' = uL. + ul ■ {ul - 2ul) + (ulf ■ ul - u\ ■ {ulf. (59b) 

The integral fi^ is a solution to the equation Dy{Q) = but does not 
follow from ^2^. 

By [20], a necessary and sufficient condition for r nontrivial inde- 
pendent solutions fi* of the equation Dy{Q'') = to exist is that K be 
the Cartan matrix of a semisimple Lie algebra q. The Toda equations 
associated with q are exactly integrable (pB]). 
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The symmetry algebra of the Toda equations. In this subsection, we 
assign classes of infinitesimal symmetries of the Toda equation to the 
functional span J7 of differential consequences to the integrals f2*, which 
were introduced in the previous subsection. Then, in Sec. |2 we find 
out which out these symmetries are the Noether symmetries of the 
Lagrangian ^Toda- 

By A = I A* I denote the vector of the conformal dimensions ( |H2] ) 

r 

A^ = J] k'^ 
i=i 

of the Toda fields exp(-u) = *(exp('u^), . . ., exp(-u'")). This notation is 
well-defined due to 

Proposition 11 ([8J. (1) The transformation 

X I— > X{x), 

y ^ yiv), (60) 

u\x, y)^ff= u\X, y) + A' log X'{x)y{y) 

is a finite conformal symmetry of the Toda equation S^oda- 

(2) The Lagrangian C^oda. = / -^Toda dxAdy is invariant with respect 
to this change. 

r 

(3) By definition, put 13 = ^Oj-A*; under diffeomorphism ()60|) . the 

i=l 

component T of the energy -momentum tensor Q is transformed 
by the rule 

Remark 5. The symmetry properties, Eq. (jHIHl and (jHT|) . of the Toda 
equation associated with the Lie algebras were considered in the pa- 
per [Hj. The present formulation of Proposition ^2 is an extension of 
the cited result to the Toda equation (jKKjl . which is now associated 
with an arbitrary nondegenerate symmetrizable r x r-matrix K. The 
coefficients A and P are already adapted to this general case. 

The infinitesimal variant of Proposition ^2 is 

Proposition 12. 1 (|75[l77j). The infinitesimal components of confor- 
mal symmetries fl(j()|l of the Toda equation are 

up to the transformation X ^^ y. Here f is an arbitrary smooth function 
and 

D = M, + A ■ D, (62) 
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is a vector -valued, first-order differential operator. 

2 (|4^). Each point symmetry (fl is a Noether symmetry of the La- 
grangian /^Toda^ 

B^/ (Lxoda da;) eimdh- 

3 ([ZZmi)- The functional T , which was defined in Eq. ()58|) . is a density 
of the Hamiltonian for the infinitesimal conformal symmetry ip^ : 

ipl = A,- E„(r • fix) dx), A, = k-' ■ D-\ (63) 

The third statement in Proposition El was formulated in local co- 
ordinates in the paper [ZZj. In Sec. El of Chapter 2 we establish the 
relation between the canonical Hamiltonian formalism for the Toda 
equation and the Hamiltonian operator Ai above. In Theorem 1521 on 
page E21 we find out that equality (jU^ is the root part of the infinite 
sequence of relations between the hierarchy of the higher analogs of the 
Korteweg-de Vries equation, see Eq. Q on page[TTJ and the commu- 
tative hierarchy within the Noether symmetries algebra for the Toda 
equation (this hierarchy is constructed in Secjl)). 

Lemma 13. Suppose K is a symmetrizable r x r-matrix and the op- 
erator n is defined by Eq. ()62|) . Then the relations 

(64) 



O ft O n = Dr O KOHO D,, 



hold. 



Proof. Here we express the first relation in local coordinates: 



oD 



^ij ^xy 



kijexpiy^^kgu^] 



\ui + A^-D, 



"'X ^xy 



U 



XX ^y 



U 



xy - 



Ux T Uxxy ~r ^ ^xxy 



y^ hjul exp f ^ kiiu^j \~\^ ^*^'^^^ ^^P W^ ^'' 



lU 



3,P 



D,. 



= D.J. o \ux + /\DJ^ o Dy, 

Q.E.D. 

The second relation is deduced from the first one by using Lemma JHl 
and the Helmholtz condition 



KE(£a 



>) 



■'E(£Toda) 



n 



that holds since the matrix k = *k is symmetric. 
Corollary 14. The vector-functions 

</. = n(0(x,f2)) (65) 

are symmetries of the Toda equation: ip G sym £^Toda for any function 
that depends on an arbitrary subset Vt of the integrals ^2* = Z)^(f2*) G 
ker i)j^. 
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Formula ()65|) gives the description of the whole symmetry algebra 
sym£:Toda: 

Proposition 15 ([75). (1) Let the assumptions of Lemma^^hold; 
then any symmetry of Eq. ()55|) is (jU3j) . 
(2) Suppose further K is such that there exist q independent inte- 
grals n* G ker Dy, where 1 < q < r. Put 

r 

f' = exp(^ kijU^) and f] = kij ■ f\ 
i=i 

Assume further that there are r constant (r x q) -matrices M^ = 
||(M(j)*H, l<i<r,l<j<q, and a constant (r x q)-matrix 
A = ||A'||, 1 < i < r, 1 < j < g, such that rank A = q and the 
equations 

A^ /^" = (M^)r /^ /I (M,)f = {Mp)t /f 

hold. Now, construct the first -order, (r x q)-matrix differential 
operator 

r 

U = J]M« ■ < + a ■ D^ (66) 

a=l 

and consider a vector = |0*(x, f2)|, 1 < i < q. Under these 
assumptions, the sections defined in Eq. (j^ exhaust all sym- 
metries of the Toda equation i^Toda- 

Therefore, in both cases we have 

sym£:Toda - iv' = O'jCJy'ix, ft) mod {x ^y)}, 

where the number of columns in the operator D is equal to the number 
q of the independent integrals ^2' and is an arbitrary vector. 

Corollary 16. Any solution tp to the equation ip{ip) = for system ()55|) 
is of the form 

^ = «:(n(0(x,n))). (67) 

We stress that the problem of finding the integrals Q is primary for 
the Toda equations, and the search for the symmetries ip as well as 
the selection of the Noether symmetries ipc follow this problem. We 
also note that each conformal symmetry ()60|) of the Toda equation is 
Noether, i.e., preserves the Lagrangian Z^Toda- Still, not each section ip 
of type (J65|l is a Noether symmetry of Eq. (J55|l . 

2. The Noether symmetries of the Toda equation 

First, we recall Example El on page fT^ and obtain an important prop- 
erty of the Toda equation. This property permits the application of the 
generating sections machinery in description of the conservation laws 
and the Noether symmetries of the Toda equation. Namely, we have 
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Lemma 17. The Toda equation f^Toda is i-normal. 

Proof. By Example El on page El it is sufficient to represent the equa- 
tion Sioda in the evolutionary form. Let 

^ = x + y, 7] = x-y 

be the new independent variables. The way to choose them is such 
that the equations 



r 

M^^ - M^^ = exp [y^ hjU^ 



i=i 



hold. Now, put f* = M^; then the equations Sev C J^{R^, M^"') of the 
form 



U^ = V , v^= u 



^5 - exp 



r 



i=i 
are the required evolutionary representation of Eq. ()55|). D 

From relation (|H7|l and Corollary El we deduce the relation ip = 
K(pc between the Noether symmetries and the generating sections of 
conservation laws for the Toda equations. This observation allows to 
specify a property common for all the integrals f2* for Eq. (j33|) . Namely, 
suppose 

d;,(fiMa;) = -Vi{F)dxAdy 

for any admissible i and consider the conservation law [t]] = [Q{x, 
ft) dx] . Then we have 

dnQix, ft) dx = -J2^. ^- ° ^*(^) d^ ^ ^2/- 
By definition, the generating section ip^i of the conservation law [rj] is 
^. = -E(-l)^(^0*°^i(^) =-E(^0*°E,.(g). (68) 

Now we compare Eq. (jHTj) with Eq. ()(J8|1 and, by using Theorem IHl and 
Lemma ini we obtain 

Theorem 18. (1) Suppose fi* is an integral for the Toda equa- 
tion (j33|) such that Dy{Q^) = Vi(-F). Then there is the operator 
Vj such that 

Vi = Vi o D* o k. 

In particular, consider the integral fl^ = T defined in Eq. (jSHJ- 
Then 

Vi = l. 
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(2) The Noether symmetries of the Toda equation are 
(^£ = no^V*oE^,(Q(x,f2)), 

i 

where fi* G ker Dy are the integrals for the equation EioAa.-, 



Eo, = > i-lVDl 



j>0 J 

zs i/ie Euler operator with respect to Q^, O is an arbitrary set of 
differential consequences to f2*, and Q is a smooth function. 

Example 12. Again, consider the Toda equation (jKH|l associated with 
the root system A2. Then we have r = 2, at = |ai|~^ = 1 for i = 1, 2, 
and 

2 -A ._i 1 /2 A r /I 



Put n = M^. + A • Dx, see ()fi2|l . The integrals Q^ and fi^ are defined 
in Eq. ()59|). One easily checks that i5y(i7^) = D* o k{F). Therefore, 
Vi = 1 and 

Dy{Q'^) = -D^oD* ok{F), 

thence, V2 = — -Dx- We emphasize that the integral ^2^ is not equivalent 

to -D^{Q^). 

Remark 6. In Theorem IT^ we established the relation Vi = 1 to hold 
for the minimal integral T, which is defined in Eq. fl58j) . Now we restrict 
all previous reasonings onto the subspace {Q{x, T)} C kerDy of the 
kernel of the total derivative. Thus, we consider the subspace that is 
generated by the integral T and its differential consequences. Then, by 
Theorem IHl the conservation laws [Q dx] for the Toda equation ^Toda 
are in one-to-one correspondence with the Noether symmetries 

ifc = noETiQix,T)). (69) 

In other words, the Noether symmetries ipc of the Toda equation ()55|) 
associated with a generic nondegenerate symmetrizable matrix K are 
()69|) up to the transformation x ■i^ y. This statement is an extension of 
the relation (jHE]) between the Noether symmetries and conservation 
laws for the scalar Liouville equation onto r > 1. We note that our 
scheme of reasonings is essentially more simple than the straightforward 
proof for the Liouville equation (ibid). 

In the rest part of this section, we generalize the method of the paper 
[SSI- W^ solve the equation 

directly with respect to ip and thus find out which of the symmetries 
if = n(0(x, T)) G sym£^Toda preserve the Lagrangian /^xoda for the 
Toda equations. We follow the paper 
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Detailed analysis of some algebraic aspects for Eq. ()55p shows (tech- 
nically, we calculate the term E2 (^Toda) of the Vinogradov's C-spectral 
sequence, see pi!) that the following two conditions are equivalent: 

©^(-^Toda) = <^=^ Eu(B^(£Toda)) = 0, 

i.e., a density is a total direvgence iff its variation is zero. We note 
that this question was not discussed in [Sn]. Further on, we get 

Lemma 19. The relation E^i(3<^(£Toda)) = is equivalent to the con- 
dition 

-E„(Dy(T)-0(x,T,...,T^)) = O, 
where ip = □(0) and the total derivative Dy is not restricted onto f^Toda- 

The proof of Lemma El is based on multiple use of the relation 

r 

i=i 

Lemma 20. The highest order m of the derivative T^ in the set of(j)'s 
arguments is even: m = 2fi, and (p = (f){x, T, . . . , T2^). 

Proof. Put Gi = 7— (B^(>CToda)) and by u-^r^i) denote the derivative 
D^ o Dy{u^). Then we have 



duU.,.) ' ' ' ' dT, 



5 

m 



whence (—1)'^ + (— 1)™+^ = and thus m is even. D 

Lemma 21. The derivatives 

8^6 



dT^m 




vanish for all I such that fi < I < m. 




Proof. It suffices to check that 




dGi dG, 


dGi 


dU2m+4: dU2m+2 


dul+e 



are equal to 0. This can be done straighforwardly. D 

Still, the derivative 



dT^dT2^ 
is, in general, nontrivial owing to 

Lemma 22. The identity 

E„(-D,(r) ■ Et{P{x, T, . . . , T^))) = 
holds for any function P. 



44 ARTHEMY V. KISELEV 

Proof. Indeed, we have 
DyiT)-ET{P)dxAdy = {Dy{T)J^P*{l)) = {iP{Dy{T)),l)+d,^ = 
(9g^j(^)P(x,T, . . . ,T^), 1) + d^7 = {DyiP), 1) + dh^ e kerE,, 

where the couphng ( , ) takes values in the horizontal 2-forms uj = f ■ 
dxAdy, d/i = ^ da;*®-Dj is the horizontal differential, d^'-f G ker E„ is an 

i 

exact form, and both the evolutionary derivation and the linearizations 
are evaluated with respect to T. D 

Proposition 23. A symmetry cp = n{(f){x, T , . . ., T^)) G keiip is a 
Noether symmetry of Toda 's equation (j33j) iff the following two condi- 
tions hold: 

(1) m = 2/i, /U > 0, and 

(2) the function = Et{Q{x, T, . . . , T^)) G imEr lies in the image 
of the Euler operator with respect to T, where Q is an arbitrary 
smooth function. 

Proof. We show that G iniEr by induction. Choose P = P{m;x,T, 



. . ., Tn) such that 



and put (j) = (j) — Er(P(m)), then we get 



0. 



d^P{^) ,_^y . d^ ^7Q) 



dT 

By Lemma |2ni the equation 

E,(D^(T)-0(x,T,...,r„_2)) = O 

holds. By using Lemma El we choose P{m) in accordance with ()7()|1 
and apply Lemma |221 Therefore, we decrease the order m = 2/i to 
with the step 2 inductively. Finally, we obtain 

= ET^(5^P(2z;a;,T,...,T,)). 

i=0 

The proof is complete. D 

Thus, we have obtained another description of the Noether symme- 
tries' class fl(i9|l for the Toda equation assigned to a generic matrix K. 

3. Recursion operators for the Toda equation 

In this section, we construct a continuum of the recustion operators, 
which are either local or nonlocal with respect to D^, for the symmetry 
algebra of the Toda equation. Although the structure of the symmetry 
algebra itself is known, see Eq. (jHSI), presence of the recursion operators 
gives us additional information about the Toda equation and permits 
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to establish the relation between ^Toda and other mathematical physics 
equations. 

The explicit method by J. Krasil'shchik that allows construction of 
the recursion operators for symmetry algebras of differential equations 
was briefly described in the Introduction on pages ITsHTSl The theorem 
below contains the result of application of this method to the Toda 
equation (jKK|) associated with a nondegenerate symmetrizable (r x r)- 
matrix. 

Theorem 24 (^3). (1) Equation fl55j) admits a continuum of lo- 
cal recursion operators R: symi^xoda -^ symi^xoda; which are 



R = noJ2f-^{x,n)-Dio 



«.j 



Here fij are arbitrary smooth functions and the linearizations 
iai with respect to the integrals fi'^ for the Toda equation are 

...,Y.^-D„..). (71) 



fcth component 

(2) There is a continuum of nonlocal recursion operators for Eq. (j33j) . 
which are constructed in the following way. Assign the nonlocal 
variables s* to the integrals fl^ by the compatible differentiation 
rules s]. = fi* and s* = 0. The linearizations Igi are defined by 
the formulas 

where ifii are calculated by Eq. ()7H1 . Then the required recursion 
operators are 

i 

where the functions fi are arbitrary. In general, these operators 
do not preserve the locality of elements ()65|) of the symmetry 
algebra sym£:Toda- 

Proof. First, enlarge the set of the dependent variables m^ by introduc- 
ing the nonlocals s* such that their derivatives are 

4 = n\ sy = 0. (72) 

(In fact, any definition of s* is allowed if it is compatible with the 
condition s'^^y = Sy^ = 0). Then, extend the total derivatives: 

D^ = D^ + J2^^-^, Dy = Dy, 
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such that [Dxi Dy\ = 0. The Cartan's flat connection is now defined on 
the equation 

^oo ^ {D^/\s') = DliSt), bl o by{s') =0, fc > 0; 

D,(F) = 0, |a|>0}. 
The Cartan's generating l-forms 

for the recursion operators -Rxoda satisfy the determining equation 

^Toda(^Toda) = 0, 

where -^Ipoda is the restriction of the hnearization of Eq. (j3^ onto 
H(.' (£). From the factorization 

^Toda O D = Dj, O n O Dy 

in Eq. ()64|) it folfows that any Cartan's 1-form 

roda = n($^/.(x,s,...,D,^-(.^))-dc(^:.(.0)) (73) 



i>0 



is an element of the kernel ker ^ipoda' whence we get the statement of the 
theorem. In particular, suppose that each /j does not depend on the 
nonlocal variables s. Then the resulting recursion operator is local. D 

Example 13. Consider the integral T defined in Eq. (jHE))- Its lineariza- 
tion is 

r 

iT= (...,^«:.X-^x-«i-^x.---) =a*oD,ok. (74) 

^ V ' 



ith component 

Now, introduce the nonlocal variable s: we set s^ = T and Sy = 1. The 
compatibility condition for the variable s is 

Sxy = 0. (75) 

Finally, we construct the recursion operator 

Apply the operator -Rxoda to the translation Ux G sym£^Toda- We obtain 
the symmetry sequence 

ipk = □(0fc-i) 

that corresponds to the sequence of functions 

0_1 = 1, 00 = Si, 01 = -pSs + |S^, 

02 = P^S5 - l(3sl - 5pSiS3 + Isl, 
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etc. In Chapter 2 we shall investigate the properties of this recursion 
operator -Rxoda together with the properties of the symmetry sequence 

in detailes. Meanwhile, we claim that the symmetries ipk of the Toda 
equation, which are obtained by multiple action of -Rxoda to the trans- 
lation ipo = u^, are local, Hamiltonian, and commute with each other. 
Also, in the next section we shall establish the relation between the 
symmetry sequence 21, the Korteweg-de Vries equations © and (|T6|l . 
and recursion operators p8j) for the latter equations. 
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Chapter 2. The Korteweg-de Vries hierarchies and the Toda 

equations 

In this chapter, we construct the commutative Hamihonian hierarchy 
21 of muhi-component analogs for the potential modified Korteweg-de 
Vries equation. This hierarchy is identified with the commutative Lie 
subalgebra in sym£^Toda composed by certain Noether's symmetries. 
Also, we discuss some aspects of the Hamiltonian formalism for the 
Toda equations themselves and establish a link between the hierarchy 21 
and the higher Korteweg-de Vries equations for Eq. Q. The exposition 
follows the papers |^ l^. 

We start our study of the relation between the Toda equation (|^ 
and classical equations and (fT^ of mathematical physics, the Kor- 
teweg-de Vries equations, with the following 

Example 14. Consider the hyperbolic Liouville equation 

^Liou = {uxy - exp{2u) = 0}. (76) 

The minimal integral, see Eq. ()58p. for the latter equation is ( P^I1UU| ) 

T = ul-U2, Dy{T) = 0. (77) 

Introduce the nonlocal variable s such that 

s^ = T, Sy = 1, (78) 

and put 

i3 = 2ui. (79) 

Then, consider the symmetry 

¥^ = (wi + px)(r) 

of the Liouville equation and calculate the evolution of the variables u, 
•&, T, and s along this symmetry. We get 

ut = -|m3 + ul (potential mKdV) (80) 

Tt = -|T3 + 3TTi (KdV) (81) 

St = -|s3 + fsi (potential KdV). (82) 

The Miura transformation (see [ZEIISI) acquires the form 

^1 = t2T t i^', ^t = ±T2 - (^Ti + ^iT) (83a) 

T = t|^i - i^'. (83b) 

One can treat relations (jH^ as Backlund transformation between the 
Korteweg-de Vries equation, see Eq. (jHTj) . and the equation 

^, = -1^3 + 1^2^^ (modified KdV). (84) 

The signs '±' and '=]=' in Eq. (|HH|l are induced by the symmetry ^ ^-^ 
— T? of equation (jHH). This discrete transformation provides Backlund 
autotransformation ([IHlinSl) for Eq. (|HT|). 



^ip_l,a_i — ^ip-\ + Cf— 1 ■ Q_ 
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The recursion operator 

^Liou = Dl - 2ui + D-\^ D, 

is common for both equations (jTUI) and (IHUj). see [IHllini- This operator 
generates the commutative Lie subalgebra 

SlLiou = spanjj(v?fc = RuonM, ^0 = Ui) 

of local higher symmetries of the potential modified Korteweg-de Vries 
equation, see Eq. (jHDI). 

Now, to the Liouville equation we assign the variable D such that 

Da; = exp(2M), (85a) 

S, = {Vy = 0^} (85b) 

and therefore the compatibility condition v^y = ^yx holds. Thence, the 
equation £^Liou is represented in the evolutionary form: 

Mt_i = V-i = f > (86) 

where the variable t_i = y is the parameter and (f-i is the shadow of 
the nonlocal symmetry 

d_ 

such that a_i = 0^. The solution to Eq. ()85b|) is = —{y + X{x))^ . 
Apply transformation ()6U|) of the form y = y{y) to the "time" y. Then 
we obtain the potential 

y'jy) 

{X{x)+y{y)) 

for the genreral solution u = \ log \3x of the Liouville equation (see jHSl 
ES]): 

^ i{X{x) + y{x)yi 

We recall that functional (f77j) is continuous on formal divergences 
u -^ ±oo of solution ()87|) . see the paper [IH] and references therein. 
We estimate the scheme of constructing the general solution for the 
Liouville equation by using the potential to be really laconic and pro- 
ductive. A similar approach was used in |l£j| to derive the Cole-Hopf 
transformation for the Burgers equation, see Eq. ()2()H|1 on page 11221 
Solution (|H7j) of Eq. (|7H|l is the mapping 

r: {Xy = 0, yx = 0}^Suou. 

Consider the evolutionary vector field Qut defined in Eq. (jHDI). This 
field can be lifted onto the inverse image of the covering r. Hence we 
obtain the equation X = and the Krichever-Novikov equation 

Xt = -2X3 + SX^/Xi = -2X1 ■ {X, x}, (88) 
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where {A',x} is the Schwartz derivative. Also, the evolution equation 

Of = ai = O3 - fOsOr^ 

holds. 

The evolution equations that appear in Example El are systemized 
as shown in the diagram below: 



The Krichever- 
Novikov Eq. (l2Ul) 




Eq. 


IS} 






Potential 
KdV Eq. (IHS) 






Eq. 


GHJ 




KdV Eq. dHU) 



Eq. 



EU 



Eq. 

JH3 



The Krichever- 

Novikov Eq. (jHH|), 


Eq. 


iSTl 


Petential 
mKdV Eq. flHHll 


Eq. 


i7^ 


Modified 
KdV Eq. (|H^ ■ 



4. Analogs of the potential modified Korteweg-de Vries 

equation 

4.1. Constructing the hierarchy 21. First, we formulate an impor- 
tant property of the integrals f2* G kerDj^| for a Liouvillean type 
equation S. Namely, we specify the evolution of these integrals along 
symmetries ip G sym£^ of the equation S. 

Lemma 25 ( |innj ). Consider a symmetry field 3^^ of a Liouvillean 
type equation S. Then the evolution 9,^(^2*) of an arbitrary integral 
Q^ G ker Dy for £ is an element of ker Dy again. 

Proof. Indeed, we have 

Dy{3^{9})) = 3^{Dy{n')) = 9^(0) = 0. 

D 

Example 15. Consider a symmetry ip = n(0(x, T)) of the Toda equa- 
tion, see Eq. (J65|) on page EHl Then the relation 



f^ = 3dw (T) = {-P Dl + TD, + D,oT)\ 



(89) 



holds. Here 



P = ^ai-A\ n = u^ + AD^, and A* = ^ k'^ . 



i=l 



Also, consider the nonlocal variable s that was defined in Example El 
by the equations s^ = T and s„ = 1. Then the evolution St of the 
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nonlocality is described by the formula 

Suppose that the sequence ipo, ipi, 02 of symmetries of the Toda 
equation is assigned to the specially chosen set of the functions in 
formula (|U3j). We set 0_i = 1 and obtain the symmetry ipQ = n(0_i). 
Then, we calculate the evolution T^_^ that corresponds to this sym- 
metry. Now we equal the next function 0o to the evolution of the 
potential s. Then, we obtain the function (pi and the symmetries ipi, 
(p2 in a similar way. The result is shown in the next diagram: 

f3Ts + 3TTi 

^KdV 



7V_i - Ti 



(90) 

The following evolution equations are met in diagram (J90|l : the Korte- 
weg-de Vries equation is 

^Kdv = {Tt = -m + 3r ■ Ti} (91) 

(we recall that its recursion operator is p8a|l on pageUHl); then, we 
also get the potential Korteweg-de Vries equation 




and the equation 



^pKdV = {St = -PS3 + |si} 



{ut = D{T{ui,U2)}. 



(92) 



(93) 



^pmKdV 

Suppose K = II 2 II is the Cartan matrix of the Lie algebra q = 3^2 (C). 
Then the Toda equation (j^ associated with this algebra is the hyper- 
bolic Liouville equation (fTUj) . and equation (j^ is nothing else than the 
scalar potential modified Korteweg-de Vries equation (jH^ . Further on, 
suppose that the matrix K in Eq. (|55|l is an arbitrary nondegenerate 
symmetrizable (r x r)-matrix, not necessarily the Cartan matrix of a 
semisimple Lie algebra q of rank r. Then we get the r-component sys- 
tem of the third-order evolution equations with the cubic nonlinearity. 
In local coordinates, this system is 



u 






a^ 



{k 



■pq 



'p,(j=l 

where 1 < i < r; we recall that 



«:« + 2{A^k 



pq 



5. 



i,q 



I'^xx'^x 



2 AX. J, 



A* 



E* 



QjIlCl CLpnJpq 



Ctq rvqp . 
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The analogs of the potential modified Korteweg-de Vries equation for 
r = 2. In this subsection, we consider evolution systems (jUHj) in case 
r = 2 and the matrix K is symmetric. 

Example 16. Let K be the symmetric matrix 

thence, 

A ' '' 



A + 2 VI, 
Introduce new dependent variables 

u = u^+u^, V = {\ + 2){u^ - u^); 

the inverse transformation is defined by the formulas 

1 1 1 2 1 1 

u = -u+ —— — f , u = -u 



(94) 



2 2(A + 2) ' 2 2(A + 2) 

Then, the normal form of Eq. ()93|) is 

2 2-A A+2 n 2-A , 

"* = -aT^"^ + (AT^"^"^ + ^"^ + 4(aT^"^"- 

A-pZo 2 — A o 

Further on, make the scaling transformation 

t = (A + 2)^-t, M = (A + 2)"^-M, v = v, (95) 

preserving the same notation x, t, u, and v. 

Proposition 26. In coordinates ()95|) . equations (jMj) become linear in 
A anc? acquire the form 

Ut = -2m3 + 2f if 2 + |ui + ^Mlf 1 + A ( -f it;2 + |mi - ^MiVi ) , 

. X (96) 

vt = -2u2Vi + lujvi + |fi + \(-U2Vi + \u\vi - \vl\. 



Remark 7. Consider the vector coefficient of A in Eq. (j96|) . Then this 
vector- valued function is not a symmetry of the fiow that stands at A° 
(and vice versa). Namely, the following statements hold. 

Proposition 27. (1) Consider the flow 

Ut = -2u-i + 2fif2 + l^f + |-uifi, 

Vt = -2u2Vi + \u\vi + \vl 

at A° in Eq. (|nHj) . T/ien its symmetries 

2 J (^, X, U, V, Ui, Vi,U2, V2, Us, V3) 
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of order < 3 are generated by 
'tut + lxu^\ /uA fuA fl\ fO 

tv, + ixvj' U;' W' W' V 

These generators are the scaling symmetry, two translations, 
and two shifts, respectively. The commutator with the scaling 
symmetry maps the translations along t and x to themselves. 
(2) Consider the flow 

ut = -vxV2 + \u\ - \uvo\, 

vt = -U2V1 + lulvi - Ivl 

at X^ in Eq. (jM|) . Then the symmetries 

2 ) (i, X, U, V, Ml, Vi, U2, V2, Us, V3) 



of order < 3 of this flow are generated by the sections 
tut + \xuA [ut\ fuA fl\ fO 



tv, + ixvj' Vv' Vv' voy' Vv' 

which are the scaling, the translations, and the shifts, respec- 
tively. 



The conservation laws for system (j96|) in the normal form are few. 

Proposition 28. Suppose A is generic {see Remark ^ below) . Then 
there is a unique generating section 

2] {t,X,U,V,Ui,Vi,U2,V2) = i A-2\.^ 



of order <2 of a conservation law for the flow in Eq. (j96jl : this section 
corresponds to the conserved density 

9 A — 2 9 

In Sec. Q we shall demonstrate that equation ()94|1 is Hamiltonian 
with respect to the operator Ai = K~^ ■ D~^ (recall that the ma- 
trix K is assumed symmetric and therefore Oj = 1). Now we get the 
Hamiltonian representation of system (jM|) by using the transformation 
rule for Hamiltonian operators with respect to transformations of the 
dependent variables. 

Lemma 29. Consider a Hamiltonian equation 

Ut = A{E^{n[u])). 

Let u = Qu be a nondegenerate transformation of the dependent vari- 
ables. Then the equation 

u = A{Y.r,{n[u])) 
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holds for the Hamiltonian operator 

A = Q-A-'Q. 
Example 17. Consider Eq. (j^Hjl . Then we get the diagonal operator 

/ t \ z— A / 

where the density hi of the Hamiltonian T^ dx for the Korteweg-de 
Vries equation acquires the following form in the coordinates m, v. 

hi = ^^{{X-2)vl-{X + 2)ul + A{\ + 2)u^:)\ 

Proposition 30. Again, suppose A is generic {see Remark IB)). The 
symmetries 

^2] {t,X,U,V,Ui,Vi,U2,V2) 



of order < 2 of equation ()96|) are generated by the translation and the 
shifts, which are 



::■ JV ri. (-) 



respectively. 



Remark 8. We emphasize that the symmetry algebra (and even the 
classical symmetry algebra) of system ()9(i|l depends on the initial ma- 
trix K essentially. The cases A = — 1 (the matrix K corresponds to the 
Lie algebra A2) and A = ±2 (the matrix K is degenerate) are special. 
In what follows, we analyse them separately. 

First, let A = —1, that is, K is the Cartan matrix for the Lie alge- 
bra A2. Then substitution (J95|l is 

u = u^ + u^, V = u^ — u^. 

After routine transformations, we obtain 



ut = -u-i + |fit;2 + |m? + ftiifi, 

Vt= - |u2fl + lujvi + |fi. 

Proposition 31. (1) The symmetries 



(98) 



V 



{t,X,U,V,Ui,Vi,U2,V2) 



of order < 2 of flow (I98j) are generated by the sections 



lu + pOgvJ' \vj' VOyl' \1 
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The first symmetry is nonpolynomial. The other two are the 
translation and the shift, respectively. The nonpolynomial sym- 
metry flow can he also represented in the form 

ul^ = exp{u + 3uu). 

(2) The generating sections 

,2 1 {t,X,U,V,Ui,Vi,U2,V2) 

of order < 2 of conservation laws for Eq. ()98|) are 

U2\ f ^{u,v,Vx) 

3v2j ' \2d'$/dv + 2t;2f r^ d'^/dvi - ui d'^/dvi 

where the function \E' satisfies the equation 
d^ 1 d^ 1 

T^ + -Vi -r ^ = 0. 

du 2 9t;i 2 

Now, for the first time within this paper, suppose that the matrix 
K is degenerate. Quite untrivially, this degeneracy can occur in two 
distinct ways. 

First, assume that A = 2 and K = (ll)- Then system ()94|1 is in 
triangle form^: 

(99) 

Vt = -U2V1 +U-^Vi. 

It consists of the potential modified Korteweg-de Vries equation and 
the auxiliary dispersionless component. We note that the variable v 
admits an arbitrary shift. Indeed, the section *(0, /(f)) is a symmetry 
of Eq. (jnn|) if / is arbitrary. In addition, there are two more symmetries. 
They are the translation and the shift of u, see Eq. ()97|1 . 

Proposition 32. If X = 2, then the generating sections 

,2 j {t,X,U,V,Ui,Vi,U2,V2) 

of order < 2 of conservation laws for the flow in Eq. ()99|) are 

U2\ /exp(M)\ /exp(-M; 

0^' [ )^ [ 

The corresponding conserved densities are 

pi = -\ul, p2 = exp(M), P3 = - exp(-M) 



The classical symmetry algebra for the Toda equations associated with this 
matrix is composed by the sections 

= aux + I3{x) 1, 

where the constant a and the function j3{x) are arbitrary. 
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respectively. The last two densities belong to the nonpolynomial con- 
served currents 

r]2 = exp(-u) dx + (2^2 — u^) exp(-u) dt, 

r]3 = — exp(— m) dx + (2^2 + ul) exp(— -u) dt 

for the potential modified Korteweg-de Vries equation 

ut = -2^3 + ul 

whose right-hand side is polynomial. 

Now, suppose A = —2, whence K = ( _2 2 ) • Then equation (j^ is 

Ut = -2m3 + 4:ViV2 + uivf, 

(100) 

Vt = v^. 

Proposition 33. The symmetries of system (jUUj) for A = —2 are de- 
scribed by Eq. (j^Tj) . The generating sections 






2 I (t, X,U,f, Ml, f 1,^2,^2) 



of order < 2 of conservation laws for the flow in right-hand side of 
Eq. (fmn|l are, then, 



Here the function \1/ is subject to the equation 

-- - 2vl ^— - 3vf ^— + vivl 6v^V2 * = 0. 

at ov ox 0V2 

Factorizations of the recursion operators. Now we return to the prob- 
lem of constructing the commutative hierarchy 21, which is associated 
with the initial Toda equation ^Toda- 

By construction, the symmetries ifi are related by the recursion op- 
erator 

i?Toda = no/^;^o£r (101) 

for the Toda equation (J55|) . We met this operator in Example fT3l on 
page Eni We fix^ this notation -Rxoda onwards. 

In this subsection, we demonstrate that each symmetry 

0k = □(0fe-i) e sym£:Toda 

defines the next function (pk explicitly by the relation 0^ = 9i?fc(s), and 
therefore, 0k+i = n(0fc)- We also prove that the sections 0^ compose 
the hierarchy of the potential Korteweg-de Vries equation ()92j) . 



As we have already noted in Chapter 1, we treat the Toda equations, as well as 
the structures on them, up to the symmetry x ^^ y. 
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First, we note that the functions 0_i, 0o; and 0i are mapped suc- 
cessfully one to another by the recursion operator {^^) 

RpKdY = -(3Dl + 2si - D-' o S2 (102) 

for equation (jHSI)- Recall that this operator was obtained in Introduc- 
tion while illustrating the Cartan generating forms method d^) by 
I. S. Krasil'shchik. 

Lemma 34. The following two decompositions of the recursion opera- 
tors hold: 

i?Toda = □ O 4, RpKdY = 4 O n, (103) 

where is = D^^ ° ^t by Eq. (j^ and the linearization ix is defined in 
Eq. dZH). 

Proof. The first decomposition holds by construction. To establish 
another one, we use the relations 



tr o (n, + A/},) 



= i ^ Kijwi D^ o u{ + I ^ tti I {^ kiju\^ D^ o u\ - 2Dl o wi } + 

i,j i j 

i,j i j 

= - ^ aiA' Dl + ^Y^ ai{kij{u\ui + mX) - 24} + 
+ 2YI (^i{kijiu\u\ + u\u{) - Aui} D^ + 

= -(3Dl + Ti + 2TD^, 
whence follows the second decomposition. D 

We draw attention to the fact that the representation of the scalar 
operator -RpKdv in the form of the product of the vector-valued operator 
n and the row is of length r seems to have been unnoticed in the 
literature. 

Now we identify infinitesimal symmetries of differential equations 
with autonomous evolution equations and, by using Lemma EH we 
extend diagram ()90|) infinitely upwards. The result is displayed in the 
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(104) 



The right-hand sides (pk of the evolution equations St^. = (pk are 
generated by the recursion operator -RpKdV; therefore they are higher 
symmetries of the potential Korteweg-de Vries equation £^pKdv- 

The times t, in the equations Ut^ = ipk and St^ = 4>k within dia- 
gram ()1U4|) are correlated. Namely, we have 

Theorem 35. The evolution B(pt.(s) of the nonlocal variable s along the 
higher symmetry ipk = D{(f)k-i) of the Toda equation (jKH|) coincides with 
the evolution (j)k that is defined by the kth higher analog st^ = -RpKdv(l) 
of the potential Korteweg-de Vries equation (jH^ : 



9r 



Rr 



^a{4>^,_l)[S) = (Pk = ^pKdVWk-1, 

Proof. From decompositions pOHj) it follows that 
^pKdv(0fc) = D-^ (9n(<^,)(^)) • 



D 



Corollary 36. 1. Each shadow 0k in the covering r^, see Eq. ^7^ . 
over the Toda equation can be reconstructed up to the true nonlocal 
symmetry Bi^^.^^ of the Toda equation. 
2. All sub diagrams 

0k+2 < 4>k+l 



^Toda 



R 



pKdV 



D 



in diagram ()lU4p are commutative: we have -RToda('^fc+i) = n(0fc+i), 
and the relations 



-Rxoda — Do -RpKdV O n , 



R 



pKdV 



□"' O i?Toda O D 



hold. 



Remark 9. In the paper jlOj, the recursion operators -Rxoda and RpKdv 
were noted to be conjugate to each other in the scalar case 
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Proof of the locality for the symmetry sequence 21. All previous rea- 
sonings do not imply the locality of the elements (pk of the sequence 
03 composed by higher symmetries of the potential Korteweg-de Vries 
equation (jHH). Indeed, these symmetries lie in the image of the opera- 
tor D~^. Therefore, the elements ipk+i = □ {(t>k) of the sequence 2t are 
not yet proved to be local, too. 

Proposition 37 ([20J). Recursion operator p()2|l for the potential Kor- 
teweg-de Vries equation (jH^ generates the local in T sequence of the 
higher symmetries 

where 0_i = 1 is the shift of the dependent variable s by a constant. 

There are several ways to prove Proposition EZI The method by 
I. S. Krasil'shchik (,£1J) is based on the weak nonlocality of the recur- 
sion operator -RpKdv (see also Eq. ()142|) on page | 



Corollary 38. The symmetries ipk = n(0fc-i) G 21 of the Toda equa- 
tion (f^3j) are local and depend on the derivatives m^, \a\ > 1 for all 
k>0. 



We also note that in the paper |1HJ we analysed the case r = 1 and 
established the locality of the higher symmetries 21 C sym£^Liou di- 
rectly, leaving apart the discussion upon the properties of the potential 
equation (jHSI)- 

4.2. Commutativity of the hierarchy 2t. A classical example of 
evolution equations that admit a commutative symmetry subalgebra is 
given in the following lemma. 



Lemma 39 ([2E])- Suppose £ is the scalar evolution equation 

Ut = Uk + f{Uk-2,---,u) (105) 

such that f is a polynomial. Consider the Lie subalgebra 

{f G sym£^ I ip = f{ua)) C sym£^ 

of its symmetries that depend on the variable u and its derivatives only. 
Then this subalgebra is commutative. 

We denote by 03 the minimal Lie subalgebra generated by the sym- 
metries (pk of the potential Korteweg-de Vries equation; here k > —1. 
From Lemma inni it follows that the algebra ?B is commutative: 

{0fc,0«} = 9?ifc(0O - 9?i,(0fc) = 0, 

therefore 03 coincides with the linear span of its generators ipk'- 

03 = spanig.(0fc | k > -1). 

Now we analyse the commutation properties of the symmetries ip = 
n(0(x, T)) of the Toda equations (|^. 
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Lemma 40. Suppose y?' = D (0'(x, T)) and </?" = D (0"(x, T)). Then 
the Jacobi bracket {if' , <f"}, which was defined in Theorem ^ of the 
symmetries ip' and ip" is 

{^',^"} = n(0{i,2}), 

where the bracket 

0{i,2} = B^'(0") - B^"(0') + D,{<P') <P" - <P' D-M") (106) 

on the arguments of the operator D is induced by the Jacobi bracket 
{ip'^if"}. Moreover, 

0{i,2} =0{i,2}(a;, T) 
by Eq. 



Further on, we denote by 21 the minimal Lie algebra generated by 
ifk for all A; > 0. 

Theorem 41. The Lie algebra 2t C sym£^Toda ^s commutative: [21, 21] = 
0; and thence 

21 = spani[j((^fc I A; > 0). 

Proof. Commute two symmetries 3^^ and 9,^^. , apply the resulting 
evolutionary vector field to the variable s and take into account the 
relation between (pk and ipk- Hence we get 



[9^fei' 9^feJ(s) = 3v.i(</'te)-9 



vk.yvk^, 



9it(0fe) - 9S>t(</>^J = {^k.^ M = 0. 
Since T = s^, thence 

[9^.,,9^.J(^) = 0. (107) 

Now estimate the evolution of the integral T by using Lemma 1301 
First, consider the bracket {(fki, ^k2} and then calculate T^^^, ^, y Re- 
call that ifki = n(0fci-i) and ipk2 = ^(^te-i)- Therefore, 

{^kl,^k2} = □(0{fcl,fc2})5 

where 

'Pik^M] = 

^<t>k^ (0fc2-l) - 9*^2 (0fcl-l) + £'x-(0fei-l) 0fc2-l - 0fcl-l 5x(0fc2-l) 

owing to Eq. pOfij) . From Example El it follows that 

9{^fci,^fc2}(^) = 9n(<^{fci,fc2})(^) = 

= {-P Dl + T D^ + D^oT)i(P{kuk2})- (108) 
Comparing Eq. pU7p and Eq. ()108p . we get 

{-(3Dl + TD^ + D^o T) <P{k,M} = 0- (109) 
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In the left-hand side of Eq. fll()9j) we obtain the operator B2 in total 
derivatives whose coefficients are elements of T. We apply this operator 
to 0{fcj_fc2}(T, . . ., T^[kiM)) ^^^ g^t in the right-hand side Thence, 
4>{ki,k2} ~ ^^^ therefore 

The indexes ki and k2 are arbitrary, therefore 2t is commutative. D 

Proposition 42. Let £^(o) = {ut^ = (po{U(^)} be an evolution equation. 
Suppose that the symmetry ipk G sym£(o) is assigned to each k > and 
assume that ipk is independent of the time to explicitly: ifk = ipk{ua). 
Then the following two statements are equivalent: 

(1) The algebra 21 = span]g((y9fc \ k > 0) is a commutative Lie alge- 
bra: {ipk,^Pi} = 0. 

(2) The evolutionary vector field 9<^j is a symmetry of the au- 
tonomous evolution equation £(k) = {ut^ = ipk} for each k, 
l> 0. 

Proof. First, identify the evolutionary vector field B<^^ with the au- 
tonomous evolution equation Ut,, = (pk- Then, consider the equahty 

In its left-hand side we have 

since cpi does not depend on any time tk- In the right-hand side of the 
equality we get 

Therefore, the comuutativity condition {(pki ^1} = for the symmetries 
(.pk and ipi is equivalent to the determining equation 

Thence ipi is a symmetry of the equation £i^k) for any k, I > 0. D 

Corollary 43. For any A;, / > 0, the following propositions hold: 

(1) The sections (^^ G 21 are not only symmetries of the Toda equa- 
tion, but also symmetries of all equations S(^i) = {ut^ = (pi}: 

0k e sym. 8'^y 

(2) The recursion operator -Rxoda is common for the whole tower of 
the evolution equations £{iy. 

-RToda e Rec £{i) . 
In particular, 

-Rxoda = -RpmKdV- 
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Remark 10. In the paper j30], the scalar potential modified Korteweg- 
de Vries equation was considered in the gauge Ut = u^, + u^. In this 
case, this equation shares the recursion operator with the sine-Gordon 
equation u^y = sinu. 

Remark 11. By using Corollary EHl we deduce that the section (p'_-^ = 
const is a central extension of the commutative Lie subalgebra 21 C 
sjmS^ of symmetries of the evolution equation S(^i^ for any / > 0, 
although const ^ sym£^Toda- 

Remark 12. By using formula p06|) . we establish a curious property of 
the Korteweg-de Vries equation (J9T|l . Recall that the proof of Theo- 
rem HH is based on the equality (p{ki,k2} — ^^^^ holds for all elements 
0fcj and 0^2 of ^he hierarchy ^. Consider the first Hamiltonian struc- 
ture Bi = Dj. for Eq. (j^Tj) . We know that the Hamiltonians of the 
higher Korteweg-de Vries equations are in involution with respect to 
this structure. Thence we have 



(Et(H,), Bi o Er(7^,)> = ^^ 0^ ■ D,(0,) G imD,. 



Now we see that the last two summands in Eq. ()106|) are of this form. 
Therefore the sum 



9<^fc,+i(0fe) - 9<^fc,+i(0fci) e imD^ 



is always a total derivative. The feature we observe is that each of these 
summands above is a total derivative by itself, therefore generalizing 
the involutivity property. Namely, by straighforward calculation we 
obtain Vtij such that the relations 



^(D.(0,))=D.(%) 



hold for the elements of the hierarchy 03. Several Vtij that correspond 
to i, J < 2, and the initial elements 



.1 = 1, 00 = Si, 01 = -/3S3 + fSi, 

02 = /5^S5 - \(3sl - 5/?siS3 + hi 
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of the sequence 03 are found in the table below: 



i = -1 


J = l 
3T 


J = 2 
-5pT2 + f T^ 


i = 


3m2 
2^ 


5T3 - 5/3TT2 + I/3T2 


i = 1 


37^2 1 3rp2 
~2^1 + 2^2 


|/32T| + f T^ - ff3T'T2 


i = 2 


- 3/5T1T3 + f T| 

+ f T^ - ffSTT^ 


- |/3'T| + f T^ + I5/32TT1T3 
- f /5T3T2 - I/32T2T2 



Obviously, Qi _i = Qi = for ahi. 



5. The Hamiltonian formalism for the Euler equations 

In this section, we consider the problem of constructing commuta- 
tive Hamiltonian hierarchies of evolution equations associated with hy- 
perbolic Euler systems (in particular, with the wave equation s^y = 
or the Toda equation (J55|l . see jHH IMl IHO]- Also, we interpret 
the canonical coordinate-momenta formalism, which is widely used in 
mathematical physics (see [01 EZl QEj and references therein), within 
nonlocal Hamiltonian operators language of the jet-bundle framework 
( [TUl 1121 ES]) preserving the distinction between the coordinates and 
momenta. Next, we consider hyperbolic Euler equations [i.e., we spec- 
ify the ansatz for the Lagrangian density) and obtain the differential 
constraint between the dependent variables m* and the momenta vcij. 
Then, we identify symmetries of these equations with autonomous po- 
tential evolution equations; also, we investigate the relation between 
Hamiltonian operators for potential and nonpotential evolution equa- 
tions that describe the evolution of coordinates and momenta, respec- 
tively. The aim of our reasonings is to assign new jet bundle tt' and 
the jet space J°°(7r') to the initial Euler equation £e-l C J^ij^), see 
Eq. (fT3j). such that the latter equation becomes an evolutionary vec- 
tor field that belongs to x(7r'), while the evolutionary fields that cor- 
respond to the momenta m become elements of the adjoint module 
k{Ti'). Finally, we relate commutative Lie subalgebras of the Noether 
symmetries for hyperbolic Euler equations and pairs of Hamiltonian hi- 
erarchies composed by potential and nonpotential evolution equations. 

5.1. Canonical formalism. Consider an abstract 2r-dimensional dy- 
namical system with r dependent variables m', momenta trij, the spatial 
coordinates x, and the time t, defined by the Poisson brackets (PITTj) 



:iii^ 
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{mi,m,}^ = 0, (110) 

{u\x, t),mj{x', t)}^ = A"-- 5{x - x')] 
where A is an (r x r)-matrix differential operator in total derivatives 

with respect to a; (a = (o"i, . . ., (T„) is a multiindex) and the skew- 
symmetric bracket acts as the derivation with respect to any of its 
arguments. Assume that Ti = [H dx\ is a Hamiltonian with the density 

H{x) = H{u{x),Ua{x);m{x),D^m{x)). 

Then the dynamics 

u= {u{x),H{u{x'),m{x'))}^, 

m = {m{x),H[u{x'),m{x'))}^ 

of the variables u and m is obtained in a standard way: 

u = {u{x),H{x')}a= i y]{M(a;),D^m(a;')}A-7T77; — T^TT ^^'• 

Here C{x) is a small contour around the point x. We emphasize that 
this language of local coordinates and the ^-distribution is spread in 
many works on the Hamiltonian formalism application in the field the- 
ory; in the sequel, we pass to the invariant exposition soon. Meanwhile, 
integrating by parts we obtain 

u = Ao-^. (112a) 

dm{x) 

Similarly we obtain the second relation 

m = -Ao——-. (112b) 

ou{xj 

These reasonings motivate the following definition of the variational 
bracket of two Hamiltonians with the densities H and H', respectively: 

we set 

5u 5m 5xn 5u 

Remark 13. Usually, all dependent variables are treated uniformly. In 
our case, the momenta can be absorbed by the additional variables 
u^^^ = m.j for 1 < j < r. In other words, we enlarge the total number 
772 = 2r of the coordinates u^ twice. The operators A and —A are also 
united in the {m x m)-matrix 

^= (a* 
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such that the dynamical equations take the form 

^ 6{H{u)) 

u = Ao — — , 

ou 

where the variation 6/6u is calculated with respect to the new vector 

u = *(m, m). Due to this reason, Remark [T^ motivates definitions EHHl 

see page ITIIl in the Introduction. 

A traditional approach (jTOl 1121 EB]) to description of Hamitonian 
PDE dynamics does not appeal to any distinction between the depen- 
dent variables, coordinates u and momenta m. Nevertheless, by using 
the double set of m = 2r dependent varuables u and m we describe 
several remarkable properties of models of mathematical physics, for 
example, of the Korteweg-de Vries equation (jHH), the modified Korte- 
weg-de Vries equation (jHH), and similar equations. These aspects are 
discussed in Sec. [3 below. 

6. Hyperbolic Euler equations 
Consider a first-order Lagrangian 

£ = /LK.„«.;..y)dxAd, 

with the density 

L = -i ^ RijU^ul + H{u; x, y), 

where k is a nondegenerate, symmetric, constant (r x r)-matrix. We 
note that the notation R is correlated with the general exposition. Sup- 
pose r = 1 and H = 0, then we obtain the wave equation (|75|l (see 
Sec. EI), li h = K = \\aikij\\ and the function H is defined by for- 
mula pi7|l . then we get the Toda equations. 

Choose the independent variable y for the "time" coordinate (the 
nondegeneracy condition implies dL/duy ^ 0), leave x for the spatial 
coordinate on the base M of the new jet bundle vr' : M x M— >]R with the 
old fiber coordinate m, and denote by iTij = dLjdu^ the jth conjugate 
coordinate (momentum) for the jth dependent variable u^ for any 1 < 
3 <r: 

r 

mi = -|^Kijui. (113) 

i=i 
The differential constraint, Eq. ()113p . between coordinates and mo- 
menta for the initial equation ()15|) is our main tool that allows to 
construct the Hamiltonian structures. 
Consider the Legendre transform 

HdxAdy= ( m, - — ) -C 
\ duy/ 
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and assign the Hamiltonian 

n{u,m) = [Hdx] 

to the Lagrangian £. We decompose its density H into the sum of two 
equal summands and use the relation 

u = -2R~^D-\m) 

in one of the components, in agreement with pi2|) : 

The hyperbolic Euler equation 

£e-l = {£„(£) = 0} 
is equivalent to the system 

(5m ' ^ 5u 



^y = TT' '^y = -^77 (H^) 



with respect to the canonical Hamiltonian structure A = 1. Owing to 
the relations 

5 15 

5u 2 "> — V — ■' dm 

Ax Ai 

the dynamical equations are separated: 




u. 



y 



A,oE^{[H[u]dx]), 



my = -lAioEr^{[H[m]dx]). 

Example 18. Suppose that r = 1, R = ||1||, and the Hamiltonian is 
trivial: H = 0. Then the Hamiltonian operators 

Bi = D-\ Bi = D., 

assigned to wave equation (f7S|) are mutually inverse. 

Now we apply these reasonings to the Toda equations ()55|1 . which 
are assigned to Lagrangian (jHSI)- It turns out that their Hamiltonian 
representation (jllSj) hints the existence of the minimal integral T, see 
Eq. (j3Hj) . The conservation of T is induced by the conservation of 
Hamiltonian density ()117|) for Eq. (j33j) . 

We start with a general 

Lemma 44. Assume that the density H of a Hamiltonian Ti = [H dx] 
for the Hamiltonian evolution equation Ut = A(E„(7i)) does not de- 
pend on the time t explicitly. Then the density H is conserved for this 
equation: 

[Dt{H)dx] =0. 
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Proof. By using the condition dH/dt = 0, we calculate the derivative 
integrating by parts, we obtain 



j:S(-^rD.{^),AoE.in)) 



since the Hamiltonian operator A is skew-symmetric, A* = —A, we 
get 

= -{AoE^{n),E^{n)). 

Again, by using the definition of the Euler operator £„ and integrating 
by parts, we have 

= -(9^.E.(H,(ff),l> = -ACK). 

Moving the result to the left-hand side of the initial equality, we finally 
obtain the required condition 2Dt{H) G imD^;. D 

Example 19 ([59j). Choose the coordinate y for the "time", define the 
momenta (see ()11H|1 ) m = dL/dUy-. 



mi - 2 



i^KijU^ (116) 

and obtain the density -ffxoda of the Hamiltonian Tixoda: 



/ 2 \ / 

-f^Toda(M, ^) = -i X] '^* ^^P ( ~ ^x\'"^i) ] - I ^ Oi exp f ^ kijl 

We have already decomposed the density H into the sum of two com- 
ponents that depend on m, and u^ , respectively. Then, the canonical 
Hamiltonian representation of the Toda equation £^Toda is 



. i l^-f^Toda 



Therefore, in terms of the dependent variables u, we have 

r r 

-f^Toda(w) = ^ai expf ^ kijun (117) 

i=l j=l 
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and 

u = AioEu{nTod.{u)), (118) 

where Ai = k^^ ■ D^^. Recall that an evolution representation for the 
Liouville equation (f7B|) was obtained earlier in Eq. (IHUj). 

In order to correlate the resulting expressions with ()58p. we apply 
the transformation x ^^ y to the Haniiltonian equation ()118|1 . consider 
the Hamiltonian density ()117|l . which is independent of y, and use 
Lemma EH Then we obtain 

r r 

whence follows expression ()58|1 . 

By using the Hamiltonian representation ()118|1 . we describe the neg- 
ative elements (fk oi the sequence 21 with k < 0: 

Proposition 45. The right-hand side of the Hamiltonian evolution 
equation 

Uy = k-^ O D-^ O E„(7i:Toda) (HIHI) 

is the inverse image of the translation ut^ = ipo with respect to the 
mapping -Rxoda; the translation Uy is the element (p_i of the sequence 21. 

Proof. By Lemma IH^ we have 

RtoUV-i) = □ (9«.(5)) = ^{Dy{s)) = 0(1) = u,. 

a 

Therefore, representations ()118|) and ()118l|) of the Toda equation 
provide the translations 

Uy = (y9„i = k'^ O Z)-l O E„(iJToda) 

and 

U^ = ^0 = k'^ O Dy^ O E„(ifToda) 

that are related by the recursion operator -Rxoda: 



^Toda 



^Toda 



ip^2 = 0{T) ^, — Uy = v?-i 7==^ Vo = Ux — °-^ Vi = □(T) 



Toda /Xm J 

loda 



The operator -R^oda obtained from -Rxoda by using the discrete symme- 
try X ■^^ y generates the symmetries ^9^ G 21 of the Toda equation with 
k < -1. 

Now we return to the discussion on the Hamiltonian formalism for 
the Euler equations Se-l = {E„(£) = 0}. We identify the symmetries 
(f of these equations with the evolution equations Ut = ip. Consider 
a symmetry ip{ux, Uxx, • • •) of Eq. (fTHjl . i.e., the potential evolution 
equation 

ut = (p(ux,Uxx---), (120a) 
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then the induced evolution rrit of the momenta is described by the non- 
potential equation 

mi = -ifi;-L)^.(V(m,m^, ...))• (120b) 

In addition, assume that the evolution ip is Hamiltonian: 

. _15H . _ 16H 

"""2 5^' """^27^- 
Then we have 

u = R~^ ■ D^^ -— , 

At 

m = -^D^- K — , 

ii 
i.e., both equations ()120|) are Hamitonian simultaneously and their 
Hamiltonian structures Ai and Ai are mutually inverse. 

We note two classical examples of the pairs of evolution equations 
that admit mutually inverse Hamiltonian operators. 

Example 20 (IM]). The potential KdV equation (j^J^ is Hamiltonian 
with respect to the operator Bi = D~^, while Bi = D^ is the first 
Hamiltonian operator (see ()l()|l ) for the KdV equation (jHJ)- O^is can 
easily verify that Eq. (j^ is compatible with the wave equation (fTSj) . 
i.e., the right-hand side 0i in the potential KdV equation St = 0i is a 
symmetry of s^y = 0, while the nonpotential equation ^T\ describes 
the evolution of the momentum T = s^ (up to the constant factor 
-1/2). 

The potential and nonpotential modified Korteweg-de Vries equa- 
tions, Eq. (|^ and ()127|) . respectively, supply another example. In 
fact, the succeeding subsections 17. II and 17.21 contain the detailed anal- 
ysis of these two examples. We demonstrate that the first pair defines 
the symmetry hierarchy 03 for wave equation ()75|1 and the second pair 
is related with the hierarchy 21 of symmetries of the Toda equation (jKK|) . 

Now we return to Eq. ()121|) and note an important property of the 
Hamiltonian equations 

u = A^o E^{n) 
and 

It turns out that we have already met the expression in the right-hand 
side of the latter equation in Lemma El on page El Lemma JSl assigns 
the generating sections ipr^ = E„(?7o) to the densities rjo of conservation 
laws [t]] for evolution equations. Formulas ()12()j) give an interpretation 
of Lemma the generating sections ip describe (up to the sign) the 
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evolution m of the momenta m along the Hamiltonian symmetries (p of 
the initial equation £. 

Moreover, we relate two pairs of mappings of different types. Assume 
that there are recursion operators Ru and R^n for different evolution 
equations ()120a|) and ()120bj) . respectively, and consider the mappings 

Ru : X — >■ X, Tu'. k ^ k 

which produce symmetries and generating sections of conservation laws 
for the same equation Eu-, respectively. The relation ( |42j ) 

% = K (122) 

holds for evolution equations; see, for example, diagram ()141|) on page 15^ 
In terms of the ()120j) , Eq. p22|) means that the skew-symmetric recur- 
sion operator 

A. = (_fl: t) (123) 

is defined for the Hamiltonian representation (jllSj) of the initial Euler 
equation S = {£„(£) = 0}. This situation is realized for the above- 
mentioned pairs of equations, see ()130|) on page ESI We also note that 
the condition Ry = R*^ is valid in this case. 

Now we study the problem of constructing a bi-Hamiltonian hierar- 
chy (to be more precise, a pair of hierarchies with respect to u and m) 
by using the recursion operator ()123|) . Namely, we find the conditions 
for the operator Ar to induce the second structure {m, m}^ on the 
Hamiltonians Tii in the diagram 



7^0 ^1 ^ 



2 



s 

5m 



J_ 
5m 



S 
(5 m 



ifQ — ^ Vl — ^ V2 



R R R 

such that the Jacobi identity ()7b|l holds. The commutativity condition 
[(y9j, ipj\ = is sufficient for the Jacobi identity to hold for bracket (0), 
supplied by the operator A^] the sections m(. = iIji also commute in this 
case. We denote by il the minimal Lie algebra generated by the sections 
ti(. = ipi. We emphasize that, in general, we require these operators 
()123j) to provide the Lie algebra structure (j7b|l on the Hamiltonians 
Tii e H'^{'k) but not on the whole horizontal cohomology group H^{tt). 
Therefore, our concept extends Definition [3 of a Hamiltonian operator 
since we treat the particular case of a hyperbolic Euler equation £ and a 
sequence of its Hamiltonian symmetries (pi whose Hamiltonians are Tij, 
respectively. Further on, suppose il is a Lie subalgebra of the Noether 
symmetries of the Lagrangian C. This assumption is sufficient for the 
existence of the Hamiltonians Tij such that 

(Pi = R —^ = 1 ■ ^, ^p^ = -R —^ — = -1 ■ ^. (124) 

dm dm ou du 
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Indeed, the existence of the conserved densities is prescribed by the 
Noether theorem (see Theorem IHl on page IT^ . 

We obtain a more usual description of the pair of the Magri schemes 
(see [70]) by using the sphtting ()121|) : 





(125) 



Here the operators Ai and Ai are defined by constraint ()113|1 between 
u and m, while the second operators A2 and A2 originate from the 
relations 

R = A20 A^^ i?* = ia o ^r^ 
respectively. In the sequel, we preserve the notation and mark the 
operators Ai^2 and -81,2 for the nonpotential equations ()127j) and (jHH), 
respectively, with the 'hat' sign. 

The notion of the i*- covering in the category of differential equations 
was introduced in [42^ , and the unifying approach towards the recursion 
operators R: x ^ x, the conjugate recursion operators T: k ^ >c, 
the Hamiltonian structures A: k ^ x:, and the symplectic structures 
A: H —^ k was elaborated for the Magri scheme ()125p technique. The 
return from diagram ()125|1 to ()124|) provides new interpretation of these 
structures and their interrelations. 

Remark 14. The correlation between the generating sections of conser- 
vation laws for the initial Euler equation Ee-l = {F = R ■ E„(£) = 0} 
(these sections are denoted by ipc) and for the same equation in its 
evolution representation (jllSj) (we denote these sections by ifj) is given 
by the diagram 

-D-' , s-i (126) 



V' 



V^£ 



<^, 



where the first arrow follows from the definition of a generating section: 

dhV = V(F) da; = V o D^{D-\F)) da;, 
and, therefore, 

^^ = V*(1), ^ = -Z},oV*(l), 
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while the second arrow is induced by Lemma IHl 

7. Properties of the Korteweg-de Vries hierarchies 

In this section, we illustrate the concept ([11112011001) formulated in 
the preceding section. In what follows, we consider the bi-Hamiltonian 
equations (jOH), (jOH), (|^- and the r-component analogs 

^mKdV = Ri=/^xO/>on(T)} (127) 

of the modified Korteweg-de Vries equation (jH^ . We also introduce 
the new variable 

d = K-u^ (128) 

that is subject to the latter equation. This dependent variable differs 
from the momenta 

by the constant factor (—2). This gauge simplifies the calculations and 
resulting expressions. 

Now we study the relations between the equations i^pKdv and £^Kdv, 
and between ^pmKdv and i^mKdv, respectively. We claim that the po- 
tentials u and s for d and T are such that the linearizations 

t^ = L>x O K, trp = Ux 

are equal to the first Hamiltonian structures Ai and Bi for the equa- 
tions £^mKdv and £^Kdv, respectively. The first Hamiltonian structures 
Ai and Bi for £^p(m)Kdv are inverse with respect to Ai and Bi. The 
Hamiltonians for the equations £^(p)Kdv are well-known. The Hamilto- 
nians Hk for the equations £^(p)mKdv are described in final Theorem IH^ 
on page EHl 

7.1. The Korteweg— de Vries equation. Let Cs be the Lagrangian 

Cs = -[^s^Sydx A dy]. 

Assign the wave equation 

see Eq. (j73|) on page HOI to Cs and consider the commutative Lie sub- 
algebra 

<B = spanij(i?^+^v(0_i), 0_i = 1, A; > 0) C sym^, 

of its symmetries. We identify this subalgebra with the hierarchy of the 
higher potential Korteweg-de Vries evolution equations, see Eq. (lU^ . 
We continue using the variable T = s^, whose evolution is described 
by the elements of the nonpotential Korteweg-de Vries hierarchy for 
Eq. (jni}, instead of the canonical momentum rris = —^Sx for the wave 
equation. By p21|) . these two hierarchies admit the pair of mutually 
inverse Hamiltonian operators Bi = B^^ = Dx for the equations i'xdv 
and ^pKdv, respectively. The coinciding elements within the initial part 



SYMMETRIES AND CONSERVATION LAWS IN INTEGRABLE MODELS 73 

of the Magri schemes (see Eq. ()125p ) for Exdv and £^pKdv are underhned 
in the diagram below: 




ipKdV 



B2 



h[T] = [i(/5T2 + T3) da;] -^ -(3s, + ^sj r, 




ho[T]= [iT'dx] 



h-i[T] = [Tdx] 




B2=-/3 Di+D^oT+T-D^ 



(129) 



The initial terms of this diagram are given by the Hamiltonians whose 
densities are known from the vastest literature: 

/i_i = [Tdx], ho = l^T^dx], h = iliPT^ + T=^) dx], etc. 



Now we consider an example to the relation ()122|) . Namely, the 
identity 

-RpKdV = -Rxdv 

is satisfied by the recursion operator -RpKdv for the potential Korteweg- 
de Vries equation (fH^ and the recursion operator 

-PDl + 2T + Ti • D;' 



R 



KdV 



for the Korteweg-de Vries equation (jnTjl . The recursion -RpKdv is fac- 
torized by the Hamiltonian structures as follows: 



B2 



D-'o{-(3Dl + D,oT + T-D,). (130) 



B2 
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The second Hamiltonian structure B2 for Eq. ()9H1 is a factor in the 
right-hand side of Eq. pHOj) . Recall that this structure appeared also 
in Remark [T31 on pagelSUl Another important property of this structure 
is the correlation between the equation £^KdV; the Toda equation £^Toda, 
and the Virasoro algebra (|2n|)- We analyse the latter property in 
Remark [T71 see the end of this subsection. 

The Hamiltonian operators Bi and -Bi are inverse, therefore the Ma- 
gri schemes ()125|) for Eq. (|U^ and ^T\ are correlated. 




hk[T] 



E, 



hk-i[s] ^—>fk-i 

id 




'■pKdV 



Tt, = $fc = ^fc (131) 



hk-i[T] -^st^^^ 



^pKdV— ^KdV 



^KdV 



^ifc_l = "^fc-1 = ^k-\ 



B2 



^Kd 



Here (^^ and $fc are symmetries and ipk and ^^ are the generating func- 
tions of conservation laws for the potential Korteweg-de Vries equation 

and the Korteweg-de Vries equation 

Tt, = $1, 



respectively. We conclude that the opposite sides of diagram ()131|1 are 
identified with the one-step vertical shift. 

Corollary 46. We conclude from diagram (jlHlj) that the symmetries $fc 
of the bi-Hamiltonian hierarchy for the Korteweg-de Vries equation are 
the gradients iIju of the Hamiltonians hk [s] for the potential Korteweg- 
de Vries equation (|^. and vice versa. 
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The equations £^Kdv and £^pKdv share the same set of the Hamiltonians 
with the densities hk\T] and hk[s\, respectively. Thence we obtain the 
following property of symmetry subalgebra 2t of the Toda equation: 

Theorem 47. The generators (fk of the commutative Lie algebra 21 are 
the Noether symmetries of the Toda equation: 

ipk e sym^Toda- 

Proof. Indeed, we have 

21 9 V9fc = n(0fc_i) = n o ExikklT]) e sym^Toda C sym£:Toda 

by Remark IHl on page |121 □ 

Remark 15. The conservation laws [77^] associated with the generating 
sections 

V'J:f = /€-noET(M 

(see Eq. ()69|) on page 1^ and diagram ()126|) ) are exactly the Hamilto- 
nians hk dx of the higher Korteweg-de Vries equations St^ = (pk- 

This fact, [hk dx] G //^^(^Toda), was noted in a similar situation in j22l 
§10] but required a nontrivial proof there. 

By Lemma ^^ the densities hk of the Hamiltonians for both Korte- 
weg-de Vries equations ^T^ and (jH^ are conserved on the correspond- 
ing higher analogs Tt,. = Dx{(pk) and St^ = (pk of these equations. We 
have 

A,(M = 5.(l^f^). 

We claim that the hierarchy 03 is composed by conserved densities 
for the hierarchy 21 of potential modified equation (|^. We need the 
following useful lemma to prove this claim. 

Lemma 48 ([in])- The relation 

E{B^iC)) = 3^{E{C))+i;{E{C)) 

holds for any (f & x and C G A"(7r). 

Proof. Suppose A G CDiff (x, A"(7r)). By multiple integrating by parts, 
we transform the expression ^A((p) to 

where the order of Aq is zero and A'((p) G CDifF(x, A'*(7r)); then by 
using Eq. (fT^ we obtain 

E(A(^))=£;(A*(1))+£1.(,)(^) 

for any section ip E h. Now, let C G A'*(7r) be a form and put A = 
£c: X -^ A"(7r). The linearization •^e(£) = ^e(/;) ^^ the image of the 
Euler operator is self-adjoint, hence we obtain the equality 

E(4M)=^;(E(/:))+£e(z:)M, 
whence follows Lemma HHl D 
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Remark 16. In particular, from this lemma it follows that any Noether 
symmetry ipc of a Lagrangian £ (we recall that the condition B,^(>C) = 
holds on J°°(7r) in this case) is also a symmetry of the Euler equa- 
tion (fT3j) that is assigned to C, i.e., 

sym£ C sym£^. 

Lemma HHl gives an explanation why the inverse statement is wrong. 

Proposition 49 ([51 ). Suppose k > is arbitrary. Then the kth term 
0fc = Er(^fe dx) of the hierarchy ^ is a conserved density for the fcth 
higher potential modified Korteweg -de Vries equation. 

Proof. First, we recall Theorem 1351 that describes the correlation of the 
times within the hierarchies 21 and *B. Second, we apply Lemma HBl 
that provides the rule to commutate the Euler operator and an evolu- 
tionary derivation. Third, we recall that the higher symmetries 0^ are 
independent of the jet variable s. Thence we obtain 

Dt.M = Q^A^Ah)) = EA^^) -i;^iBAhk)) = 

, ^ / 

=0 

= -J2i-^yDiof^.EAh,))eimD^. 
j>o V i / 

Therefore the density (f>k is conserved on the equation £(^k) = {utk = 
(fk}- □ 

Remark 17. The second Hamiltonian structure 

B2 = -pDl + D,oT + T-D„ 

for the equation SKdv endowes the Fourier coefficients t^ of the energy- 
momentum component (see Eq. (jSH))) 






k& 






for the Toda equation ()55|) with the structure of the Virasoro algebra 
([SH]) with the central charge c = — 1/3. We have 

27ri [in, U] = 2(?2 - m) !„+„ - /? ■ (n^ - n) (5„+m,o, 

We emphasize that the known property jHl EH] of the integral T is 
adapted to the case of the Toda equation ()55j) associated with a non- 
degenerate symmetrizable matrix K that is not necessarily the Cartan 
matrix of a semisimple Lie algebra. The central charge c is expressed 
by using the constant [3 = ^^ OjA* that depends on the symmetrizator 
a. See Remark El on page IHHl for a similar situation. 
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In the paper J52], we considered a class of generalizations of the 
Virasoro algebra, see Eq. ()132j) . such that the determining relations for 
those generalizations are given by an iV-ary skew-symmetric bracket 
and A^ > 2. 



7.2. The analogs of the modified Korteweg— de Vries equation. 

In this subsection, we investigate the relation between the potential and 
nonpotential modified Korteweg-de Vries equations that correspond to 
the symmetry subalgebra 2t C sym /^xoda C sym £^Toda of the Euler type 
Toda equations (|55|) . We analyse the properties of the Miura transfor- 
mation T = T{'&, '&x) that maps solutions of the modified Korteweg-de 
Vries equation p27j) to solutions of Eq. ^T^ and thence we estab- 
lish the invariant nature of the operator D, see Eq. fl(i2|l on page EHl 
that appeared in description of the structure of symmetries for the 
Toda equations. Also, we construct the Hamiltonian structures for the 
hierarchy 21 and prove that modified equations ()93|) and ()127|) share 
the Hamiltonians [hk dx] with the Korteweg-de Vries equations ()91|) 
and dni. 

We recall that the dependent variable {} introduced in Eq. ()128|1 
satisfies the modified Korteweg-de Vries equation ()127j) and its higher 
analogs SmKdY{k)- By Corollary IHHl on pageEHl all right-hand sides in 
the equations £^mKdv(A:) are local with respect to {}. 

Similarly to Sec. 17.11 we establish identifications between the higher 
symmetries of the potential equation ^pmKdv and the generating sec- 
tions of conservation laws for the equation S^Kdv, and vice versa: 



Hk[u 




R 



^mKdV 

mKdV 



pmKdV 



^t.+. = ^rr = t 



^fc[^] — ^^ '"tfc = 0k = i', 



mKdV 



id 




pmKdV 



Ai=A- 



k+1 



RmKdV— RpniKdV 



^t, = ^l 



mKdV 



^k 



mKdV 



'Hu-m-^n,,_,=0k-i=^T-l 



KpmKdV-RmKdV 



R„ 



K-. = VkT = i^k-i 



R„ 



A2 
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The Lie algebras 21 and 03 are commutative, see Lemma EHl and 
Theorem 1^ Therefore, Jacobi's identities ()7b|) are vahd for the oper- 
ators Ai 2 and Bi2 that are either induced by relations ()12H) for the 
subscript 1 or obtained from the recursions p23p by using the Ma- 
gri scheme ()125j) for the subscript 2. By construction, the recursion 
operators -Rp(m)Kdv admit the standard decompositions 

RpraKdY = ^20^^1=00 D^^ O i^, 
-RpKdV = B2 O Bi , 

where the skew-symmetric factors Ai 2 and -81^2 are 

A, = k-'oD-\ (Br=D~\ 

A2 = noD-'o n*, \B2 = -f3D, + si- D-^ + D~' o s,. 

We denote by i^ the linearization of functional (jKH|) with respect to 
the dependent variables u and by £f. its linearization with respect to 
i). The operator £^ is given in Eq. ()74|1 on page ESI and we have 

Lemma 50. The relation D* 



Proof. The verification of the first statement is straightforward. Ex- 
pressing Ux via i), Ux = K~'^'&, we write down functional (J58|) in terms 
oi ^: we have 

T=U2 '^'"'^'^'" - 5^ ^' ■ ^^ 



(••• 


.1: 

a 


dT 


■„...). 


n* = 


- «y 


holds. 


Moreover, 


i^ 


= n 


•oil 





n 

Lm 



' XI 



where k ^ = ||/^'™|| and the identity 

r 

is used. Consequently, we obtain 

The second statement follows from the definition of the linearization 
on page |H1 □ 

Several decompositions follow from this lemma. For example, we 
obtain 

iT = n*o A^^ = n* oAi = n*oDxok = n* o £«. 

The recursion operator -RpKdv is factorized to the product 

-RpKdv(0fc-i) = 9n(0fc_i)(s) = is{0{(j)k-i)) = 
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of the Hamiltonian operators glued together by the operators D and 
D*. 

Proposition 51 dSUED])- Each Noether symmetry 

ipc = n o EriQix, T)) e sym^Toda 

of the Toda equation associated with a nondegenerate sym,m,etrizable 
m,atrix K is Hamiltonian with respect to the Hamiltonian structure 
Ai = K,~^ ■ D~^ and the Hamiltonian Ti = [Q{x, T)]: 

The proof of Proposition |^ follows from the definition of the Euler 
operator, E{H dx) = ^^(1), and Lemma IHHl 

Remark 18. The propetry of the Noether point symmetries 

4 = noET{T-f{x)dx) 

of the Toda equation to be Hamiltonian, 

ipl = AioE^{T-f{x)dx), 

was established in |77]. One easily verifies that equation (j^ is also 
Hamiltonian with respect to the cohomology class Hq, see diagram ()131|) 
on page EH 

ipi = Ai oEu{hodx). 

Now we extend Remark ^J onto the whole hierarchy 21. Here we 
formulate the most remarkable relation between the hierarchy 21 for 
the potential modified Korteweg-de Vries equation ()9H|1 and the hier- 
archy 03 for scalar equation (jH!)). 

Theorem 52 ([SDinn])- The Hamiltonian [h^dx] for the kth higher 
Korteweg-de Vries equation is the Hamiltonian for the Noether sym- 
metry (y9fe G 21 for any integer k > 0. 

This theorem follows from Proposition 15 II and Remark ITHl on pagelTKl 
that assign the Noether symmetries 99^ G 21 to the conserved densi- 
ties /ifc. 

Remark 19. The Lagrangian representation of the Hamiltonian equa- 
tions S(^k) within the hierarchy 21 is 

£{k) = |AioE„([-i^fi:iXM^^ -hk^i] dx) =o|, k> 0. 

id 

Here we complete the description of the Magri scheme for the equa- 
tion (jHSI)- The Toda equation ()118j) is the first nonlocal term of the 
symmetry tower below: 
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r: 



■pmKdV 



R, 



■pmKdV, 



D 



ipKdV 






St, 



hk 



i^k 



R. 



pmKdV 



R, 



pmKd\> 



Ut, = (fk 



D 



B2 



^pKdV 



-RpKdv V'fc-i ^ hk-1 



Eu 



hk-1 —^ if. 



k-l 



A, 



Bi 



Si,_i = (pk-1 



B2 



R, 



pmKdV, 



Ut,_, = 0k+l ^- Si,_, 



R. 



pKdV 



^pKdV i'k-2 "—^ h 



k-2- 



Bi 
'>k-2 



Finally, we have the following assertion 
Proposition 53 ([^IHO]). (1) The substitution 

T = T(^, ^^.) : £^mKdv '-^ ^Kdv (133) 

is the Miura transformation between the higher equations 

^mKdV(fc) = {^tfc = -Dx ■ k{(fk)} 

and 

^Kdv(fc) = {Ttk = -Ds(0fc)}. 

(2) The operator 

maps the Lie algebra syni£^mKdv 3 f of symmetries of the equa- 
tion £^mKdv iiT'to the symmetry Lie algebra sym SkcW for the Kor- 
teweg-de Vries equation (jHH), see Example fTHl 

(3) The operator D = D** maps the generating sections of conser- 
vation laws which are dual to the symmetries, 

D : (^fc = ^Ahk dx) f-> v^fe G 21, 

in the opposite direction. 
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Remark 20. From diagrams ()1U4|) and ()13H) we see that the successive 
apphcation of the mappings D* = £|, and D = D** does not preserve 
the index of a higher symmetry: 



i'k 



,^mKdV 

Yk 



W 



i'k 



$t > 



Ar 



0k = ^, 



mKdV 



n=n* 



l-pKdV 



(134) 



>k-l 



and the recursion operator -RpKdv measures the difference in their ac- 
tion. 
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Part II. Group-theoretic properties of the mathematical physics 
equations: methods and apphcations 

Within Part II of the present paper, we illustrate the cohomologi- 
cal concepts and algorithms in the geometry of differential equations 
([iniESESlinil)- Namely, we apply these methods to the study of the 
properties of the dispersionless Toda equation, the multi-component 
nonlinear Schrodinger equation, the Liouville equation, and related 
systems. Invariant solutions, Noether's symmetries, local and nonlo- 
cal conservation laws, weakly nonlocal recursion operators, parametric 
families of Backlund (auto)transformations, and zero-curvature repre- 
sentations are the structures, which are obtained for these equations of 
mathematical physics. 

Chapter 3. Symmetries, solutions, and conservation laws for 
differential equations 

In this chapter, we consider two examples of applications of the PDE 
geometry methods; namely, we analyse the geometric properties of the 
dispersionless Toda equation (jnOlEE]) and obtain the structures (fHKj) 
for the multi-component analog of the nonlinear Schrodinger equation, 
which is related with the former one (see [T^l. 

8. Nonlinear Schrodinger equation 

In this section, we study the properties of the m-component analogs 
of the nonlinear Schrodinger equation (123111122]) 

^t = i^.. + ifm)^, (135) 

where ^ is an m-component vector {m > 1), i = \^—l, and / G C^(R) 
is a smooth function. 

The scalar {m = 1) nonlinear Schrodinger equation describes the 
propagation of light pulses or wave packets in the linear dissipation 
and nonlinear self -focusing media, e.g., in optical fibers, nonlinear 
cristals and gases, multicomponent Bose-Einstein condensates at zero 
temperature, etc. The collective evolution of a complex of space- 
incoherent solitons in a nonlinear media with a Kerr-like nonlinearity 
is described by system ()135|1 of the nonhnear Schrodinger equations, 
where \& = *(\l/^, . . ., ^1/™) are the amplitudes of the light wave, 

•m 

is the density of energy, t is the coordinate along the direction of the 
wave propagation, and x is the coordinate along the front of the wave. 
Suppose / = id, then this equation is known (jHE!) to admit a com- 
mutative bi-Hamiltonian hierarchy of higher symmetries and an infinite 
number of conserved densities in involution. Still, for an arbitrary / 
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this is not so. In the sequel, we obtain the symmetry algebra of the 
latter equation in the physically actual case of a homogeneous func- 
tion / of weight A; also, we desribe the set of rn^ conserved currents, 
which generalize the previously known m separate conservation laws 
of energy of an ith mode and the conservation law of total momen- 
tum of the system. Those conservation laws correspond to two types 
of Hamiltonian symmetries of the Schrodinger equation, the scaling 
transformation and the translation. 

So, consider the m-component analog of the nonlinear Schrodinger 
equation, the multi-soliton complex (P) 

-. -I -T -, - (136) 

F^ = ^^ + i^^^ + if (I) ■ ^'^ = 0, 1 < A; < m, J = ^^. 

This equation is Hamiltonian for an arbitrary function /: 

'^\ _ / 1\ f6nNLs/S^ 



where the density .^nls of the Hamiltonian TYnls = [-f^NLS dx] is 



We also note that -f^NLS is a conserved density for Eq. ()136|) , quite sim- 
ilarly to Eq. |TT| for the dispersionless Toda equation or to Eq. ()119|1 
on page IHHl for the Toda equation ()55|1 . From the Hamiltonian rep- 
resentation ()137|1 we deduce that the complex conjugate variables ^ 
are the canonically conjugate variables (the momenta) for the dynam- 
ical coordinates \Ef; in Sec. El we investigated a generalization of this 
situation. 

Equation ()135|) is known to possess the following property (see [T] 
and references therein): in addition to the conservation of the total 
energy ^^, there are m separately conserved densities 

Qi = -^'- ¥, DtiQi) = 0. (138) 

Physically speaking, there is no energy transfer between the modes ^*. 
Still, this observation about the properties of the nonlinear Schrodinger 
equation is uncomplete, since these m integrals of motion p38|) are only 
particular instances within the set of m^ currents 



ry.^. = ^^^i dx + i (^; ^^' - ¥ ^{.) dt, dhVij = 0. (139) 



These currents are conserved on Eq. ()135|1 with an arbitrary nonlin- 
earity /(X). They reflect the correlation between the intensities of the 
light modes with different indexes i and j, 1 < i, j < m. The set of 
conservation laws ()139|1 remained unnoticed in the papers [H 102] , and 
the arising conservation restrictions seem to have not been taken into 
consideration by the authors when making computer experiments. 
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The generating sections of the conservation laws [rjij] defined in 
Eq. dnni) are 

where 

and 

^fe) = '^(ij) = f°^ ^' ^ ^' ^" ^ ^'• 
The canonical Hamiltonian structure 



Ti 



1 
-1 



assigns the point symmetries 

to these sections. Suppose i = j, then we get the gauge symmetry (|nHl) 



Axi^iii)) : ^' ^ exp(A)^\ ¥ ^ exp{-\)¥. (140) 

The Hamiltonians of the symmetries (^(m) are the conserved densities 
Qi in Eq. (fTT^ . 

In the paper [98j, the case /(X) = X of cubic nonlinearity in Eq. pH5|l 
was analysed. Then, the cubic ??i-component equation ()135|) admits the 
recursion operator 



R 



NLS 



'D, 
D, 



-^ 



+ 2 ^ ■^."°(^'^)■ 



o5Z^fe■)■^-'°^fe■)• 



«J 



This operator generates an infinite sequence of the commuting local 
higher symmetries of Eq. p35p . starting with scaling symmetry 



Qidx 






'^^NLS 



TIa 



E* 






^1 



V'a 



-Rt 



V^3 



; T -^NLS 



Ti 



ri 



^ ^NLS ,t) ^NLS ,t, -RnLS -^ -RnLS -» -RnLS 
V^{n) ^ y^x ^ ^t ^ V^3 ^ ¥^4 ^ 

(141) 
This diagram is an example of a sequence of Hamiltonian symmetries 
such that there are no Hamiltonians for a half of these symmetries, 
02k+i- we have ^a/c+i ^ imE^. 
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Remark 21. We have obtained the formulae for the generating sections 
ijj(^ij) of conservation laws ()139jl and the Hamiltonian symmetries 0{ij)- 
Now we conclude that the recursion operator -Rnls is weakly nonlocal 
([HT]). i.e., this operator is decomposable to the form 

i? = local part + > ipa ° D~^ o ip^, (142) 

where ^a are symmetries and t/'q, are generating sections of conservation 
laws. We recall that in the preceding chapter we used the recursion 
operators' decompositions similar to p42j) . and then we generated the 
commuting symmetry sequences for the Toda equations and established 
nontrivial properties of these sequences. 

We finally note that the symmetries (^ that are independent on x 
and t (only such symmetries were considered in [^H]) do not exhaust 
the whole algebra of classical symmetries for the nonlinear Schrodinger 
equation ()135p . Moreover, the resulting symmetry algebra of this equa- 
tion is noncommutative. 

Example 21. Suppose there are no restrictions for the nonlinearity func- 
tion / (X) in Eq. ()13fi|l . Then the point symmetry algebra of the latter 
equation is generated by the fields whose sections are 

Further on, suppose additionally that / is homogeneous, 

/(AJ) = A^-/(J). 

Then equation ()135|) acquires the scaling symmetry 

ip^\ _ f2^^ + Ax% + 2At^ 
(^7 ~ \^2^^ + A X ^^. + 2 A t ^j 

9. The dispersionless Toda equation 

The dispersionless Toda equation is an analog of Eq. (fS^ with the 
continuous variation of the index j that enumerates the dependent 
variables uK In this section, the classical symmetry algebra of the 
dispersionless Toda equation is computed and 5 classes of conservation 
laws for the latter equation are reconstructed; also, some questions of 
the Lagrangian formalism with higher-order derivatives are discussed. 

Consider the hyperbolic Toda equations 

u^y = exp{Ku), u'^y = K ■ exp(tt'), (143) 

associated with the type Ar-i Lie algebras with the Cartan matrices 
K, introduce an additional independent variable 2 G M, and extend the 
values of the discrete index j G [1, t], which enumerates the dependent 
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variables u\ onto the whole line M. Also, let the rank r tend to infinity: 
r — i> oo. Then for any section u of the bundle vr we set 

u^ = u{x,y,z)\^^.^, 

where e is the cell step. 

As r — i> oo and e -^ +0, the continuous limit equation, which is ref- 
ered as either the dispersionless Toda equation, or the heavenly equa- 
tion, or the Boyer-Finley equation (|II]), also appeares in many mod- 
els, e.g., in gravity ([HZIIII]) as a reduction of the antiselfdual vacuum 
Einstein equation (ASDVEE). The heavenly equation ^heav exists due 
to the special structure (see Eq. (HH))) of the Cartan matrices K = \\kij\\ 
of the type Aj. algebras. As a matter of fact, instead of the Cartan 
(r X r)-matrices K one should consider (|HZIIHH1) the Cartan operator 
K = —D^. Then the scalar equations, which originate from Eq. ()143p . 
are (0111]) 

S = {F = u^y- exp{-u,,) = 0} (144a) 

u',y = -Dloexp{u'y, (144b) 

the Cartan operator K defines the sign "— " in the one-dimensional 
equations 

Urr = exp(-M2^), u'^^ = -DI o exp(M') (145) 

that appeared in the Introduction. 

Following the papers J87| 116] , we consider the limit 

lim lim Lxoda = D^^ {\uxUyzz - exp(-M^2)) 

e— >+0 1 — ►oo 

of the Lagrangian density /^xoda, defined in Eq. (]K?)]l . as r — > oo and 
e — i> +0. Then the Lagrangian Z^Toda itself is mapped to the functional 



C = dxdy Ldz 

with the density 

L = -| u^zUyz - exp{-Uzz) 

that depends on the second derivatives of the sections u = u{x,y,z). 
Applying the Euler operator E„ to £, we obtain the equation 

^heav = {-^heav = U^yzz - D^ O exp{-Uzz) = O} ; (146) 

we see that double integrating in z maps Eq. ()146|) to Eq. ()144a|l . while 
the substitution u' = —Uzz niaps it to Eq. ()144b]) . 

9.1. Symmetries and exact solutions. Computing the symmetries 
if e Sym £ of Eq. ()144aj) by using the analytic transformations software 
Jet ([Z2]), we get the following solution: 
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Proposition 54. The point symmetries 

ip{x,y,z,u,Ux,Uy,Uz) 
of Eq. fll44a|) . which solve the determining equation 
Dxy{^) + exp(-M2^) ■ Dl{ip) = 0, 
are 

<^iM = ("y - l^^ Dyj g{y), (147a) 

(147b) 

(147c) 

<^4[g]=g(2/)2;, (147d) 

^5[r]=r{y), (147e) 



^ilf] 



^2 

ip^iq] = q{x) z, 
cp^[r] = r(x), 



u^- \z^ D^j f{x) 



-^ZUz + U 



2^ ' 



where f , q, and r are arbitrary smooth functions of x and g, q, a nd f 
are functions of y. The commutation rules for symmetries ()147|1 are 
given in the skew -symmetric table below: 





Vi f 


^1 9 


(P2 


V3 


ip4 q 


(^4 9 


(P5 r 


>P5 r 


mif] 











^4-f'] 


^4-fq'] 





^^o[~fr'] 





viig] 










^4-9'] 





vil-gq] 





V'si-gr'] 


V2 









ip3+ip42] 


'P'iiq] 


'fM 


^5[2r] 


V5[2f] 


<^3 











'P5[-q] 


<^5 [-q] 








f4 q 






















(f4 q 





















Vb r 




















>f5 r 




















Remark 22. 1. The operator D = u^. — |z^ D^ in Eq. ()147ap is an ana- 
log of the operator D, which was introduced in Eq. (jH^ on pagelHHl 
2. The symmetries of Eq. ()144bjl that correspond to ()147aj) - ()147c|l 
were obtained in [2j; the symmetries (p^, ip^ G kerD^ of the equation 
S demonstrate the (inessential) distinction between the geometries of 
Eq. ()144a|l and Eq. (jl44b|l . this distinction is provided by the substi- 
tution u' = —Uzz- 

Remark 23. In fact, equations p44|l and ()14(ij) are huge inspite of an 
apparent compactness of the notation. Indeed, one can estimate the 
number of A;th order internal coordinates for these equations or the total 
number of nontrivial relations that define the kth. prolongation S^J^^^- 
Therefore the calculation of the higher symmetry algebra sym ^^^^v f*^^ 
Eq. (fri6|l is really difficult. 

The passage from difference-differential equations pH|l associated 
with Cartan's matrix pH|l to their dispersionless limit p44aj) provides 
the additional constraint 



"•zzz 



0. 
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The latter equation must be taken into consideration when describ- 
ing the Lie algebra sym f^^^v composed by all symmetries of the hmit 
dispersionless Toda equation. 

9.1.1. Constructing exact solutions to the dispersionless Toda equation. 
By using the structure of the point symmetry Lie algebra SymS for 
Eq. ()144ap and applying various geometric methods (^3), the problem 
of constructing exact solutions to the heavenly equation was considered 
in j2] and [TTj. In the former paper, invariant solutions of the equa- 
tion Uxy = ±(exp(M))^^ were obtained; they correspond to Eqs. ()152|) 
and ()153j) of the present article. A class of solutions that are non- 
invariant with respect to the whole Lie algebra SjmS was pointed out 
in [21]. The approach of this paper is wider than in |2j in the follow- 
ing sence: any time an equation under consideration being reduced to 
an auxiliary equation for a function that depends on fewer arguments 
than the initial one, we construct the point symmetry Lie algebra for 
the new equation and obtain its invariant solutions that are parame- 
terized by arbitrary functions. Therefore, what we get is much more 
than a list of particular solutions. Also, we note that the method of 
solving Eq. ()144a|l together with the constraint ipi = is similar to the 
scheme (|ZI]) used to construct non-invariant solutions of differential 
equations. Indeed, we treat the variable z in the constraints (pi = as 
a formal parameter in order to write down these additional equations 
(fii = in total differentials. 

The symmetry ipi G SymS. Consider the generator 

^1 = ^i[f] + ^ilg] =u^f{x) - \z^ f {x) + Uy g{y) - \z^ g' {y) 

of an infinitesimal conformal symmetry of Eq. ()144a|l : this generator 
depends on two arbitrary smooth functions / and g. Now we represent 
the invariance condition y^i = in the characteristic form and obtain 
the first integral 

/x ny 

dx/f{x) - / dy/g{y) 

that is independent on the variable u. In order to construct another 
integral C2 we treat the base coordinate z as a parameter: 

^z"^ log fix) + \z^ \ogg{y) -u = C2{z). 

Then, a solution u{x,y,z) to the equation ipi = is defined by the 
condition 

U{t,C2{z))=0, 

where 11 is an arbitrary (for a while) function of two "integrals" , one 
of which depends on z. Solving this relation with respect to u, we get 

u = W\ogf{x)g{y) + ^t,z). (148) 
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Moreover, if m is a solution to Eq. ()144a|) . then it is necessary that the 
function $ satisfies the one-dimensional equation 



$ii + exp(-$^^) = 0. 

The reduced equation is nothing else than the one-dimensional analog 
of Eq. p45|) . In a similar situation, the hyperbolic Liouville equation is 
reduced to solving its one-dimensional analog when one constructs so- 
lutions that are invariant with respect to the conformal symmetries of 
the initial Liouville equation. We emphasize that the second "integral" 
C2{z) of the equation ipi = was obtained by using a nontrivial inter- 
pretation of the coordinate 2; as a parameter which is auxiliary with 
respect to the independent variables x and y. This approach seems 
to have not been used in j2|; in the sequel, we apply it to Eq. ()145l|) 
again. We use the Jet ([Z2]) software and compute the generators of the 
point symmetry Lie algebra of Eq. () 145111 ; then, we obtain the invariant 
solutions of this equation. The result is given in 

Lemma 55. The Lie algebra of classical symmetries of Eq. 1)145111 is 
generated by eight evolutionary vector fields whose generating functions 
are 

(j)l = t$t - Z\ 02 = $t, 03 = Z<l>, + Z^ - 2$, 04 = $^, 

05 = Zt, 06 = t, 07 = Z, 08 = 1. 

Further on, we point out the famihes of solutions to Eq. ()145l|l and 
assign the solution class ()148|1 for the dispersionless Toda equation (see 
Eq. ()144a|) ) to each of them. 

The symmetry 0i. Consider the equation 0i = 0; again, we require 
the variable z to be a formal parameter, then we have z'^ log |t| — $ = 
C{z). Hence we express the solution $, substitute it into Eq. ^451)) . 
and solve the latter with respect to C; we finally get 

$(t, z) = z^ log |t| - \z^ log \z\ - \z^ + CxZ + C2, 

where C\ and C2 are arbitrary constants. 

The symmetry 02 . One easily checks that the real 02-invariant solu- 
tions of Eq. p45l|l do not exist at all, since the equation exp(— $22) = 
is insoluble over M. Still, we consider a more general case: namely, we 
seek solutions to Eq. ()145l|l that are invariant with respect to the linear 
combination 

0(a:fe) = 02 + (a : &)04 

of the symmetries, where {a : b) E MF'^, i.e., we seek the propagating 
wave solutions. Substituting the function $(2; — (a : b) t) into Eq. ()145H) . 
we get the equation 

(a:6)2^ = -exp(-^). (149) 
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Here by \l/ we denote the second derivative $" of $ with respect to its 
argument 

w = z — {a : b)t. 
Suppose the condition 

—y/e < {a : b) < \^ 

is satisfied, then the obstacle for the low-velocity 0(a;fe)-invariant so- 
lutions to exist is clear from Eq. ()149p . Indeed, there is the critical 
(minimal) velocity \a : b\ = y/e and the wave solution 

$ = -^{z ± ^tY + a{z ± ^t) +p, 

where a, P eM.. Quite remarkably, if the velocity -y/e < |a : 6| < cxd is 
greater then the minimum velocity, then there are distinct roots \E'i^2 
of Eq. p49|) and Eq. ()145l|) admits two solutions 

$ = 1^1 2 ■ ti'^ + au; + /3, w = z-{a:b)t, a,/3eM, (150) 

similtaneously. This pair of roots \I'i 2 appeares after the obvious bi- 
furcation: no solutions of Eq. ()149p and (jl45l|) are assigned to a low 
inclination (a : 6)^ of the straight line y = (a : 6)^ \E' on the 0\E'|/-plane; 
then, at the inclination (a : 6)^ = e this straight line is tangent to 
the graph of exponent y = — exp(— \E') at the point (—1, — e), while 
at a greater incline e < (a : 6)^ < 00 the point of tangency is spht 
to the pair of distinct intersection points (^1^1,2, (a : &)^^i,2), where 
— oo<\E'i<— 1<^2<0. Each admissible value of \E' determines the 
second derivative of a polynomial solution p50|) of Eq. ()145l|) . A unique 
root \E' = and the 04-invariant solution of Eq. ()145l|) (see below) is 
assigned to the point a : b = 00. 

The symmetry 03. First, we solve the auxiliary ordinary differential 
equation: 

xy\x) — 2y = —x^ =^ y{x) = {l — \ log x^^ ■ x^. (151) 

We shall use formulae ()15H) twice, in order to construct the solutions 
of Eq. ()145l|l that are invariant with respect to the symmetry ^3 and 
to get the (y92-invariant solutions of the initial equation ()144aj) . In these 
two cases, the function 7 depends on two distinct sets of variables: the 
first case is 7 = 7(t) and the second case is 7 = 7(x, y). 
Substituting the expression 



$ = (7(t)-ilog2;2 



z 



2 



that appeared in Eq. ()15H) into Eq. ()145l|l and paying attention to the 
sign of the constant of integration, we get the following expressions for 



SYMMETRIES AND CONSERVATION LAWS IN INTEGRABLE MODELS 91 
the function 'jit): 



7i = § - log 



■\/esinh2 artanh exp(±-\/e(t — to)) 



72 = | + log[±t-to]^ 



73 = I - log 



2 



v/ecoshlogtan(v^(±t — to)) 



where e > and to G M. Now we return to formula p48j) and assign 
the solution class 

u = \z^ log f{x)g{y) + (7,(t) - \ log z'') ■ z\ % = 1, 2, 3, (152) 

of the dispersionless Toda equation ()144a|l to each of the functions 
$(t, z). 

The symmetry 04. Consider the constraint ^z = 0; then, a solution 
to the equation ^tt = 1 is a polynomial 

$(t) = 1^2 + Cit + C2, 

where Ci and C2 are constants of integration; the corresponding solu- 
tions of the dispersionless Toda equation ()144a|) are again defined in 
Eq. 



The symmetry ip2 G Sym S . Now we construct the y92-invariant solu- 
tions of Eq. ()144ap . First, we consider the auxiliary ordinary differential 
equation in Eq. (JISII) : then, we obtain the hyperbolic scal^-Liouville 
equation (see Eq. (jH^ on page lHUjl . Now, subsitute the expression (7(x, 
y) — |logz^) ■ z"^ for u in the equation u^y = exp(— m^^) and make the 
transformation 

X = x exp(3/2), y = y exp(3/2). 

Hence we get the Liouville equation 

-fxy = exp(-27) 

whose solutions (1001) ^^^ easily obtained from formula (|H7jl on pageHHl 
by a trivial coordinates change; we have 



7(A', 3^) = -I log[-fiX)g'iy){Q{[f{X) + g{y)Y) 
where the mapping Q is one of sin, id, or sinh, whence we finally obtain 



^2 }Q([f(e'/^x)+gie^/^y)r) 
u{x, y,z) = - log _,2;.(e3/%),^(e3/2y) ' ' ^''^^ 

This class of solutions of Eq. ()144a|l is in fact present in paper j2] , and 
their physical interpretation is known (ibid.): these expressions pro- 
vide the instanton solutions (^Hj) of the antiselfdual vacuum Einstein 
equations (ASDVEE). 
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The symmetry (f^ G Sym£^. In order to obtain the 993-invariant solu- 
tion of Eq. ()144a|l . we solve the equations Uz = and u^y = exp{—Uzz) 
together. The solution is 

uix, y, z) = xy + f{x) + g{y), 

where / and g are arbitrary functions. 

The symmetries ip4 and ip^ of Eq. ()144ap depend neither on the 
unknow function u nor its derivatives. Therefore, if the equation 
Uxy = exp{—Uzz) is overdetermined by the condition ip^ = or ip^ = 0, 
then the problem of search for exact solutions to the dispersionless 
Toda equation is not simplified. 

9.2. The Noether symmetries and conservation laws. In this 
subsection, we consider the problem of constructing the conservation 
laws for the dispersionless Toda equations ()144j) and ()14fi|l . First, we 
discuss the general correlation between an Euler equation (fTH|) and the 
fixed set of conservation laws which reflect the conservation properties 
of the energy-momentum tensor. 

9.2.1. On the Lagrangian formalism involving higher-order derivatives. 
Suppose that C = j L da; is a Lagrangian with the second order deriva- 
tives and g^v is the corresponding metric; then we can construct a con- 
servative analog of the energy-momentum tensor T'^^ that satisfies the 
equation 

Y,D,{T^n = ^. D,^D,\^. (154) 



M 



The equations of motion £ are 






dL \ dL 



dui, / du^ 



(155) 



and the tensor T'^^ coincides with its classical definition ([Hj) if the 
higher-order term 

dL 



is trivial. Again, we use the notation m!^ = D^t^iv}) and u]^.^^ = D^t^ ■ 
Dx-^iu"). We raise and lower the indexes ; /x, ; i^ by using the metric 
g^j^y. By a straighforward substitution into Eq. p54j) one checks that 
the tensor 

satisfies the required restrictions ([M"), and therefore condition ()154|) 
provides conservation laws for system ()155p . Still, we note an important 
feature of Eq. ()146|) : as r ^ 00 and e -^ +0 in the Toda equations, 
the metric gf^^ becomes degenerate and the Lagrangian £ looses its 
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covariance. Indeed, it is impossible to find a nondegenerate metric g^^ 
such that 



L ~ iu;^'"M'^;^ + exp(-up. 



For the same reason we cannot define the energy-momentum tensor 
T^^ for the equation £^heav as the variation of the Lagrangian £heav 
with respect to the metric g^^ since the nondegenerate metric does not 
exist. 



9.2.2. Constructing the conservation laws. From the preceding section 
it is clear that the search of conservation laws for Eq. ()146|) is less triv- 
ial than finding the symmetry Lie fields by their generating functions 
and requires additional reasonings. In the sequel, we find the symme- 
tries sym£^heav of the dispersionless Toda equation, select its Noether 
symmetries sym £heav in a fixed coordinate system, and reconstruct the 
conservation laws by using the homotopy method, which is described 
in [El EH and ^. 



Lemma 56 dSOj)- The generating functions ()147a|l . ()147c|l - ()147ej) are 
the Noether symmetries of the equation £heav, while ()147bjl is not. 
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Proposition 57 ([50]). The conserved currents rji G A^(£^heav); (as- 
signed to the Noether symmetries (pi of Eq. p46|) . are 



11 Z I 

- Y^/'(^) ^y~ 12 f^^^^^^^y^ + 6'^'*'^'' "^^ ^ ^f{x)uxUyzz 



v2 /•! 

^yzz 



Z 



+ f{x) / tuu^zzz exp{-tuzz) dt\ dy A dz 

+ |-g/'(a;)MMx^^- -/(x)mm^^^^ + -/"(x)m + — /'(x)m^m^^ 
1 z 1 1 

+ TTt/(^) ^zUxxz - -^f"{z) U^ - —f'{x) UxU^^ - —f{x) UxxUzz 

VI o z4 24 

^2 , 1 1 1 

-f{x) UxyUxz + -/(a;) MxMa;j/2 " T^/'(^) ^^"j/^ ~ I'^''^^^ ""^^ 

Z Z^ 1 1 

fi ^'^' ^y ^r ^"^^ "^xyz + TT^/v'^j Uyllxxz T/v-^J UUxxyz 

1 1 Z 1 

T TTTjl^a^j U^Uxxy T7)/\*^J '^xx'^yz + Ti'/ V'^J ^xj/ + Tq/ l"^/ '^y'^xz 

+ T^/'(^) UzUxy+-;T;f"i^) Uyz+f{x) ■ / [-tMM^^^^-tM^^n^^+tw^n^ 

^0 



12 ' ' ^ 12 

^tU^Uxzz+t^UUxzzUzzz-t^UxUzzUzzz I expl-tw^.^^) dt + f {x) ■ I \ -U-^ 



^] exp(-tM;,;2) dt + /'(x) • / [-1 

Jo 



....4..-....4....Je.p(-..J.}d..d„ 



^73 = { ^u^Uy^^ - -Uy^Uzz + -^rrUyUzzz - -sUUyzzz > dy A dz 

J — -- — -- \ r] r] 

' 1 r\A '-^Z'-^XZZ '^^XZ^ZZ ' r\A '"'X"'ZZZ '^"'"'XZZZ f "Z /\ (IX 

11 1 1 1 1 

+ <! ■^UUxyzz + -UzUxyz--UxyUzz--UxzUyz + -UxUyzz + -UyUxzz 

+ Uzzexp{-Uzz) + exp{-Uzz) + UzUzzzexp{-Uzz) \ dx A dy, 
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r]4 = I --q{x) Uyz + -q{x) Uy-^z \ dy ^dz 



,.-,-., ,,-.-. ,,-.-.,i- 



+ \ o9(^) "^J/2 ~ «9(^) ^^J/ ~ «^'(^) "j/2 + «^'(^) S 



^g(a;) M^^^ exp(-M2^) + g(x) exp(-M22) ^ dx A dy 



rix) , , f r'ix) rix) 

j^w,..dyAdz+<j-^^^- ' 

r'(x) r(x) 



?75 = — ^ Uyzz dy A dz + <j -^^ Uzz + — ^ «x^^ !> d2; A dx 



+ { :^ Uyz + ^— u^yz + r{x) Uzzz(iw{-Uzz) \ dx Ady. 



Any divergence dh{Vi) equals on the equation Shea,v', calculating the 
divergences does not require precise evaluation of the integrals in the 
homotopy parameter t in rational functions , because differentiating with 
respect to x or z under the integral sign is allowed. 

Proving Lemma EEl and Proposition EH is based on Lemma EHl and 
involves the general scheme of reconstruction of conservation laws by 
their generating sections. 

Method of reconstruction of conserved currents. Now we describe a 
method ([H^) ^^^ constructing the differential (n — l)-forms that are 
exact on the equation £ = {F = 0}; we follow the paper [45J. 

Consider the function / = u/t and the evolutionary vector field 3/ 
whose flow is 

Ar-. {X\ui) ^ {x\Tui), 

then we have 

^^ = Al{3,{u)) = Al{L{f)) (156) 

for any differential form u. Let this form u be 

uj={F,ip) = Fipdx^A...A c/x" 

and suppose that rj is the desired current which is assigned to the 
generating function ip on the equation {F = 0}: dh{ri) = V(F) and 
ip = V*(l). We note that the right-hand side in Eq. ()156p contains the 
term 

C(/) = (C(l),/) + d.G(Co/), 

where the first summand is cohomological to since the image of the 
Euler operator is trivial if the Lagrangian density is a total divergence, 

C(l) = ^(F,^>(1) = E((F,^)) = E(d,77) = 0. 
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Next, the mapping G: CDiff(jF, A") -^ A"~^, which is restricted onto 
the equation £, is defined by the rule 

o- |cr|>0 jea ' 

where 

uj^_j) = (-1)^+^ dx^ A . . . da? . . . A dx" 

and cr — Ij is the resuh of excluding an index j from a multiindex 
a = (cTi, . . ., an). 

Now we integrate Eq. ()156|) in r from to 1 and get 



"^ Aliu) dr = Aliuj) - A*iuj) = {F, ^) - A* {{F, ^)) 



d'^ 



AX (C(/)) dT = dH f a; (G(C of))dT = dhV, 



'0 JO 

whence we obtain the required formula for the current rj. 

Proof of Lemma EHl and Proposition EZl First we note that symmetries 
fjl47|) of Eq. ()144a|) are also symmetries of the Euler equation ()146p : 

Sym£^ C sym£^heav5 because the relation ^F^cav — ^1° ^P holds due to 
Lemma [7| on page 1181 By using Lemma EHl we check that symmetries 
p47a|) and ()147c|) - ()147e|) of Eq. ()146p are Noether, while for symmetry 
()147b|) we have 

^F,_(</'2) + ^;,(i^heav) = ^U^yzz ' |l),'(exp(-M,,)) ^ 0. 

From Theoreminiit follows that the Noether symmetries of Eq. ()146|) 
are in one-to-one correspondence with the generating sections of its 
conservation laws. By Corollary El on page IHEl this correspondence 
for the equation £^heav is provided by the identity mapping. Following 
the reasonings scheme above, to each of four Noether symmetries ipi 
we assign the conserved current rji] the arising expressions are rather 
large and therefore omitted. Nevertheless, the result is correct, for the 
conditions d/,(?7j) = may be checked straightforwardly. D 

Remark 24. Suppose that the arbitrary function / in the conserved 
current ?7i is / = L Then rji is a countinuous in z analog of compo- 
nent ()58|1 of the energy-momentum tensor B = T dx+T dy for the Toda 
equations (jKHjl . We recall that the conserved current T dx for Eq. (jKHjl 
corresponds to the Noether symmetry ipl = D(Et(T dx)) E sym/^xoda- 
Another (nonlocal) conserved current 



Ua 



: exp{—Uzz) dx A dy + {^u^^ — Ux^) dx A dz, 



for Eq. ()14fjj) (see Eq. (fTTj) on page fT^ could be also treated as an 
analog of the integral (|3Hj) for Eq. (|^3j) due to the following reason. 



SYMMETRIES AND CONSERVATION LAWS IN INTEGRABLE MODELS 97 

Consider the Hainiltonian representation 

Uy = D-^ o Df o E„ (i/heav dx A dz) , (157) 

of Eq. dUni), here 

is the continuous hniit of the Hamiltonian, see Eq. ()117|) . for the Toda 
equations. By Lemma EH the density -ffheav is conserved on the related 
Hamiltonian equation p57|) . therefore the current t] is also conserved. 
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Chapter 4. Backlund transformations and zero-curvature 

representations 

In this chapter, we investigate the relationship between Backlund 
transformations and zero-curvature representations for the hyperbolic 
Liouville equation, the wave equation, and the Liouville scal'''-equation. 
They are 

Uxy = exp{2u), Vxy = 0, Tr,y = exp(-2T), (158) 

respectively. In Sec. [T] on page EHl we we discussed a natural geomet- 
ric scheme that provides Eqs. ()158|) and gives their interpretation. In 
what follows, we construct one-parametric families of Backlund trans- 
formations and zero-curvature representations and consider exmaples 
of integrating Backlund transformations and zero-curvature in nonlo- 
cal variables. 

10. Backlund transformations and their deformations 

In this section, we study the properties of one-parametric deforma- 
tions of Backlund transformations for Eq. ()158|1 . Consider the covering 
Tf-. St ^ Suou defined by the extended total derivatives 

D, = D, + u,—, Dy = Dy + Uy—, [D , , D y] = . (159) 

We assume that the partial derivatives of the nonlocal variable u with 
respect to x and y are 

Ux = Ux + exp(— t) ■ exp('u -|- u), (160a) 

Uy = —Uy + 2exp(t) ■ sinh('u — u), (160b) 

respectively. Recalling Remark 121 on page ^1 let the diffeomorphism 
/i be the swapping u ^^ u oi the fiber variable u and the nonlocal 
variable u combined with the mapping x i— *> —x and y ^^ —y. Then 
diagram (j^ supplies Backlund autotransformation B{£t, Tt, Tt o /^, 
^Liou) for Eq. p87p . and for the Liouville equation. The equations St 
of this Backlund autotransformation (J^) are, obviously. 



{u — u)x = exp{—t) ■ exp{u + u), (161a) 

{u + u)y = 2exp{t) ■ smh.{u — u). (161b) 

Put 

for any A; G N. Consider the scaling symmetry 
OX oy 
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of the Liouville equation and extend X° onto entire S^^^: 

^ k>l '^ k>l 1^ 

Proposition 58. The symmetry X cannot he extended to a symmetry 
of the covering equation St- 

Proof. Assume the converse. By 9^ denote the evolutionary vector 
field 

(J 

on St^ where (p G C'^(£t). By definition, set 

Suppose there is a smooth function a G C°°{St) such that the linearized 
system 

ipif) = 0, D^^{a) = B^,a(MxO = \'^v + "'-Q^] (^xO (163) 

holds. This means that the field 9^,a is a local symmetry of the covering 
equation St and X is extended onto £t constructively. Nevertheless, 
system p63|) is not compatible since 



D^ o Dy{a) 7^ DyO D^{a). 
Indeed, D^ o Dy{a) — Dy o Dx{a) does not depend on a at all and equals 

-2yu^e2*+"+" + 2a;u^.e2*+"+" - xuy^^-'' + xe^^w^ 

-ye^'^Uy + 2e^xul + yUye^"" - 2e^yUyU^ ^ 0. 

This contradiction concludes the proof. D 

Therefore, the scaling symmetry X is a Tf-shadow only {i.e., a solu- 
tion to the equation ipi'p) = 0) and generates the family of the covering 
equations St over Suou'i these £t are parametrized by t G M. 

Consider the Cartan distribution Ct on the equation St and by Ut 
denote the Cartan's connection form on St- In local coordinates, the 
form Ut is 

d f 

Ut = 2J (^c{ua-) (S> -^ h idu- (mx + exp(M + u-t)) dx + 



\ d 
+ {uy - 2exp(t) sinh(M - u)) dym—, (164) 



100 ARTHEMY V. KISELEV 

where dc is the Cartan differential. In coordinates, we have 

i 

Let Q G D(A^(£^)), then fi is the degree of the derivation Q. By 
[■, ■] we denote the Frolicher-Nijenhuis bracket f |37| l62p: 

in, ei™(/) = Ln(e(/)) - (-1)'^^ ■ Le(n(/)), (les) 

where fl,Q e D(A*(£^)), / G C°°{£), and the degrees /i, i/ are /i = 
degf2, 1/ = degG, respectively. We also assume that 

Ln = [in,d]:A'i£)^A''+''''^''i£) 

is the Lie derivative and 

is the inner product. The Frolicher-Nijenhuis bracket is a natural geo- 
metric structure in differential calculus. The de Rham differential d and 
the Richardson-Nijenhuis bracket [■, ■] (see jS2] for its applications) 
are other examples of natural structures. 

Theorem 59 (jSZj). Let t: £ ^ £ be a covering and At: £ -^ £ 
be a smooth family of difjeomorphisms such that Aq = id and Tt = 
T o At: £ ^ £ is the covering for any t G M. Then the evolution of the 
Cartan connection form Ur^ is 

^ = lX„f/J™, (166) 

where Xt is a Tt-shadow for any t G R. 

Suppose £^ is a finite-dimensional manifold, then there is the isomor- 
phism 

D(A*(^))~A*(^)®D(^). 

Thence, any derivation Q G D(A*(£^)) is decomposable to the finite 
sum such that the summands are Q = uj ^ X, where u G A*(£^) and 
X G D(£^). The Frolicher-Nijenhuis bracket of such elements is 

lu^X,e® Yj™ = u A e (^ [X,Y] + u ALx{e) ^ {¥) + 

+ {-ly dujA{x^e)®Y- {-ly^ ^ a Ly (cj) ® x - 

-{-l'f+^)^deA(Y^uj)(^X, (167) 

where X, F G D(^), oo G A'(^) and 6 G A^{£). If the dimension of £ 
is not necessarily finite, then there is the embedding 

A*{£)^'D{£) cB{A*{£)) 

defined by the rule 

(a7®X)(/) = X(/)^ 

for any function / G C°°{£). 
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Consider the coverings Tt: £t ^ £ defined in Eq. fll59p . Then we 
have 

IT J r\ 

—— = (expf-u + u — t)dx — 2 exp(t) sinhftt — u) dy) t— • (168) 
at ou 

We claim that the scahng symmetry X is the r^-shadow such that 
the evolution of the connection form Ut ()lfi4j) in the covering Tf, see 
Eq. ()159|) . is given by ()168p in virtue of equation p66|) . We need 
Lemmas innHnni to prove this. 

Lemma 60. [X, Ut\ _: dn = [dUt/dt) _: dtt. 

FN - FN - FN 

Lemma 61. |X, Ut\ ^ dx = |X, f/J ^ dy = |X, Ut\ ^ dw = 0. 

Proof. The proof of Lemmas 1^01 and |^ is based on successive applica- 
tion of formula ()167|) . The coefficient of d/dx is 

- ^ <\cUa A ljd/du„ {-x) - dcM A ljdidu{-x) = 0; 

the coefficient of d/dy is calculated analogously. Now consider Eq. (jl(')2|l 
and decompose X in the coordinates {x, y, Uk, u-^. We get X = 
^uja^ Xa, where Ua is a 0-form and Xa is a derivation for any a. We 

a 

have i = and j = 1 in Eq. ()167p . Thence we get 

d 

^{uJa A Lx^idcu) + dua A (X^, _: dcu)) (g) — , 

a 

where the first summand is 
y uja A d(XQ, _i dcu) + 2^ ^a A (Xq, _i d(dcn)) = —u^ dx + Uy dy, 

a a 

and the second summand equals u^ dx — Uy dy, whence their sum is 
also trivial. Finally, we calculate the coefficient 

^(wc. A Lx^idcu) + duja A (X„ _: dcu)) 

a 

of d/du by using the explicit formula for dcU- The first summand is 
equal to —u^ dx — Uydy, and the second summand is 

[ux + e~* exp(-u + m)) da; + [uy — 2e* sinh('u — u)) dy. 

Consequently, we obtain the required expression 

^ d 

(e * exp(M + u) dx — 2e* sinhfw — u) dy) Tr:r. 
^ ou 

This concludes the proof. D 

FN 

However, the calculation of coefficients of d/duk or d/du-^ in [X, Ut\ 
is untrivial for k> 1. First, we have 
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Lemma 62 ([47 ). Let u{x) and f{u) be smooth functions and D^ be 
the total derivative with respect to x. Denote u^ = D^{u{x)), k > 0, 
Uq = u. Then the relation 

n . DlU\u)) = Y,^^rn -j—D:U{u)) (169) 

m=l '^"™ 

holds for any integer n>l. 

The proof of Lemma in21 is based on Lemma 1^21 and Corollary IMl 
below. 

Lemma 63. Let the assumptions of Lemma\6^ hold. Take an integer 
n > and a positive integer I < n — 1. Then the relation 

is valid. 

Corollary 64. Let the assumptions of Lemma IMl hold. Then we also 
have 

(n + !)«„+! ^D:+\fiu)) = {n + !)«„+, ^D^ifiu)) = 

OUn+l OUn 

= (n + l)n„+i-f(n). (171) 

Proof of Lemma W^ i [47' ) . We prove ()169|) by induction on n. For 
n = 1, relation (I169D holds. For ra > 1, we have 



(n + 1) D:-^\f{u)) = D^{nDlU{u)) + DlU{u))) 
by the inductive assumption, 

= D, [Y, ^n^ ^ Dl{f[u)) + Dl{f{u)) 



by the Leibnitz rule 
du 



J2mu^^^,^^D:ifiu)) + 

m=l 



+ J2 ^^m D,—- D^ifiu)) + D, D^ifiu)) = 

rra=l 

by Eq. ()170|) applied to the second sum, 

" /) ^ Pi 

m=l m=l 

- E ^"™ 7[. mfiu)) + D^D^ifiu)) = 

___ 1 C''Um— 1 
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by the definition of D^ and tlie subscript shift in the latter sum, 

= E '^^^ ^^ ^r'(/M) + Y.^m+ i)u^^, —- Di{f{u)) - 

m=l "* m=0 ™ 

m=U 

since almost all summands in the latter two sums coincide, 
= E ^^™ ^ DT\f{u)) + (n + l)«„+i ^ Dl{f{u)) = 

m=l 

bydnH), 

^^ dUm OUn+1 

n+1 rj 

= E^«-7^^^'(/W)- 

m=l 

The proof is complete. D 

Here we offer another proof of Lemma 1^21 This proof is based on 
the weights technique. The idea of this proof was communicated by 
V. V. Trushkov. 

Proof of Lem,m,a W^ ( j58j ) . By definition, put the weight wt{uk) = k, 
wt{uk ■ ui) = k + I, and wt{uk^ + u^J = ki if ki = k2. We have 

n 
D::{f{u)) = Y,Pn,m-f'-\u), (172) 

m=l 

where Pn,m = Yl const (ra, m) ■ Uj-^ ■ . . . ■ Uj . We claim that Pn,m is a 

f 
differential polynomial such that 

ji + • • • + ji{n,m) =n VJ, V n, m. (173) 

We prove this fact by induction on n. Indeed, if wt(P„^m) = n, then 
wt{Dx{Pn,m))) = n + 1 due to the Leibnitz rule. Besides, 

n 
m=l 

and therefore the weight wt(D"^^(/(M))) is well defined and equals 
n + 1. 

Now, consider the weight counting operator 

^=E"^-«-9^ (174) 

m>l 
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that acts on the right-hand side in Eq. ()172|) : 

^m-u„,-—o^Y^ const (n, k)-Uj,-...- u,-,^^^^ ■ f^\u) = 

m>l k=l J 

n 

= Y^Y^coTisi{n, k)-n-Uj,-...- m^- ^^ ^, ■ f^^'^u) =n-D2{f{u)), 
k=i J 

since condition ()173p holds for all multiindexes j. We finally conclude 
that the functions D^{f{u)) are eigenfunctions for operator p74|) and 
the integers n G N are the eigenvalues. Therefore, relation ()169|) is a 
solution to the problem X ■ (f = W{<f). The proof is complete. D 

FN - FN 

Lemma 65. |X, t/J _i duk = {X, Ut] _i duj: = for any k > 1. 
Proof. Let k eN. Taking into account Eq. ()167j) . consider the 1-form 

lX,Utj™ ^duk= ((A;-l)5,M,-5^/Mi— D^i(exp(2M))) -dy. 

We see that all coefficients of dx, du, dui, duj vanish for all / > 1. Fi- 
nally, note that DyU^ = D^~^{exp(2u)). By Lemma 1^21 the coefficient 

FN 

of dy is trivial. Arguing as above, we conclude that [X, Ut] _i duj^ = 0. 
This completes the proof. D 

Theorem 66. The Tt-shadow in Eq. p62|) satisfies the relation 

FN _ _ 9 

[X, UA = (exp(-u + u — t) dx — 2 exp(t) sinh('u — u) dy) ® 7—, 

ou 

i.e., the abelian group of difjeomorfisms At = exp(tX) provides smooth 
one-parametric family ()160|) of one-dimensional coverings over Liou- 
ville's equation Suou- These coverings correspond to Bdcklund autotrans- 
formations (jlfjlll for the Liouville equation, see diagram (J2S1)- Evolu- 
tion of the connection form is given by ()l(i8|l . 

FN Q 

Proof Decompose the bracket |X, Ut] with respect to the basis (^, 
^, g^, ^, ^). By Lemmas inU and inni all coefficients of the deriva- 
tions ^, ^, g^, -£^ equal for /c > 0. Thence, Lemma lUUl supplies 
the required expression, which is given in Eq. ()168j] . D 

Remark 25. Consider Backlund transformation (|14j) 

(f — u)x = e~* exp(M + v), (f + u)y = — e* exp(M — v), (175) 

t G M, between the Liouville equation (f76|) and the wave equation 
Vxy = 0, as well as Backlund transformation (|14j) 

(T - u)^ = 2e-* cosh(T + u), (T + u)y = -e* exp(M - T), (176) 
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t e M, between £^Liou and the Liouville scal^-equation T^y = exp(— 2T). 
Then these transformations possess the property similar to one de- 
scribed in Theorem IBHl Namely, scaling symmetry ()162|) of Suou is 
the required r^-shadow in both cases. The proof of these theorems is 
quite analogous to the proof of Theorem IHSl above and appeals to useful 
identity p(J9j) in total derivatives. 

11. Integrating Backlund transformations in 

NONLOCAL variables 

In this section, we study nonlocal aspects of integrating Backlund 
transformations between PDE. The aim of this section is to illustrate a 
natural scheme that provides nonlocal variables associated with a cer- 
tain PDE and to obtain nonlocal structures for the hyperbolic Liouville 
equation. We reconstruct the coverings Tj from Eq. p(il|l and ()175j) - 
p7(ij) and demonstrate the nonlocal variables to be potentials for the 
fiber variables u, v, and T. In what follows, we use the notation Su as 
a synonim of Suou, meaning that the Liouville equation u^y = exp(2M) 
is imposed on the fiber variable u. 

Now we construct the one-dimensional non-abelian coverings such 
that we can integrate Backlund transformations (HHH), p75|) - ()176p in 
the corresponding nonlocal variables. We fix an arbitrary t G M and 
consider the extended total derivatives 



d ~f =f ,„o ^ „ . (9 



^x ^x ^ pr^ ' ^y ^y ^ l"oo ^^ ^"y^oo; o^ ? 

C/i — ^oo L/i — ^oo 

Dj" = Dl- - e^" — , bl- = Dl- + (Sf + 2«.S; + e-'*)^, 

Df" = Dl^ + e^^-^, bl^ = m^ + {2vyZ^, + e^*^ — 



g^v^ y y ' \--y~t • - ^ Q'z:v 



I>-^Z>- + e-^, Z,?^4- + ((Hn-2T- + e-^ 



(m) 

We see that the extended derivatives commute in all cases: 

[b.,,by] = o. 

Therefore the coverings 

T ■ F ^ coo -T ■ P -^ P°° t' ■ F' ^ P°° 

't-'-'t ' '-'u 1 'oo • '-'oo '-'u 1 't ■ '-'t '-'u 1 

Tl.£",^£^, tJ:£J^£^ (178) 

are well defined. The explicit form of the covering equations £^t, £ooi 
£1, £^, and £^ is discussed in Remark 071 
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Remark 26. Coverings (jl78|) with nonlocal variables ()177|) cannot be 
reduced to local conservation laws for the underlying equations S^, Sy, 
and Sf] these coverings are called non-abelian. We also claim that the 
t-parameterized coverings {e.g. r^ at the points ti and ^2) are equivalent 
(xf^ ~ Tfj), i.e., there is a functional dependence between the nonlocal 
variables (S^^ and Et^ in our case). We have 

Sij + y ■ exp(2ti) = Et^ + y- exp(2t2) = Si=_.oo Vti, t2 G M; 

similar relations hold for other coverings ()178j) . 

Remark 27. The covering equations can be obtained explicitly since 
each nonlocal variable in p77j) is a potential for at least one of the 
dependent variables u, v, and T. For example, 

Now we derive the covering equations imposed on the fiber variable S^ 
and its limit Sqo at the point t = —00: 

£, = <( fi = H? + =^ - exp(2t) } , (179) 



dy 



dx 



where t G M U {—00}. The equations Sj., £^, and £^ are obtained 
analogously. 

11.1. Integrating in nonlocal variables. We emphasize that trans- 
formations (HnH), (!175|) . and ()176|) cannot be integrated in local vari- 
ables. Therefore we consider one-dimensional non-abelian coverings 
fll78j) and extend the sets of (local) variables with new nonlocal ones, 
see p77|l . and then we integrate the transformations successfully. The 
results can be summarized as follows. 

Theorem 67 ([SHI)- Consider Bdcklund (auto) transformations p61|) . 
fll75|) - ()176|) for the equations S^, £v, ond £-^. These transformations 
are integrated in nonlocal variables explicitly. 

(1) Bdcklund autotrans formation (11611) for £u- 

u = u + t — log St u = t + u — \ogEt[u]{—x, —y), 

i.e., to inverse the transformation and obtain u[u], the inversion 
X 1-^ —X and y t— > — y is required. 

(2) Bdcklund transformation ()175|1 between Su and the wave equa- 
tion Vxy = 0: 

V = u + t — logSoo and, conversely, u = v + t — logS^. 

(3) Bdcklund transformation ()176|) between £u and the scaX^ -Liou- 
ville equation T^y = exp(— 2T): 



T = —u + t + logS^ and, conversely, u = — T — t — logS 
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Proof. We consider the case u[u]{x,y) within Backlund autotransfor- 
mation ()161|) . By definition, put U = exp{u) and T = exp(— -u). From 
Eq. ()161|) we obtain the Bernoulh equation 

i^x = Ux -U + exp(M — t) U^, 

whence U^^ = T = exp(— m — t) • S, where the nonlocal variable S is 
such that -Da;(S) = — exp(2M), and also get the Riccati equation 

7:y = Uy-T + exp(u + t)T^ - exp(t - u). (180) 

Substituting exp(— n — t) ■ H for T in pSOj) . we get Dy{E) = S^ + 
2uy'E — exp(2t). Now we refer p77|) and compare the result with the 
derivatives Dxi^t) and Dy(Et). 

The proof of other five cases is quite analogous. Assume /(x, y) G {u, 
it,, V, T} is a known solution to the PDE Sj. Then we obtain either 
two Bernoulli equations for Eq. p75|) or one Bernoulli equation and one 
Riccati equation for Eq. p61|) and Eq. p76|) after a proper change of 
variables. Solving these ordinary differential equations for the solution 
g{x, y) G {m, u, t>, T} to the PDE £g that is related with £f by one 
of Backlund (auto)transformations ()175|l - ()17(')j) . we finally obtain the 
rules to differentiate the nonlocal variable, which appears in one of the 
coverings p78|) . D 

Consider diagram ()24j) that arises in the definition of a Backlund 
transformation and apply Eq. pK|l to Theorem IHZI Consider the cov- 
erings Ti and T2 in Eq. p77|) . We see that in all cases one of these two 
mappings is a first order differential operator that depends on the non- 
local variable only, while another projector is a zero order morphism. 
Thence we have 

Theorem 68. Consider equations p79|) . Then the following relations 
hold: 



u = \\og{-Et)^, u = t + log[S,-y(-H,)J, 

M = i log (-Hoo)^, v=t + \og [H^^ -^(-Hoo) J 

ilog(-S;) . T=t + log[^^^], 



u 



|log(SX ^ = i + log[(Hn"^V^]^ 



T = -l^og^y. u = -t + log[(-^)-y(S7)J. 

In other words, the nonlocal variables that satisfy Eq. p79j) are po- 
tentials for both solutions /, g to the equations Sj, Sg within the list 
Su, Sv, and £r- The property of all the coverings in ()178|1 to be nonlin- 
ear differential operators of order not greater than 1 is a special feature 
of these equations. 
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11.2. On nonlocal symmetries. Let E be an equation and r: E ^^ E 
be a covering. By definition, a T-shadow ipis a solution of the linearized 
equation i£{(p) = 0. The shadow fields 3^ G Dc(^) are not true 
nonlocal symmetries since they do not describe the evolution of the 
nonlocal variable and, generally, not all of them can be extended up to 
nonlocal symmetries. 

We have proved that Backlund transformations ()1(J1|1 . p75j) . and p7(ij) 
themselves do contain certain information about nonlocal variables 
such that these transformations can be integrated successfully. The 
corresponding non-abelian coverings in p78|) provide nonlocal conser- 
vation laws for the underlying differential equations. Still, the struc- 
tures on the covering equations are "too close" to the initial ones, so 
that the point symmetries of the initial equations and the classical 
symmetries of the covering equations are in one-to-one correspondence 
such that we get no nonlocal symmetries except the liftings of local 
t r ansf ormat ions . 

Now we introduce new nonlocal variables such that the symmetries 

we are in search of depend on them. Let Et= Sj + Uy be the new 
nonlocal variable such that 

D^y^ilJt) = (St) + Uyy -u\- exp(2t). 



^— * 



Consider the limit of Z"* as t — ^ — oo. We see that at the point t = — oo 
there appeares the automodel variable 

_ , exp(2n) 

^oo — Ux -\ — 



such that i)^"(Soo) = 0. We claim that 

Uoo= hm Et 

t— ►— oo 

and Soo differ by the discrete symmetry x ^^ y. Indeed, consider the 
derivatives D^ and Dy of Sqo and Eoo- We get 

Also, we have 



DxiSoo) = 0, DyiSt) = (Zoo)' + Uy 



ul. 



"yy "'y 

Therefore we use the nonlocal variable Sqo only and treat all relations 
up to the symmetry transformation x ^->- y ioT the Liouville equation. 

By definition, put T,t = {x '^ y) ■ (Et)'- we have 

Dxi^t) = Sj + Uxx -ul- exp(2t) and DyiT.t) = 0. 
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Nonlocal variables enable us to find shadows of nonlocal symmetries for 
Eq. (jTBj) and to extend these shadows up to true nonlocal symmetries 
of the Liouville equation. 

Proposition 69. (1) Let f{t,x,T,t) be a smooth function. Then 
the generating function 

1 df Idf 

if = -(Si^ + u,^ -ul- exp(2t)) . ^ + -^ + n, . / (181) 

is a Tf-shadow of a nonlocal symmetry of the Liouville equation. 
(2) Let f{x, Soo) be a smooth function. Then the second order r- 
shadow (f{x,IJoo,u,Ux,Uxx) for the Liouville equation is 

1 df 1 df 

ip = -(S^ + uxx -0-^+2'£ + ^^f = °(/(^' ^oo)). (182) 

Nonlocal shadows i\lSl\\ and fll82j) belong to the class {this class was 
considered in the paper jnHI) of solutions (jF)H|) to the equation ipif) = 
0. This class is now provided by the operator 

that contains the extended total derivative D^. 

Recontruction of nonlocal symmetries. In order to extend the r^-shadows 
Q^ up to true nonlocal symmetries 

d 



where a G C°°(£^) and t G R U {— oo}, we solve the equations 

Dxia) = 9^,,p,(S0), Dyia) = 9^,,(£),(S0) 

for the function a. 

Proposition 70. (1) Let f{t) be a smooth function and the func- 
tions if and a{t,T,t,Ux,Uxx) be defined by the relations 

^ = Ux- fit), a = (S^ + uxx -ul- exp(2t)) • fit). (183) 

Then the field 3^ + a ■ d/dHt is a true nonlocal symmetry of 

Eq. dzni). 

(2) Let fix) be a smooth function and the functions ip and a(Soo, 
Ux, Uxx) be given in 

V' = «./(a^) + ^£,a=(SL + «..-^D/(^) + Soo£ + ^0. (184) 
Then the field 3^ + a ■ d/dTjao is a true nonlocal symmetry of 

Eq. dzni). 
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The proof of Propositions inHl and [701 is very extensive and cannot 
be completed without apphcation of the Jet ([72j) software for analytic 
transformations. The Jet environment allows to set the determining 
equations in a dialogue mode, select the simplest differential conse- 
quences to these equations, and then specify the expressions for the 
unknown symmetries. 

Nonlocal symmetry ()183p is defined up to elements of CD(£^) 3 g-Dx, 
g G C°°(£^). Thence, the nonlocal symmetry class (jl83p is [9i^,a] = 
[— /(t) ■ d/dx], where f{t) ■ d/dx is the translation. Symmetry ()184|1 
is the lifting of a classical point symmetry (^q, see Proposition 1111 on 
pagelHHl As usual, a class of nonlocal symmetries 3ip,a can be obtained 
from Eq. p83j) and ()184j) by using the discrete transformation x ^^ y. 

11.3. On permutability of Backlund transformations. Now we 

illustrate the permutability property of Backlund (auto)transformations 

(Hnn), (HzsD, and diini). 

Proposition 71. (1) Let u\ j = i, ii, be solutions to Eq. (f7B|) such 
that Bu{u,u^;tj) = 0, tj G M. Then there is a unique solution 
u"'{x,y) to the system 



Buiu',u"';t2) = 0, 
B^{u",u"';t^) = 0. 

Namely, the solution u'" satisfies the relation 



;i85) 



, ,„. , . k2exp{u') -kiexpju") 
exp M = exp -u ■ z — ^ ; ; — -, 186 

where kj = exp{tj). 
(2) Let j = i, ii and tj G M. Suppose both v^ are solutions to the 
wave equation Vxy = such that Buv{u,v^]tj) = 0. Also, sup- 
pose that u^ are solutions to the Liouville equation such that 
Buv{u^ ,v]tj) = 0. Then there are unique solutions u'" and v'" 
to the systems 

Buviu"',v';t2) =0 J ! Bu.{u',v"'-t2) = 

Bu,{u'",v";ti) =0 ana <^ B.,,{u",v"'-t,) = 0, 

respectively. Denote kj = exp{tj), then the following relations 
hold: 



exp{u") = exp(M) ■ 



.'i'\ 



exp[v ) = expl^fj 



^2 exp(v') — fci exp(t>") 
/c2 exp(f ") — fci exp(t>') ' 
fci exp(M") — k2 exp(u') 



k2 exp(u") — ki exp(u') 

(3) Let j = i, ii and tj G M. Assume that T^ are solutions to 
the scal'^ -equation Sr such that Buriu^T^^tj) = 0, and sup- 
pose that u^ are solutions to the Liouville equation such that 
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Bur{u-^, T;tj) = 0. Then there are unique solutions u'" and T'" 
to the systems 



Bur{u"',r-M) = 
But{u"',T'-M) = 



and 



Bur{u',r"-M) 
Bur{u",T"-M) 




0, 



respectively. We also have 



exp(-u"') = exp('u) 
exp(T"') = exp(T) 
where kj = exp(tj). 



A;2exp(r)-A;iexp(T") 
k2 exp(T") — ki exp(T') ' 
ki exp(-u") — k2 exp('u') 
k2 exp(-u") — ki exp('u') ' 



Proof. We consider the case of Backlund autotransformation pfilj) only; 
cases El and El are treated analogously. Consider the subsystem in 
Eq. ()185j) that consists of relations pCilj) which contain the derivatives 
with respect to x only. Then the solution u"' defined in Eq. (I18(i|l fol- 
lows from the linear dependence of the left-hand side within Eq. ()185|) 
and is unique. Consider another subsystem composed by the relations 
in Eq. p61|) that contain the derivatives with respect to y. Then there 
are two solutions to this subsystem; we denote them by u'" and u'". 
The solution u'" is defined in Eq. (I186|l and 
relation 



u'" is obtained from the 



exp(M 



exp(— uj 



ki exp(M') — k2 exp{u") 
k2exp{—u') — kiexp{—u") 



The solution u"' is irrelevant. Therefore, the function u"' is a unique 
solution to the whole system ()185|) . D 

Remark 28. Proposition ITTl means that the diagrams 



u 



> u 



u 



> V 



*3 // 
> U 



te 



t2 



t2 



U 



-^ u 



-* u 



-> V 



tl 



i3 



and the diagram 



u 



t2 



-y// 



Tf ^3 // 
> U 



t2 



-y u 



are commutative for any ti, t2, t^ G R. 



t2 



-, T'" 
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12. Zero-curvature representations 

In this section, we illustrate the relationship between the parametric 
families of zero-curvature representations and Backlund transforma- 
tions for Eq. ()158|) . In Sec. [l] we assigned the Liouville equation to the 
Lie algebra q = sl2{C), and now we benefit from the use of different 
representations of this algebra. 

There is a natural equivalence (flBJ) between g- valued zero-curvature 
representations of a PDE S and the special type coverings over the 
equation £. Further on, we study the case g ~ sl2{C). It is essential 
that there is a representation of sl2{'C) in vector fields. We use it in 
order to construct the required coverings over the hyperbolic Liouville 
equation 

^Liou = {uxy = exp(2M)}. (187) 

Here '2' is the Cartan 1 x 1-matrix of the Lie algebra Ai while x and 
y are the coordinates in the standard two-dimensional extension [z, 
z) ^^ C'^ 3 (x, y). By (e, h, f) we denote the canonical basis such that 

[h,e]=2e, [h,f] = -2f, [e, f] = h. m) 

Consider the representation p: sl2{C) -^ D(C2[H]) of the Lie algebra g 
in the space of polynomial- valued derivations: 

P(e) = l-^, pih) = -2E.-^, p(/) = -S2.A (188a) 

such that the Lie bracket is the commutator of vector fields: [A, B] = 
AoB — BoA. This representation of the Lie algebra sl2{C) was used in 
the paper [52j, where a class of A^-ary analogs for the Lie algebraswas 
constructed. Also, consider the matrix representation 

^(e)=(o J), m=(l _°i), ^(/) = (? o) (188b) 

such that the Lie bracket is the matrix commutator: [A, B] = A ■ B ~ 
B-A. 

Given an equation £, consider the flat connection form ()47|) in the 
bundle C°°(£'°°) ® G — > C^, where G is the Lie group of g. Suppose 
A^ B E C°°{£°°) g. The zero-curvature condition pHj) . which is 
equivalent to the relation 

[D, + A, Dy + B]=0, 

is satisfied on the differential equation S. Finally, we obtain the matrix 
equation 

DyA-D^B -[A,B]=0. (189) 

Now, decompose the matrices A and B with respect to the basis in the 
representation g: q ^ {M e Mat(2,2) | trM = 0}: 

A = ae®Q{e)+ah®Q{h)+af®Q{f), B = be®g{e)+bh®g{h)+bf®g{f), 



(190) 
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where a^, bi, G C°°{S°°), and construct the one-dimensional covering r 
over S such that S is the nonlocal variable, the extended total deriva- 
tives Dx and Dy are 

D^ = D^ + tte ® p(e) + ah(S) p{h) + aj (g) p{f), 

Dy = Dy + be0 pie) + hh® p{h) + hf® p{f), 

respectively, and the rules to differentiate the variable S are 

D.^(E) = dx_i (ae(8)p(e) +ah® p{h) + aj ®p(/)), 

Dy{E) = dy^ {be ® p(e) +bh® p{h) + 6/ ® p(/)). 

We see that the Maurer-Cartan condition ()48|). which is satisfied on S, 
is equivalent to the compatibility condition [D^, Dy] = that holds in 
virtue of the equation £^. 

Example 22. First, we obtain Backlund transformation between the 
Liouville and the wave equations. Consider Eq. (j3^ on page OH and 
choose the gauge 

tte = al = exp(KM), bf = b}f = exp((2 — k)u), 

for an arbitrary constant k. Then the covering equation S is 

X = {k, — 2)ux + exp{K,u — V 



X f • (191) 

Vy = KUy — exp((2 — k)u + v) \ 

Here the variable v = log S is a transformation of the nonlocal variable 
S, see Eq. fjl9Up . The compatibility condition for Eq. ()19ip is 

Vxy = (/« - l)exp(2u). 
Finally, suppose k = 1, then Eq. p9H) is the Backlund transformation 

(v + u)x = expiu — v), 

( 1 ^ cm) 

(y — u)y = — expi^M + V) 

between the Liouville equation Eq. p87p and the wave equation 

Vxy = 0, (192) 

while the coordinate S is exactly the one that allows integrating equa- 
tion ()191t ) in nonlocal variables, see Sec. ^2 

Remark 29. Transformation p91l ) is the particular case t = 0, k = 
exp(t) = 1 within the family of Backlund transformations ()175p be- 
tween Eq. ()187|) and Eq. p92|) . We note that the mapping k h-> 
—k is the replacing of representation ()188ap with the representation 
p: 5[2(C) ^ D(C2[S]) such that 

p(e) = -1, pih) = -2S, p(/) = E'. 



114 ARTHEMY V. KISELEV 

Zero- curvature representations constructed by using Bdcklund transfor- 
mations. Consider Backlund autotransformation ()161|) for the Liouville 
equation ()187p and Backlund transformation ()176p between the Liou- 
ville equation and the Liouville scal'''-equation 

£r = {T^y = exp(-2T)} 

(we recall that the latter equation is the hyperbolic representation of 
the Gauss equation for the conformal metric of constant curvature +1, 
see Example IIUI on page EH])- Next, consider the coverings over the 
Liouville equation that provide these transformations. Then, these 
coverings exceed ansatz ()5H). 

Namely, the flat connection form for Backlund autotransformation 

(Hnn) is 

exp[u — t) 2'^xJ y—exp{t — u) ~2'^y J 

The s [2- valued fiat connection form 6 for Backlund transformation ()176|1 
is 

9 = ( --f^ ,, ^^P^r' " ^^1 da: + C^^y - "M' + ^^l d^. 
\^-exp(M-t) ^Ux J V -^Uy J 

By construction, these s[2-valued forms are zero-curvature representa- 
tions for the Liouville equation. 

Bdcklund transformations constructed by using zero -curvature repre- 
sentations. The problem of constructing multi-parametric families of 
Backlund transformations by using known zero-curvature representa- 
tions for Eq. p87|) was discussed in the paper [22]; later, this problem 
was studied thoroughly by V. Golovko in j5^ . 

Now the exposition follows [BQ', and |H2]. Let us find three s[2-valued 
classes of zero-curvature representations 9 G A^(£^£?,^) (8>s[2(C) for the 
hyperbolic Liouvile equation p87j) . Assume that A = A{ux), B = 
B{u), and suppose [A,B] ^ 0. Then Eq. (jlHI) on page 11121 is reduced 
to 

1 dA , , dB r 1 ■ / s -, 
u^^-T exp(-2n)- [M^M,exp(-2n)5] = 0. 

Applying d"^ /dudux to this identity, we get the equation [M, A^] = 0, 
where 

_ d{A/u,) d{B/eM2u)) 

OUx OU 

There are three possible cases: 

(1) M = 0, 

(2) AT = 0, and 

(3) M = r{ux) -C, N = s{u) ■ C, where C ^ is a constant shiR)- 
valued matrix nd r, s G C°^(]R) are smooth functions. 
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Hence we get three gauge non-equivalent clases of zero-curvature rep- 
resentations: 

Case 1 (M = 0). There are no nontrivial solutions to Eq. ()49l ). 

Case 2 (A^ = 0). There are two classes of zero-curvature representa- 
tions. They are 

/ 2au^ + 2/3 2a 

\ulil-2a)- 4(3u^ + 27 -2au^ - 2/3^ 

« = (exp(2«) o) • d''^'^) 

A_( aul + 2P 26exp{—2aUx 
~ \^27 exp(2aM^) -2aul - 2(3 

^^/aexp(2«) \ 

y — aexp(2M)y ^ ' 

Case 3 (M ^ 0, A^ ^ 0). Another solution to Eq. (gni) is 

where a, /?, 7, and 6 are arbitrary constants. 

By jIH!, any 5 [2- valued zero-curvature representation for S provides 
some special type covering over the equation S. Consider the represen- 
tation p': s [2(C) ^D(C[H]): 

p'(e) = exp(-.)|-, p'if) = -exp(.)|-, p'ih) = -2^, 

of the Lie algebra sl2{'C) in the space of differential operators on C, i. e., 
I' G C is the nonlocal variable and p' : sl2{C) -^ D(C[[u]]). Then we con- 
struct one-dimensional coverings over the Liouvile equation ^Liou that 
provide Backlund transformations between £^Liou and some differential 
equations that depend on the initial zero-curvature representation. 

Proposition 72 (|32I). Representations (ITTIHajl . (IT93b|) . and (fTMj) cor- 
respond to Backlund transformations between Euou o-nd the equations 

+ 



+ 2aexp(— f') — 4/5 — 27exp(t>'), 

v''\-2^ex^{^+v" 

^.y -exp(-2.;'")«^ + 4a). 

If a = in representations fll94j) . then we get Backlund transformation 
between Suou ^^t? the wave equation Eq. 



< = 


i(2«-l)exp(t;')(^ 






2{(3exp{v')-a)(^-^- 


-<) + 


K = 


av''^ I 

\ 4/5 + 25exp 
2vf \ 
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Finally, we analyse the removability of the parameters a, /3, 7, and S 
in zero-curvature representations ()193|) - ()194p with respect to the gauge 
transformations 

A ^ SAS-^ - {D^S)S-\ B ^ SBS-^ - {DyS)S-\ 

The result is described in 

Remark 30 ([^)- The parameter /3 in Eq. ()193ap is removable by the 
gauge transformation 

that depends on an arbitrary constant a G C; under this transformation 
7 I— > 7 + p'^/a. All other parameters a, P, 7, and 6 in zero-curvature 
representations ()193p - ()194|) are non-removable. 
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Appendix A. Geometric methods of solving 

BOUNDARY-VALUE PROBLEMS 

In this appendix, we consider several methods for solving boundary- 
value problems (mainly, the Dirichlet problem) for nonlinear equations 
of mathematical physics. Then we compare and analyse the results 
of a computer experiment in applying the described algorithms, see 
also jlH]. Again, the methods which we study are based on the geom- 
etry of the jet spaces (PH ESI EI]) and treating differential equations 
£ (and their prolongations £°° as well) as submanifolds £ C J'^(vr) of 
the jet space of order k for a certain fibre bundle vr (respectively, we 
have £°° C J°°(7r)). By using this approach, we pay attention to the 
analysis of the following objects: they are 

• differential equations and sections s such that the jet jfc(s) C £ 
defines a solution of £, 

• boundary-value problems V = {£,V, f = s|^^) in a domain V, 

• and their deformations V. 

We emphasize that the geometrically motivated technology of solv- 
ing the boundary-value problems is discussed within this appendix. 
We do not stop on particular properties of some solutions in practi- 
cal situations. Basic definitions and concepts were formulated in the 
Introduction, see also jTUl ESI 113 EI]- We use the elhptic Liouville 
equation 

^Liou = {uzz = exp(2M)} 

as a basic example. 

The appendix is organized as follows. In Sec. lA.lf the method of mo- 
notonous iterations for the solutions m* of the boundary-value problem 
V is described. Here t G N. This method is based of the theory 
of differential inequalities ([Znj)- In Sec. IA.21 we construct a method 
for solving the boundary-value problems that involves the evolution 
representation of the equation £ under study. We shall see that this 
interpretation is admissible in a sufficiently general situation (|in[l81j). 
Then, in Sec. IA.3l we consider various methods based on the deforma- 
tions / of the boundary conditions / and simultaneous invariance of 
the equation £: 

V = {£,Vj). 
In Sec.lA.4lwe describe the relaxation method based on the substitution 



of the boundary-value problem such that solutions u of the problem V 
are stable stationary solutions for the new problem V' with respect to 
the evolution equation £' . The deformations 
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which are considered in Sec. IA.51 are in some sence antipodal to ones 
described in Sec. IA.3I Now, the boundary condition / is invariant and 
the equation S{t) moves. The "planting" of a nonlinearity could be an 
example. In the final section of this appendix, we discuss the results of 
a computer experiment in practical application of all these methods. 

Remark 31. In what follows, we assume that there is a unique classical 
solution to the boundary- value problem at hand (see [Ml)- Also, we 
recall that the transformation laws of solutions with respect to trans- 
formations of the independent variables are known for some equations, 
e.g., the conformally invariant Toda equations, see Proposition 1111 on 
page I38|). Therefore, the domain D, where the problem V is solved, 
can be chosen relatively regular (we set T> r^ Bq for Eq. (j33j) ). 

Remark 32. First, we point out a practically useful way to construct a 
set of solutions (that can be large enough) of a boundary- value problem 

V = {S, V, /). Assume that the equation S is invariant with respect 
to a vector field X. Denote its generating section by ipx and suppose 
that At is the fiow of the field X such that 

A{dV) = dV and A(jfc(/)) = Jkif), 

that is, the boundary condition is unvariant with respect to the sym- 
metry X of the equation S at hand. Then problem V can be reduced 
to finding solutions of the boundary- value problem 

v' = {sn{^x = o},v,f) 

that are invariant with respect to X at any point in V. If the problem 

V has a solution, then, in general, this solution may be not unique. 

Example 23. Consider the Dirichlet boundary- value problem 

{u,, + Uyy = eM^u), u|,=i = /, /gC°'^} (195) 

in the unit disc Bq and assume that the boundary condition is homo- 
geneous and trivial, / = 0. The solutions for the Liouville equation 
^Liou that are invariant with respect to its point symmetries are 

(196a) 
(196b) 
(196c) 







u = 


2 log 


vl + ^l 
sinh^f ' 






u = 


2 log 


vl + vl 

y2 ' 






u = 


^log 


vl + vl 

■ 2 ' 

sm f 


where v 


is 


a harmonic function. 










Av = 


= 0. 



Consider the radial symmetry 



^ d d 

^ = ~yTr^^^■ 
ox oy 
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Then the solutions ■u(r) = u{^^/x^^-y^) of problem p95|) are the fol- 
lowing: 

'W2(^) = — logr — log(l — logr), (197b) 

1 a^r"^ ±1 

^3(^) = -log . ^ -, fj = — arcsma, a 7^ 0. (197c) 

2 sm (ap — alogr) a 

A unique classical solution ui has no singularity at the point 0. The 
solution U2 has an integrable singularity at r = 0, and its gradient on 
the border vanishes. The solution u^ has the cardinal set of logarithmic 
singularities along the radius r. These singularities accumulate as r — ;> 
0. 

Remark 33. Suppose that the equation £ at hand is Euler. Then, 
obviously, one can search solutions of the Dirichlet boundary- value 
problem by using the direct minimization of the action functional. We 
emphasize that the projective methods can be used in addition to the 
lattice discretization methods. 

A.l. The monotonous iterations method. The monotonous itera- 
tions method is a useful instrument in the theory of differential inequal- 
ities (pOl)- This method allows to contruct solutions of the Dirichlet 
boundary-value problems 

{Am = h{u, x), u\g^ = /, dVe Cl+^ f,he C°'^} (198) 

and then check the local uniqueness of these solutions. 

Definition 13. A function a e C(V) n C^iJ)) (resp, /3) is a lower 
{upper) solution of problem p98|) is the following two conditions hold: 

(1) Aa - h{a, a;) > in "D and 

(2) /(x) > a{x) on dV 

(respectively, <). Suppose further a{x) < [3{x) for each x eV. Then 

we set a ^ p. 

Proposition 73 (|Z0])- Assume that a lower and an upper solutions 
of problem i]198\i are ordered, a ^ (3. Suppose there is a nonnegative 
constant C eM such that 

h{ui,x) — h{u2, x) < C • {ui — U2) 

for any x E V, Ui, and U2 provided that a < U2 < Ui < j3 . Consider 
two sequences U_ = {u!'} and U = {u^} of solutions of the problem 

/\u^ -Cu^ = h{u^-\x)-Cu^-\ u\^ = f, xeV,ken. 

(199) 
Here the inital values are v^ = a and vP = (3, respectively. Then the 
monotonously nondecreasing (nonin creasing) sequence U_ (U ) converges 
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to a solution u (resp., u) for problem fll98p . Moreover, we have u<u 
in V and the hound u<u^.<u holds for any solution u^. G [a, [3] . 

Suppose further that h is monotonous with respect to u\ h{ui,x) — 
h{u2,x) > provided that a ^ U2 ^ ui ^ p. Then the solution for 
problem p98|) is unique on [a,/3]: u = u. 

The notions of a lower and upper solution, which were introduced in 
Sec. lA.H and the method of proving the local uniqueness for solutions 
of problem ()198|) will be used in Sec. I A. 41 The relaxation method will 
be described there. In that case, the evolution of the superscript k is 
continuous unlike in problems 



A. 2. Evolutionary representation of differential equations. The 

following remarkable result was in fact obtained in the papers concern- 
ing the formal theory of differential equations (see ^21 IM] and refer- 
ences therein): any equation that satisfies the assumtions of the '2-line 
theorem' for the Vinogradov's C-spectral sequence (IHIlinSI) admits a 
representation in the form of an evolution equation. An illustration is 
given in Example IHl on page 1121 In practice, this means that any equa- 
tion that does not have gauge symmetries (unlike the Maxwell, the 
Yang-Mills, and the Einstein equations and similar systems) admits 
a set of coordinate transformations and differential substitutions that 
map it to an evolution equation (possibly, the new equation is imposed 
on a larger set of dependent variables) . The class of equations subject to 
the assumptions of the '2-line theorem' is really wide. Of course, fixed- 
precision integrating of evolution equations is simpler that integrating 
of arbitrarily chosen equations of mathematical physics. Therefore, 
solving the initial boundary- value problem is divided to two stages. 
First, we find an evolution representation of the equation 8. Then 
we reconstruct the boundary conditions for the auxiliary dependent 
variables if there appeares a necessity to introduce them. 

Example 24. Consider the boundary- value problem 

U^:y,±Uyy = SX]) {2u) , V= {(X,?/),|x| < 1 , | ?/ | < l}, 

u{l,y) = fi{y), u{x,l) = /2(x), u{-l,y) = /ad/), u{x,-l) = fi{x) 

(200) 
for the elliptic (resp., hyperbolic) Liouville equation £^Liou- Introduce 
the additional dependent variable v = Uy. Then the equation Euou has 
the form 

d_ fu\ _ f V 

dy [vj ~ \TUxx ± exp(2M) 
and the conditions on the bound &D are split to the initial and the 
boundary conditions: 

u{x,-l) = fiix), u{l,y) = fi{y), u{-l,y) = f^iy) 

v{x, -1) = v[f2], v{l, y) = dfi{y)/dy, v{-l, y) = df3{y)/dy. 
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Thence, solving boundary- value problem ()200|) is reduced to recon- 
struction of the initial section v{x, —1) by the terminal section u{x, 1) = 
/2(x). In practice, this can be done by using the conjugated gradient 
method. 

A. 3. Evolution of boundary conditions. Recall that any symmetry 
if G sjmS°° of an equation S = {F = 0} is an element of the kernel 
ker ip of the linearization 

T^D .1.. 

restricted onto S. The evolutionary vector field 

d 



3^ = j2da^' 



dui 

commutes with the total derivatives Di and is tangent to the infinite 
prolongation 

£°° = {D„{F) = 0, |a| > 0} C J°"(7r). 

The latter is the projective limit of smooth epimorphisms by definition. 
The manifold S°° that is defined by the closed algebra of smooth func- 
tions is in fact infinite-dimensional, therefore, formally, the field Q^p 
has no flow. Indeed, there is no Cauchy theorem for the initial-value 
problem 

u„ = D^{!f),\a\>0, u{t = 0)=u (201) 

composed by the cardinal set of equations. The integral trajectories 
for problem ()20H) exist provided that {p{x,u,Di{u^)) is a contact sym- 
metry. In particular, the point summetries that are linear with respect 
to the derivatives suit well. From Eq. ()201|) it follows that the func- 
tions ip describe the correlated evolution of the dependent variables and 
their derivatives. This evolution preserves solutions of the equation S 
if V? e kerfi?. 

A. 3.1. Suppose a solution uq of some boundary-value problem Vo = 
{S, V, /o) for an equation S is known. Then, deform the problem 
V = {£, V, f) such that fit = 0) = /o and f{t = 1) = /. Here by f{t) 
we denote the deformation velocity of the boundary condition on dT> 
and by ip{t) the corresponding deformation of the solution. One easily 
checks that ip G keiip- Suppose further that the problem 

^e ker If, <^la© = / (202) 

is soluble at each t G [0, 1). Then, 

U = Uq + / ip{t) dt 

Jo 
is the required solution for the problem V. 
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The condition ip e ker i.p is the decomposition of the deformation ii = 
(/? to a power series in t. The analysis of the deformation equations that 
are coefficiets of the higher powers of t is interesting by itself, see [HI] • 
For example, insolubility of the equation at t^ implies impossibility to 
solve the problem V by using this deformation method. 

Remark 34. The Liouville equation £^Liou has the following peculiar 
feature. Its general solution, see Eq. (jHTj) . that depends on an arbi- 
trary holomorphic function v{z) is always invariant with respect to a 
point symmetry X whose coefficients are related with v{z) by the Abel 
transformation. The curve X(t) in the space sym£^Liou is assigned to 
the deformation / of the boundary condition. This curve is such that 
the solution u{t') is invariant with respect to X{t') for any t' G [0, 1]. 
Therefore, the solution of boundary- value problem ()195|) can be re- 
duced to analysis of the equations that define the curve X(t). 

A. 3. 2. Consider a particular covering strucuture S'^ — > S°° over the 
equation S at hand, that is, suppose that the substitution u = u[v] 
maps solutions v G Sol£^ to solutions u G Sol£^. Several examples are 
well known: the Cole-Hopf substitution 

n = — , (203) 

V 

the Miura transformation m = f^. — f ^, and the formula 

u = -log[4dv-dv/{l-vvy] (204) 

by Liouville ([01]) that relate the Burgers and the heat equations, the 
Korteweg-de Vries equation and the modified Korteweg-de Vries equa- 
tion, and the Liouville and the Cauchy-Riemann equations, respec- 
tively. We recall that all three expressions ()19fi|l are transformed to 
Eq. ()204j) by an appropriate change of variables. Also, suppose the con- 
dition / is fixed. Then there can be several substitutions u[v] that solve 
the boundary-value problem V. Therefore, the reconstruction problem 
for the function v in the whole domain V is incorrect by Hadamars. 
Still, assume that a class of the substitutions f ^— *> m is fixed. Then one 
easily obtains the coordinate expressions for the equations 

which are analogous to Eq. ()202|1 . plus the condition 

B^{u[v]) = f 
defined by the boundary functions / and /o, and the quadrature 



v{l) = f(0)+ / ilj{t)dt. 
Jo 

Thence we construct the solution u of the initial boundary- value prob- 
lem V in the whole domain V by using the exact formula u = u[v] 
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and the solution w(l) of the covering equation £. Still, we see that the 
condition B^(M[f]) = / may not even be defined by an operator with 
directional derivative if the function u depends on the gradient grad v 
explicitly. 

Nawadays, there exist regular algorithmic methods ([ID], see also 
Chapter 4) that allow obtaining and classification of the coverings over 
equations of mathematical physics. These algoriths are already avail- 
able as environments (^72]) for the symbolic transformations software. 
We also recall that the concept of coverings over differential equations is 
closely related with the theory of Backlund transformations and zero- 
curvature representations (see the preceding section) for PDE. The 
latter structures also provide some classes of the substitutions u[v] for 
a given equation £. 

A. 3. 3. Now we consider in more details solving the Dirichlet bound- 
ary-value problem, see Eq. ()195|) . for the Liouville equation £^Liou within 
the class of substitutions ()196b|) by using the methods described above. 
Recall that the harmonic functions v{t) admit the representations P{g\ 
via the Poisson kernel P (see P^) by their boundary value g = v\g.j-,. 
Therefore, the initial boundary- value problem V which is solved by 
using the homotopy / of the solution u such that f{t) are its values on 
dV is reduced to the integral equation 

with respect to the deformation g = iplgj^ oi the bounder value for the 
harmonic function v. Suppose the boundary value for v is obtained, 
then an approximation of w on a lattice in V can be obtained by using 
multiple averaging of the values of v in the neighbouring interpolation 
points (see also Sec. lA. 41 below). 

In Remark lA. 6. 21 on page 11251 the results of a computer experiment 
in application of this method for solving boundary- value problem ()195|) 
are discussed. 

A. 3. 4. Finally, we note that solutions of some other equation S' (in 
general, distinct from the equation S at hand) can be used for construc- 
tion of the initial sections uq that correspond to the boundary values 
/o on dV. Here we assume that the solutions u(t) satisfy a degenerate 
equation £' as t — > 0. 

Example 25. Suppose / =^ — oo such that max/ — min/ < const < oo. 
Then a solution u of boundary-value problem ()195|1 in the disc Bq 
approximated with respect to the norm || • ||co(b1) by the harmonic 
function P[f] as closely as desired. Let 

Au{t) = exp(2M(t)), f{t) = f + logt, < t < 1, 
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then in the same notation we have 



t 
-1 



A(^(t)-2exp(2M(t)V(t) = 0, f{t)=t-\ u{t)=ut^,^+ ip{T)dT, 

where tmin -^ +0 and wt^;^^ = P[/ + logtmin] is the initial approxima- 
tion. Now suppose t is small. We expand the solution uif) by using 
the Hadamars lemma and obtain u(t) = logt-|-P[/] +t'^ -U. If t — > +0, 
then the function U satisfies the homogeneous problem 

A^ = exp(2P[/]), U\g^.=0 
for the Poisson equation. Its solution is 

f/(p)= / exp(2P[/])(g)-G(p,g)dg, 



where G is the Green function ([34 ) for the Laplace operator in the 
unit disc. 

In fact, the method described in Sec. IA.3.4l is based on the simulta- 
neous deformation V of the boundary condition / and the equation S 
itself. Indeed, we have S{0) = £' and £(t) = S ioi t E (0, 1] by con- 
struction. In Sec. IA.4[ we consider an opposite situation: an auxihary 
equation £' is used at any value of the parameter t. Here we assume 
that solutions of S' tend to solutions of the problem V as t increases. 

A. 4. The relaxation method. Complement the mixed boundary- 
value problem, which is a generalization of problem p98|) . by the re- 
laxation term du/dt. Assume that a and (3 are the lower and the upper 
solutions of the initial stationary problem, respectively. Then, extend 
the boundary value / onto &D x R_,_ and fix an initial approximation 
Mo such that a ^ uq ^ /3: hence we obtain 

— Au = —h(u,x), a'Tz^ + hu = f, uix^Q) = Uq{x). (205) 

at an 

Proposition 74 ([ZH])- Let the above assumptions hold. Assume fur- 
ther that h G C^[a,/9]. Then the following two statements are equiva- 
lent: 

(1) the stationary solution Ug of problem ()2()5|1 is unique in [a, (3]; 

(2) the solution Us G [a,(3\ is asymptotically stable such that the 
stability domain is {a,(3\. 

Here we discussed a method based on the replacement of the elliptic 
equation £ in problem ()198j) by the evolution equation £' . This method 
is an alternative to the evolutionary representation method for the 
initial equation £ (see Sec. IA.2|) . 
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A. 5. Deformation of the equation. Now we consider the deforma- 
tions V of the boundary-value problem V{t) that are induced by the 
evolution £ of the equation £{t) = {F(t) = 0} at hand. Here we as- 
sume that the boundary value / remains invariant. Namely, consider 
the homotopy 

S{t) = {l-r{t))-£o + r{t)-S, 

where r(0) = and r(l) = 1. Then the velocity ip of deformation of 
the solution u{t) is subject to the boundary- value problem 

l^(i)(v.) + r'{t) ■ {F - Fo) = 0, iflg^ = 0. 

Therefore, the solution of the problem V is defined by the quadrature 
u = Uq + J^ (pit) dt. We emphasize that this method provides reliable 
solution approximations in computations. 

A. 6. Discussion and practical hints. 

A. 6.1. Consider the deformation {l—r{t))-VQ+r{t)-V of the boundary 
values (see Sec. IA.3|) or of the equation itself (see Sec. IA.5|) . Then the 
smooth step-like homotopy function 

r{t) = exp(-ctg^(7rt/2)), < t < 1, 

is more preferable than the linear motion V = V — Vq. Then, the final 
stage of solving the boundary- value problem V is in fact equivalent 
to the Newton method that starts with the approximation obtained in 
the initial computations. 

A. 6. 2. The comparative analysis of practical computations by using 
the algorithms described in this appendix was carried out in 2002 in 
the dimploma papers by A. V. Punina (Chair of Higher Mathemat- 
ics, ISPU; Master thesis "Comparative analysis of the methods for 
solving the elliptic Liouville equation by using the homotopies") and 
N. P. Cheluhoeva (Chair of Higher Mathematics, ISPU; Master the- 
sis "Integral equations on the border method studied for solving the 
elliptic Liouville equation"). The following assertion is in order. 

• The method of deforming the equation S (see Sec lA.5|l is twice 
more precise with respect to the absolute deviation from known 
exact solutions than the method based on the deformation of 
the boundary function, see Sec. IA.31 although the latter method 
converges 3-5 times faster. 

• The relaxation method is the simplest among the algorithms 
such that the evolutions of t is continuous. The interative 
method is also preferable for solving the Laplace equation that 
appeared in Sec. lA.3l since the use of the Poisson kernel requires 
greater time and the Gauss method cannot be speeded up ow- 
ing to the strongly sparse matrix of the discrete Laplace opera- 
tor A. In Sec. IA.31 we described the reduction of the Dirichlet 
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border-value problem to the one-dimensional problem and the 
integral equation. This method is very sensitive with respect 
to smoothness of the boundary conditions f{t) since the Pois- 
son kernel has a singularity on OBq and the harmonic functions 
are less smooth on OBq than in T> (where they are infinitely 
smooth). We also see that a solution of the Liouville equation, 
when moving along t, can be "attracted" to a nonclassical (of 
type ()197b|l or ()197cj) ) solution which satisfies the same bound- 
ary conditions. This can happen if the time step is suffucuently 
large. 
Summarizing, we conclude that all these methods demonstrated com- 
parable precisions. Therefore, the choice of an algorithm should be 
based on the actual properties of the problem at hand. 

A. 6. 3. Conservation laws for the analysed equation S serve an impor- 
tant instrument for the precision control in computations. Nowadays, 
there exist regular methods (jTUl 113 IMl) of reconstruction of the ex- 
haustive set of conservation laws for the equations subject to the '2-line 
theorem'. These methods are realized in the form of the computer soft- 
ware ([Z2]) for systems of analytic transformations. 

A. 6. 4. Application of the methods of soving boundary-value problems 
described in this appendix does neither neglect nor underesteeme the 
standard check of continuity of the resulting solutions, their Holder 
or the Sobolev space properties an so on. Meanwhile, we hope that 
these practical ideas will compliment the everyday set of instruments 
for numerical analysis of nonlinear equations of mathematical physics. 
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Final remarks 

1. Recently, Demskoi and Startsev (|19J) analysed the correlation be- 
tween the integrals fl and symmetries if for the Liouvillean hyperbolic 
systems. They assigned the operators D that factor symmetries of 
these systems to the product 

to the linearizations i^^i of the integrals Q. These results are reported 
in the note [19 , which is found in this issue. In fact, they generalize 
the statements of Lemma El on page 1391 and Lemma IHUl on page U^ to 
the case of arbitrary integrals f2*, i > 1 (we recall that Q^ ^ T). 

2. The commutative Hamiltonian hierarchy 21 of the local Noether 
symmetries ipk G sym^Toda, where k > 0, was constructed in Part I. 
We identified this hierarchy with the sequence of higher r-component 
analogs for the potential modified Korteweg-de Vries equation and re- 
lated 21 with the bi-Hamiltonian hierarchy 03 for the scalar potential 
Korteweg-de Vries equation, see Eq. ^^ on page EH The degeneracy 
of the matrix coefficients of the higher-order derivtives is an immanent 
feature of these evolution equations. 

A comment is in order. In the fundamental paper j^^J; the integrable 
bi-Hamiltonian hierarchy €. of the Drinfel'd-Sokolov equations was as- 
signed to the hyperbolic Toda system associated with a semisimple Lie 
algebra g or a Kac-Moody algebra g. The hierarchy (L was related with 
the sequence 2) of the multi-component Kortewegde Vries equation's 
analogs. Unlike in the case of the hierarchy 21, the symbols of the 
Drinfel'd-Sokolov equations are always nondegenerate. Therefore, it 
is quite reasonable to check whether there is an interrelation between 
these two systems 2t and €, as well as 03 and S). Also, the necessary 
and sufficient conditions for 21 to be bi-Hamiltonian remain unclear. 

Suppose the hierarchy 21 is assigned to the Toda equation (j^ which 
is associated with a nondegenerate symmetrizable matrix K. We con- 
jecture that the hierarchy 2t is bi-Hamiltonian if and only if K is the 
Cartan matrix of a semisimple Lie algebra. We also conjecture that the 
pair {Ai, A2) of the operators is Hamiltonian in the sence of Definitional 
on page HUl under the same assumptions. 

3. Chapter 3 contains the exposition of the geometric structures for the 
scalar heavenly equation, which is a continuous dispersionless limit of 
the r-component Toda systems as r tends to infinity. We observed that 
local structures for the limit equation are relatively few, and therefore 
a nonlocal setting must be introduced. Meanwhile, local Noether's 
symmetries, conservation laws, and recursion operators for the Toda 
equations themselves were considered in Part I. Therefore it is logical 
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to describe the permutability properties of the diagram 

<^Toda ^ local structures for £^Toda 

) — >oo, 

Uzzz=0^ 

^heav ^ nonlocal structures for £^heav, 

that links together the local geometry of £^Toda and (yet undiscovered) 
nonlocal geometry of £^heav 

4. In Chapter 4, we applied the cohomological schemes developed by 
I. S. Krasil'shchik and constructed one-parametric families of Backlund 
transformations for Eq. ()158|) only, leaving apart the general case ([3]) 
of Backlund transformations for the Toda equations associated with 
the semisimple Lie algebras. Still, from the resulting expressions in jH] 
it is clear that the scaling symmetry is the required generator of defor- 
mations for any r > 1 and any g of rank r. 



The author hopes that the reasonings of the present paper demon- 
strate the profits one obtains by using the invariant coordinate-free 
approach towards the mathematical physics equations. The reader will 
enjoy an excursion to the world of the PDE-related algebraic structures 
in the recent paper jH2], where a natural class of iV-ary generalizations 
for the Lie-algebra structures (in particular, of the symmetry algebra 
sym ^Toda for the Toda equations) was considered. 
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